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ABSTRACT 


This paper discusses the electrical and magnetic design and analysis of a permanent 
magnet generation module for naval applications. Numerous design issues are addressed and 
several issues are raised about the potential improvements a PM generation system can offer. A 
proposed 16 MW PM generation module design is presented along with a detailed design 
methodology. 

Eighty different machines and power conversion modules are sized, designed, and 
analyzed with a final design selected. Specifically, sizing and detailed machine design and 
analysis is performed examining the effects of numerous parameters including number of phases, 
number of poles, magnetic geometry, machine dimensions, and material types. Analytical 
models are developed to study rotor losses caused by stator winding time and space harmonics 
and slot space harmonics. Power electronics and conversion modules to connect the high-speed 
generator to a DC distribution system are designed and analyzed. In depth simulation of the 
eighty complete systems is performed using the software programs MATLAB (Version 12.0, 
Mathworks) and PSIM (Version 6.0, Powersim, Inc.). 

The 16 MW permanent magnet generation module, consisting of the generator and 
associated power electronics, provides an excellent alternative to traditional wound rotor 
synchronous machines. The final design offers significant reductions in both weight and 
volume. Specifically, it is estimated that the PM generation module has a 7x reduction in 
volume and a 10x reduction in weight compared to similarly rated wound rotor systems. These 
reductions can provide flexibility to naval architects since power, weight, and volume are 
integral parts of the design and construction processes. However, further study is necessary to 
verify the PM generation modules thermal, structural, and mechanical performance. 
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Chapter 1 Introduction 


1.1 Purpose 


The purpose of this thesis is to design and analyze a permanent magnet generator and 
power module for naval applications. When deciding whether to implement an electrical 
technology or component onto a naval vessel, the size, weight, and cost are the major factors for 
successful integration. Significant performance improvements and cost reduction of power 
electronics, coupled with the availability and decreasing cost of high energy permanent magnet 
(PM) materials makes PM generators attractive for naval usage. These machines offer numerous 
desirable features, including light weight, small size, simple mechanical construction, easy 
maintenance, good reliability, and high efficiency [1]. 

Before analysis of a generator can begin, it must be properly designed for typical naval 
power requirements. This involves sizing the generator along with designing the associated 
power electronics for connecting the machine to the distribution system. A specific concern 
associated with PM generators is possible inefficiencies and excessive heating; in particular rotor 
losses caused by space and time harmonics during the energy conversion processes. The 
optimum machine design is one that delivers the required power through a matching process 


between the generator and the power electronic converter [2]. 


12 Problem 


The Navy's commitment to develop an integrated electric power system for the next 
generation warships offers the expectation of using the installed generation capacity to power 
ship propulsion, advanced weapons components, and high power combat control systems [3]. As 
these electrical loads increase, it becomes increasingly important to efficiently utilize installed 
power as well as develop smaller, effective power generation systems. Navy ships are extremely 
high performance systems and therefore power and weight considerations are integral parts of 
the design process. 

The life cycle of a navy ship is on average 2-3 times longer than a commercial ship and 
therefore navy ships undergo excessive modernization and upgrades throughout their service life. 


Many of the newer components have significantly higher power requirements than the originals 
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putting a burden on the power generation system. Therefore, the Navy is moving toward 
designing all electric ships with integrated power systems (IPS) and increased power generation 
that can be efficiently managed to meet future demands. 

The effective integration of electrical power in future naval ships requires the 
development of technologies that ensure volume and mass reduction in critical mechanisms. 
Military ships require higher power density components, impose more stringent acoustic and EM 
signature requirements, and subject systems to harsher environments than commercial 
applications [4]. Rotating generators, coupled with prime movers, need to be lighter in weight 
and higher in power density. High-speed PM generators provide a substantial reduction in size 
and weight making them a logical choice for naval applications. Currently, the Navy has not 
designed or built a high power (megawatt) PM generator and therefore the need exists. 

In conjunction with constructing a PM generator, a DC bus architecture is one of the 
preferred schemes for the future [5]. DC power distribution systems can offer a size and weight 
reduction compared to high-power AC systems [3]. With a DC bus distribution, the PM 
generator can be optimized independent of producing 60 Hz frequency as required in the past. 
However, conversion of the high-frequency generator AC output to DC requires power 
electronics. Rectification of AC to DC presents the problem of creating harmonics in the input 
current which are then reflected back onto the generator causing rotor losses. In addition, the 
generator produces space harmonics which also produce losses in the rotor. Therefore, the PM 
generator and power electronics module (PEM) need to de designed and optimized to deliver 


constant DC power while minimizing machine losses. 
13 Background 


1.3.1 History 

The Navy has designed and built electric ships since the early part of the 20 century. 
The original advantages perceived for electric ships, superior performance, reduced manning, 
arrangement flexibility, and fuel efficiency, are still relevant today [6]. In the early part of the 
20" century, diesel-electric submarines, small surface ships, and some battleships and carriers 
had electric propulsion. By the late 1940s, mechanical drive systems became popular because of 


improvements in metallurgy and manufacturing. However, the capability of mechanical 
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transmissions reached its limits in the late 20" century spawning a renewed interest in electric 


drive and integrated power systems for military applications. 


For the last thirty years, the commercial industry has designed and operated integrated 


power systems (IPS). An IPS is a ship architectural paradigm in which the ship's power and 


propulsion are provided by a common electrical distribution system instead of having a separate 


mechanical drive for propulsion. For military applications, an IPS provides numerous benefits: 


Decreased life cycle costs because of increased fuel economy and efficiency (a Navy ship 
with IPS may consume 10-25% less fuel than a similar ship with mechanical drive [7]) 
Increased ship design and arrangement flexibility since the ship is not limited to having a 
long mechanical shaft line 

Reduced system complexity 

Higher degree of modular design using power components 

Broad industrial base for implementing IPS design 

High levels of automation and control 

Increased power available for non-propulsion uses since a Navy ship spends a large 
portion of its time operating at low propulsion levels (approximately 95% of the time) 


Increased stealth, survivability, and payload 


Since an IPS provides power for both the ship loads and propulsion, larger generation capability 


is required. As the Navy builds its new ships with IPS architectures, compact, high-power 


generation systems must be examined to help facilitate implementation of the new designs. 


1.3.2 


Power Generation & Distribution 


Almost all naval core power generators to date are air-cooled 50/60 Hz machines that are 


military derivatives of commercial generators and are therefore size and weight excessive [8]. A 


typical turbine generator system is shown in Figure 1 and Table 1 contains nominal 


characteristics for several current commercial and naval generators. 
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Typical Turbine-Generator System 


Figure 1: Typical Turbine Generator System 


Table 1: Examples of Current Generator Characteristics 


With the advent of high-power, cost-effective power electronics, it is no longer necessary to | 
generate power at 50/60 Hz so generators can be optimized independent of frequency. High- 
speed, power-dense generators become the logical choice for naval purposes. Permanent magnet 
machines are ideal for this high-speed application due to their simple structure and high power 
density [9]. 

Since the generator can now be designed to produce higher frequencies, distribution 
architectures are not limited to being 60 Hz. Either a high frequency AC system or a DC system 
can be designed through the use of power electronics modules, with the DC distribution being 
preferred because of its advantages in size and weight. In this thesis, a high-power DC zonal 


architecture is assumed using solid state converters to generate AC where needed. Each zone is 
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electronically isolated from the other zones with automatic fault detection and reconfiguration to 


provide continuous power during damaged conditions. The PM generator sets and power 


electronic conversion modules serve as the backbone thereby providing a reliable power system 


for navy ships. 


14 Scope 


The scope of this thesis is limited to the PM generator and associated power electronics 


AC-DC conversion module. The following is accomplished: 


Determine the electrical power requirements for a Navy IPS ship in order to properly size 
the generator 
Compare typical wound rotor machine design to a permanent magnet design to determine 
applicability for IPS applications 

Conduct material analysis and selection for the generator design 

Perform initial PM generator detailed design 

Design the power electronics conversion module to perform high-power AC-DC 
conversion 

Conduct detailed analysis of rotor losses of the PM generator, in particular those caused 
by time and space harmonics | 

Perform numerous iterations of machine and power electronics designs to develop 


optimized generation scheme 
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Chapter 2 Power Requirements and Machine Selection 


2.4. Machine & Module Requirements 


On board Navy ships, electricity is used to provide power to virtually all components, 
including mission systems, support systems, combat systems, and communications systems. In 
addition, as ships continue to be upgraded and modernized, more power is needed for newer 
combat systems and weapons components. Most current naval platforms have some form of 
mechanical propulsion system with separate ship service electrical generators supplying the 
ship's power. With an IPS ship, the ship's generators provide power for propulsion and the 
ship's service loads, and through proper utilization, power is efficiently managed. 

To properly size the PM generator, the power requirements must be identified and 
therefore a typical load list is developed for an IPS naval ship and is included in Appendix A. 
The overall power requirement for the generator is 16 MW. Since size and weight are important 
factors and the generator can be optimized independent of frequency, high-speed operation and 
maximum power-density are desired. Therefore, the highest possible speed is selected while 
ensuring the PM generator is compatible with both gas turbines and steam turbines. 

Traditionally, gas turbines run at much higher speeds than steam turbines causing the 
steam turbines to be more limiting. From information collected from the Elliot Turbomachinery 
Company, Inc., 16,000 RPM is approximately the highest speed steam turbine that can 
reasonably be constructed at the megawatt power level [13]. Therefore, to provide a degree of 
conservatism, 13,000 RPM is selected for the nominal design speed for the PM generator. 

The power electronics module (PEM) converts the AC voltage from the generator to 
700 VDC. Overall, the PEM and generator must be designed so that losses suffered by the 
permanent magnets on the generator rotor are minimal. Table 2 lists the general requirements for 


the entire system. 


22 Machine Selection 


Military ships require high power density components and improved acoustic and 
electromagnetic signature requirements while subjecting systems to harsh environments [14]. It 
is therefore important to ensure the power generation system is capable and efficient. 
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Table 2: General Machine & Module Requirements 


2.2.1 Permanent Magnet versus Wound Rotor 


Reducing the size and weight of ship's turbine generator sets offers significant 
advantages to naval architects. Replacing older generators with lightweight ones could make it 
possible to decrease the size of some generator sets by as much as 5096 [15]. PM generators 
therefore become an attractive alternative compared to wound rotor machines because of the 
availability and decreasing cost of high energy PM materials along with improved power 
electronics. 

A wound rotor generator normally consists of armature windings on a stationary stator 
frame with field windings on an inner rotor. The rotor is turned by a prime mover, usually a gas 
or steam turbine, and current is supplied to the field windings through brushes or a brushless 
exciter. As the current-carrying field windings rotate past the stator windings, current is 
produced in the stator windings through Faraday's Law. An example of a wound rotor machine 


is shown in Figure 2 (8] and Figure 3 [16]. 


Figure 2: Example of Wound Rotor Generator 
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Figure 3: Cross Section of Wound Rotor Generator 


Wound rotor generators have been the backbone of power generation for the U. S. Navy because 
they are a proven technology that is battle tested. They offer steady voltage regulation using the 
field windings, a large air gap for producing the rotor magnetic flux, low fault currents, and high 
power capabilities. | 

However, the machines tend to be complex, weight excessive, and they require field 
windings which limit design alternatives. There are several drivers which cause these problems 
to occur. First, to generate the necessary magnetic flux levels, wound rotor generators have large 
pole pitches to support the required field windings. These pole pitch windings in turn require 
larger end turns and thick back iron to support the magnetic flux, both of which contribute to 
increased size. Second, because of the winding losses in the rotor, large cooling systems can be 
required thus increasing the number of support components. 

High speed generators offer a reduction in machine size and weight because as a 
machine's speed increases, its size decreases for a given output power. The PM generator is 
ideal for high-speed applications because of its simple structure and high power density [9]. In a 
PM generator, the rotor field windings are replaced by permanent magnets which do not require 
additional excitation. As the permanent magnets are rotated by the prime mover, current is 


produced in the stator windings. An example of a PM generator is shown in Figure 4 [8]. 
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Figure 4: Example of PM Generator 


PM generators offer several advantages: they have no rotor windings so they are less 
complicated; they have high efficiencies; the gap field flux is not dependent on large pole pitches 
so the machine requires less back iron and can have a greater number of smaller poles; and they 
usually require smaller and fewer support systems. Assuming the same flux density and 
circumferential arc, doubling the number of poles produces the same radial flux but requires half 


the stator core thickness, as shown in Figure 5. 


Half Stator Core Thickness.- 
“Same Radial Flux — 


Figure 5: Flux vs. Number of Poles 
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However, PM generators also have some disadvantages. They do not possess field 
excitation control and therefore voltage regulation can be problematic. This can be corrected by 
using external voltage control such as large capacitor banks or power electronics, as well as 
choosing the turns on the stator winding properly to produce the anticipated required nominal 
voltage. Additionally, since the permanent magnet fields cannot be turned off, there exists the 
risk of excessive currents in the event of an internal fault. This problem can also be solved 
through the design of the turbine governor and controller or dynamic braking. Overall, the 
advantages of the PM generator over the traditional wound rotor generator make it a better 
alternative for high-speed navy applications. A summary comparing the different designs is 


given in Table 3. 


Table 3: Comparison of Wound Rotor and PM Generators 


Steady voltage regulation with Weight excessive 
field windings 


High power capabilities Large size 


Wound Rotor Rotor windings & associated 
. Large air gap for flux losses 


Low fault currents 
Proven, robust design T 
Lack of inherent voltage 
Less complicated ' 8 
regulation 


Reduced size and weight Potential fault currents 
Permanent Magnet High efficienc 


No excitation supply or field 
windings 


High speed applicabilit 


2.2.2 Type of Permanent Magnet Machine 

There are numerous layout possibilities for permanent magnet machines and only the 
most common are discussed here. These include radial flux inner rotor, radial flux outer rotor, 
and axial flux designs. In most PM machines, flux crosses from the rotor to the stator in the 
radial direction [17]. The first type, the radial flux inner rotor design, is the closest configuration 


to the classical AC synchronous generator. An example of this design is shown in Figure 6 [17]. 
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Permanent 
Magnets 


Figure 6: Example of Inner Rotor PM Machine 


In this type of machine, the windings are placed on the stator, either in slots or in a 
slotless ring, and the magnets are surface mounted on the rotor or buried in the rotor. Buried 
magnet designs often result in rotors that are larger than equivalent surface-magnet machines 
with high-energy magnets [18]. Buried magnet machines can also have significant structural 
issues in high-power applications [19]. When the magnets are surface mourited and the machine 
is operated at high speed, the magnets are often secured with a retaining device made of either 
alloy steel or carbon-fiber. Overall, the inner rotor machine possesses high torque/power 
capability and good heat conduction and cooling properties making it ideal for high-speed, 
higher-power applications [18]. i i 

The radial flux outer rotor machines are commonly used in hard disk drives, small 
computer ventilation fans, and some blowers. This type of design is very efficient, low-cost, 
easy to manufacture, and applicable for low-power applications [18]. Itis the opposite of the . 
inner rotor because the stator in on the inside with the rotor and magnets on the outside. A cross 


section of an outer rotor machine is shown in Figure 7 [17]. 
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Permanent 
Magnet 


Figure 7: Example of Outer Rotor PM Machine 


The axial flux machine is significantly different than the previous two because flux flows 
in the axial direction vice radial direction and the windings are oriented radially vice axially (see 


Figure 8 for an example diagram [20]). 


intermediate rotor disc 


Figure 8: Example of Axial Flux PM Machine 


The main advantages of this design are their low cost, flat shape, and smooth rotation. However, 
if axial-flux machines are operated at high speeds (above 1000 RPM), eddy-current losses and 
heating can become excessive [18]. Also, stator construction is difficult because it must be 


laminated circumferentially. An example of this design is the turntable for a record player. 
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Overall, because of its inherent advantages in heat removal and cooling, the abundance of 
manufacturing capabilities, and its high-power, high-speed applicability, the radial flux inner 


rotor with surface mounted magnets is selected for the 16 MW PM generator design. 
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Chapter 3 Material Selection and Machine Initial Design 


3.1 Material Selection 


One of the key considerations during the electromagnetic, structural, and thermal design of 
a permanent-magnet machine is the selection of the magnet, stator, and rotor materials [21]. 
Machine output, heat rise, weight, and cost are a few of the characteristics which are directly 


influenced by selection of the machine materials [22]. 


3.1.1 Permanent Magnets 

The size and performance of high-speed PM generators depend on the permanent magnet 
material properties [9]. The magnets must be selected to provide the necessary air gap magnetic 
field and ample coercive force to compensate for possible damaging effects while minimizing 
the volume of material because of cost and weight considerations [23]. 

Ferromagnetic materials are the most common substances used in the construction of 
machines and their properties are normally described using B-H curves and hysteresis loops. 
These curves represent an average material characteristic that reflects the non-linear property of 
the permeability of the material but ignores the multi-valued properties [17]. An example of a B- ` 


H curve is shown in Figure 9. 


. B (Tesla) 


Load Line 


-Initial 
‘Magnetization: 


Recoil Line H (kAlm) 


Br: Remnant Flux Density Hc: Coercivity 
Bm: Magnet Flux Density Hk: Limiting/knee magnetizing 
force 


Figure 9: Example of B-H Curve 
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Several basic magnetic properties are of critical importance for the permanent magnets in a PM 
machine: 

e Remnant Flux Density (B,): It is the value of the flux density remaining after 
magnetization and it directly influences the air gap flux and magnet sizes. 

e Coercivity (Ho): It is the value of magnetizing field needed to reduce the flux density in 
the magnet to zero and it gives a first order estimate of a magnet's resistance to 
demagnetization. 

e Energy Product (BHmax): It is the maximum energy product of the magnet and it is 
inversely proportional to the total magnet volume required. | 

è Recoil Permeability (uwe): It is the gradient of the B-H curve and it gives the magnet's 
ability to return to its initial magnetization after subjected to damaging forces. If the 
magnet goes below Hk, then it will recoil along a lower line resulting in a lower magnet 
flux density. | 

e Load Line: It is a line drawn from the origin to the magnet operating point on the 

hysteresis curve (Bm). The magnitude of the slope of the load line is the permeance 
coefficient. 
Permanent magnet materials come in many varieties and the four most common types for 
machine applications are Alnico, Ferrites, SmCo material, and NdFeB material. Table 4 and 


Figure 10 show the characteristics and typical B-H curves for these materials [18]. 


Table 4: Magnet Material Properties 


Remanence (B,) 0.6 — 1.3 0.35—0.43 | 0.7— 1.05 1.0 - 1.3 
Coercivity (H.) 40 — 130 180—400 | 800—1500 | 800-1900 


Recoil ud | 19-17 | 1.05-1.15 | 1.02-1.07 | 1.04-1.1 
TEC, 


0 
Energy Product (BHmax) 20 — 100 24 — 36 140 — 220 180 — 320 
500 — 550 250-350 | 100-200 


B: Temperature — | ec | 001t0-002 | -02 0.05 — | -0.08 to -0.15 | 
Coefficient 


The rare-earth magnets, SmCo and NdFeB, have become more popular for high performance 


applications because of their greater power density, high coercivity, high flux densities, and 


linearity of the demagnetization curves [24]. 
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= Alnico = SmCo NdFeB 


Figure 10: Typical Magnet B-H Curves 


Between the two rare-earth permanent magnets, NdFeB is preferred because it is cheaper and 
more readily available. It does possess some adverse characteristics such as moderate corrosion 
and lower resistance to temperature effects, but these can be controlled using surface treatments 
and proper cooling [25]. Therefore, NdFeB magnets are selected for use in the PM generator 


with the conservatively assumed values listed in Table 5. 


Table 5: Selected Magnet Properties 


— Remanence (B) 
Coercivity (H.) 


| Recoil Permeability (pe) | Å 
| Maximum Temperature | — C — | 
[Resistivity | — um | 
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3.1.2 Stator and Rotor Material 

The type of material chosen for the stator and rotor is important because it impacts the 
machine losses and efficiency. The rotor is usually built from the same material as the stator for 
ease of construction but is can be made of any economical steel provided it is strong enough for 
the given function [18]. No one material is optimum for every application and the normal 
criteria for selection are cost, permeability, core losses, and saturation flux. It is important that 
the material act as a flux guide and absorb the minimum amount of magnetomotive force (MMF) 
so that the flux is concentrated in the air gap. In addition, the material should minimize core 
losses including hysteresis and eddy current losses. 

High-quality, non-oriented, electrical grade lamination steels are typically used in most 
machines because the laminations help minimize losses. The four main materials are low carbon 
steels, silicon (Si) steels, nickel (Ni) alloy steels, and cobalt (Co) alloy steels. Low carbon steels 
are the lowest cost and are used in high volume applications where high core losses are 
acceptable. Silicon steels usually have 3% silicon which increases the resistivity to reduce eddy 
current losses. They are selected and specified based on core loss, with each grade (M19, M27, 
M36, and M43) having higher core losses and lower cost [22]. The lamination thickness is a 
tradeoff between cost and performance and the most common sizes are 0.014 in, 0.0185 in, and 
0.025 in (29 gauge, 26 gauge, and 24 gauge). 

Nickel alloys are either 49% or 80% nickel and they have lower losses than the silicon 
steel but are much more expensive. In addition, they require careful handling and not suited for 
high flux density environments (above 0.8 T) because of saturation. The cobalt alloys are only 
used in extremely high-performance situations such as military aircraft and space applications 
because of the high cost. Table 6 summarizes the different stator materials and the M19, 29- 
gauge electrical silicon steel is selected for the PM generator because it is economical, its thin 


laminations minimize losses, and it has a saturation flux density of about 1.8 T [2], [18], [22]. 
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Low Carbon | Fair E Good | 
Steel 
| Good | 


Good 


32 Machine Design Parameters 


3.2.1 Stator Mechanical Design 


The stator is an important part of the machine because it serves as the main structural 
component, it provides the housing for the armature windings, and it completes the flux path for 
the magnetic circuit. The main consideration in the mechanical design of the stator is whether to 
make it slotted or slotless. A slotless stator has the armature windings located in the air gap of 
the machine as shown in Figure 11 [19]. 


Rotor Core 
Stator Winding (Shaft) 


Figure 11: Slotless Stator Design 


One of the advantages of the slotless construction is unique winding layouts are possible to meet 
specific performance goals. Another advantage is that the space available for the armature 


windings increases by a factor of about two since there are no stator teeth. This produces lower 
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conductor losses since less current flows in each winding. The flux density is reduced, however, 
because the effective air gap is much larger since the windings are in the air gap. Overall, there 
exists a higher electrical loading and a lower magnetic loading. 

One disadvantage of the slotless design is there are no good conduction paths to remove 
the heat generated from the windings. This reduces the allowable current density in the windings 
and lowers the power output. Another disadvantage is that the windings are directly exposed to 
the rotating flux which raises the possibility of additional eddy-current loss in the conductors and 
further losses due to circulating currents in the windings [18]. Overall, the performance of a 
slotless stator is almost always lower than that of an equivalent slotted stator design and 
therefore slotless stators do not appear often in high-power applications [17]. 

Slotted stators are the traditional stator design and consist of openings around the stator 
for the armature windings as shown in Figure 12 [26]. The openings provide rigid housings for 


the conductors and associated insulation. 


Back iron 


Figure 12: Slotted Stator Design 
Stator slots vary in size and shape with the most common configurations being rectangular or 
trapezoidal. In this paper, the slots are assumed to be approximately rectangular as shown in 


Figure 13 and contain form-wound windings so that the depression width is the same as the slot 


top width. 
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Web - slot bottom width 
^wst - slot top width: 

T chs -slotheight. — 
"wd - depression width 
hd: «depression height’ 


+Wd 


airgap 


Figure 13: Stator Slot Geometry | 
Slotting is used because is provides a good way to achieve a narrow air gap length while keeping 
the winding conductors close to the magnets to maximize the flux linkage. The slots also greatly 
increase the surface contact area between the windings and stator steel providing a path of low 
thermal resistance for good heat conduction which is important for keeping the windings and 
magnets cool. 

The resulting narrow air gap from the slots makes the permeance greater and therefore 
the air gap flux density greater producing a more powerful machine. In addition, the depression 
in the slot tops help control parasitic losses in the rotor by improving the uniformity of the air 
gap field. The limits of the size of the slots are twofold: the magnetic teeth must be able to carry 
the air-gap flux without saturating and the slots must be large enough to support the necessary 
current density in the windings. Typical limits for stator current density are shown in Table 7 


and in this páper it is assumed that the limit on current density (J) is 3000 A/cm? [18]. 


Table 7: Stator Current Densities 
UNUS ho TE 
450 — 550 


Fan Cooled 800 — 1200 
Liquid Cooled 2300 — 3200 


The disadvantages of the slots are that cogging torque may be a problem and it can be 
costly to insert the windings if proper construction techniques are not used. Overall, however, 
slotted designs are preferred in high-power applications and therefore a slotted stator is selected 


for the 16 MW PM generator. 
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The number of slots in the machine is usually a function of the number of phases and 
windings and can vary based on the application. The initial design of the generator assumes a 
three-phase machine but the number of phases will be examined to determine an optimum 
design. In order to allow for high-power operation and the possibility of a high number of 


phases, 36 slots is chosen for the initial generator design. 


3.2.2 Rotor Mechanical Design 

For high-speed applications, the rotor aspect ratio, defined as length-to-diameter (L/D), is 
a critical parameter. If it is relatively low, then the rotor has high stiffness and good dynamics 
but a large diameter which increases the weight and makes magnet retention extremely difficult. 
Additionally, the centrifugal force on the surface-mounted magnets is directly proportional to the 
rotor diameter so the rotor radial size must not be excessive. 

Permanent magnet machines offer flexibility in selecting pole sizes which allows for 
smaller diameters. They are therefore ideal for high-speed applications because they can have 
higher L/D ratios. This is because they do not have rotor field windings which have end turns 
necessitating big pole pitches and large diameters. A normal L/D ratio for a wound rotor 
machine is 0.5 — 1.0 compared to 1 — 3 for a PM machine [27]. Staying close to these ranges 
usually ptovides a first order estimate of satisfactory machine dynamic performance and 
acceptably low vibrations or oscillations. | 

The rotor radius and the rotational speed also determine the tip speed of the machine 


which is the surface velocity of the rotor (as defined by Eqn 3-1). 


Vtip = Rm 


where or = angular speed (rad/sec) 
R = rotor radius (m) 
Eqn 3-1 
For most rotating machines, the upper limit on tip speed is between 100 — 250 m/s depending on 
the design. For surface magnet PM machines, retaining sleeves are sometimes used to help keep 
the magnets in place and allow for higher speeds. These sleeves can be constructed from alloy 
steel, carbon fiber, or other materials. The metal sleeves usually provide increased mechanical 


performance but have eddy current losses. 
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The carbon fiber and graphite composite sleeves have high strength-to-weight ratios 
which produce a thin sleeve and the sleeve's lower conductivities yield reduced eddy current 
losses. However, the carbon fiber and graphite composite sleeves have lower temperature 
ratings and lower thermal conductivities making heat removal and increased cooling for the 
magnets and sleeve important issues [28]. Overall, the use of a retaining sleeve is necessary for 
the 16 MW generator since it is operating at high-speed and this allows the maximum tip speed 
limit at the rotor surface to be 200 m/s. The actual material for the retaining sleeve is examined 


later when detailed rotor loss analysis is performed. 


3.2.3 Number of Poles and Magnet Pole Design 

The optimum number of poles is a complex function depending on a number of factors 
including the magnet material, the speed of rotation, the desired output frequency, and the 
mechanical assembly of the rotor. An even number of poles is always used because this provides 
a balanced rotational design. As the number of poles increases, the individual pole pitch goes 
down which reduces the amount of stator back iron needed to support the magnetic flux. In 
addition, for a given power/torque, as the pole number rises, the required magnet volume 
decreases. 

Assuming a constant mechanical rotation speed, the generated electrical frequency is 
proportional to the number of poles as shown in Eqn 3-2. 

N-(2p) = 120-f 


where N = speed (RPM) 
p = number of pole pairs 
f = electrical frequency (Hz) 


Eqn 3-2 
If a PM generator is going to be the source for a DC bus through a rectifier system, a high pole 
number is desirable because as the electrical frequency increases, support components such as 
filter capacitors and inductors can be much smaller. Therefore, for a given rotational speed, one 
cheap and efficient solution is to have a higher number of pole pairs and frequency [27]. 
However, as the frequency increases, higher stator losses result because core losses are 
proportional to frequency squared. In addition, as the pole number gets larger, the number of 
slots per pole per phase decreases and can cause the voltage waveforms to become less 


sinusoidal so all factors must be considered. 
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The pole arc of the magnets can also be varied. Magnets seldom span the full pole pitch 
because the flux at the transition between north and south poles leaks between poles without 
linking the coils in the stator. The gaps between the poles usually contain non-magnet pieces, 
such as soft-iron, so that no flux crosses over the air gap between magnets. A full pole arc is Ome 
= 180^ E and produces a full voltage waveform but has increased harmonic content. As the pole 
arc is reduced (up to 20 — 30 96) and those areas are filled in with soft-iron pieces, the resulting 
flux waveform is more sinusoidal and has fewer harmonics and therefore lower rotor losses [29]. 

The magnet poles are sometimes skewed to reduce cogging torque and smooth out 
variations in air gap reluctance, flux, and voltage waveforms. Skewing of the magnets occurs 
axially along the length of the rotor to provide a constant rotational torque and prevent pole 


pieces from exactly lining up with stator teeth. A skew factor is used to account for this effect 


and is shown in Eqn 3-3. 


kan = 


where @,= skew angle, radE 
= n = harmonic number 


sin(n8,) 
0, 
2 
| Eqn 3-3 
As the pole number is increased, the stator conductors-per-pole decreases so that the per- 
unit inductance and synchronous reactance decreases with higher pole number. This can 
sometimes result in improved performance of the machine since the reactance is lower. Overall, 


the initial 16MW generator has 6 poles but this is examined later to determine an optimal design. 


3.24 Magnetic Dimensions 

The primary magnetic dimensions that affect a PM machine are the air gap and the 
magnet height. These two parameters play a major role in determining the air gap magnetic 
field, the air gap flux density, and the induced voltage in the machine. To a first order 
approximation, the air-gap flux density (Bg) can be represented by Eqn 3-4 [30]. 

The radial air gap is usually made as small as possible to maximize the air gap flux 
density, minimize the flux leakage, and produce a lower reluctance value since the air gap 
constitutes the largest part of the machine permeance/reluctance. However, the use of rare-earth 
permanent magnets (NdFeB or SmCo) with their higher flux density and coercive force permit 


some flexibility in the size of the air gap. 
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where ha = magnet height (mm) 
g = air gap (mm) 


B, magnet remnant flux density (T) 
Eqn 3-4 

Once the permanent magnet material is selected, the desired air gap flux density and 
induced voltage help determine the magnet height needed. If the magnet height is too large, the 
air gap flux density might be significant enough to cause the stator core material to saturate ` 
which reduces machine performance. The goal is to use the minimal amount of magnet material 
to achieve the desired effect because this reduces the size and weight of the machine and 
decreases the magnet material cost. Also, losses in the magnets can be reduced by using smaller 
magnets. In order to provide uniform magnetic fields, the magnet height is usually larger than 


the air gap by a factor of 5 — 10. 


3.2.5 Number of Phases 

In general, the number of phases affects a machine's power, current, and voltage ratings 
as shown in Eqn 3-5. If the power is fixed, then as the number of phases increases, the phase 
voltage and/or current decreases, assuming the total number of turns is constant. 


|P + jQ] =q- VI 
where P = real power (W) 
Q = reactive power (VAR) 
q = number of phases 
V = RMS phase voltage (V) 
I = RMS current (A) 
Eqn 3-5 
Most motors and generators are three-phase machines because it is the industry standard, it is the 
most common form of power, and it is the lowest number of phases that produces balanced 
torque with out pulsations in rotating machines. However, higher utilizations in generators can 
be achieved with higher phase numbers especially if the generator is connected through power 
electronics to a DC bus distribution. This is because the higher number of phases produces 


lower ripple in the DC bus voltage. 
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However, the AC line current harmonics are more substantial in increased phase 
machines because the triple-n harmonics are higher order as the phase number increases. For 
example, a 3-phase machine suppresses harmonics of order 3n, a 5-phase machine eliminates 
order 5n, and a 7-phase machine removes order 7n. T herefore, in higher phase machines, a 
greater number of large harmonics result in the AC line current. Also, as the number of phases 
increases, the phase inductances and reactances change since there are a greater number of 
windings influencing each other. 

Most machines are usually designed with the phases balanced meaning that they have an 
evenly-spaced phase distribution around the stator of the machine. This produces voltage 
waveforms that are identical in shape from phase to phase but differ by a phase offset angle. In 
order to initially size the PM generator, it is assumed to have three phases but this will be 
optimized later in conjunction with the power electronics module. It is also assumed that the 


- 


phases are always balanced. 


3.2.6 Slots per Pole per Phase 

The number of slots per pole per phase (m) is an extremely important design parameter 
when considering generator design and it is calculated using Eqn 3-6. It is used to help 
determine the relationship and interactions between the rotor poles and the stator windings as 
well as shape the generated back voltage of the machine. When m is an integer, the machine is 


an integral slot machine and when m has a fractional part, it is a fractional slot machine. 


Ns 


m= 
2-p-q 

where N, = number of slots 

p = pole pairs 

q = number of phases 

| Eqn 3-6 
In an integral slot machine, the back EMFs of all of the coils making up a phase winding 

are in phase with each other and add up so that the final voltage amplitude is the direct sum of 
the individual coil voltages. In a fractional slot machine, the back EMF of all of the coils are not 
in phase so the net voltage has a different shape than the individual winding voltages. Varying 
the number of slots/pole/phase is one method used to produce a more sinusoidal voltage 


waveform and reduce the harmonics generated by the machine. 


38 


| | 


3.2.7 Stator Windings 

The stator windings are the location where the generator voltage is induced due to the 
time varying magnetic flux caused by the permanent magnets on the rotor. In a slotted machine, 
the winding arrangement is used to help shape the back voltage to produce a more sinusoidal 
waveform. The windings can be distributed by three methods: pitch, skew, or 
breadth/distribution. 

The pitch of a winding (a) refers to the angular displacement between the sides of a coil, 
usually expressed in electrical degrees or radians. When the individual coil pitch differs from 
180° E, the winding is said to be short-pitched or fractional-pitched. This causes angular 
segments where the back voltage is zero because the flux linkage is constant and can help 
produce a sinusoidal waveform when multiple coils are connected. It also has the advantage of 
lowering the coil resistance and making the stator end windings more manageable. 

Windings in the stator can also be skewed axially along the length of the machine. This 
requires the stator slots to be more intricately designed which complicates the mechanical 
construction of large machines. Therefore, since the generator being designed is a large, high- 
power machine, skewing of the stator windings is not used but skewing of the rotor is employed. 

The breadth of a stator winding results from the coils occupying a distribution or range of 
slots within a phase belt. A stator winding normally consists of several coils each separated by 
an electrical angle y. The distribution of the coils causes each to link the rotor flux slightly out of 
phase with each other so when they are added together, they produce a more sinusoidal 
waveform. | 

Within each stator slot, there are geometric size constraints which determine how many 
conductors can be placed in a slot. In smaller machines, coils are composed of round insulated 
wires that are placed in the stator slot along with insulation material. A slot fill factor (As) is used 
to determine how much of the slot cross-sectional area is occupied by winding material as shown 
in Eqn 3-7. 

WindingArea 
5' TotalSlotArea 
Eqn 3-7 
In larger machines, form-wound windings are used for ease of construction and for better 


performance. A sketch of what a form-wound winding looks like is shown in Figure 14. 
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Typically, machines contain two coil sides per slot making the winding a double-layer design 
[17], [18], [19]. Overall, slot fill factors vary in value from 0.30 — 0.70, depending on the 


number and size of the conductors in the slots as well as the amount of labor utilized. In this 


paper, a slot fill factor of 0.50 is assumed. 
Insulation Material 


Conductors 


Cooling : 
Channels EE Lamination 
T= Material 


Figure 14: Example of Form-Wound Winding 


In conductors that carry high-frequency currents, skin effect can become an issue and 
affect the operation of the machine. Skin effect is caused by eddy currents in the windings 
themselves due to the changing magnetic field. These eddy currents force the current flowing in | 
the conductor to crowd to the outer edges of the conductor. This in turn causes the current to 
flow through a smaller cross-sectional area and increase the resistance of the conductor. 
However, the generator under design is expected to operate at less than 2 kHz and for 
frequencies below 12 kHz, RAc/Rpc « 1.01 so skin effect can be neglected [32]. 

Within a phase, stator windings can be connected in wye or delta patterns as well as 
series or parallel. Almost all machines use series, wye-connected windings because they provide 
the safest alternative. This is because in a delta or parallel connection, the back EMFs can 
produce circulating currents which can result in addition losses, heating, or damage. Therefore, 


wye series connected windings are selected for use in the designs in this paper. 
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3.3.1 Basic Model 


Since the machine is assumed to balanced, parameters can be determined on a per-phase 
basis and then applied to all of the phases. Each phase of the machine can therefore be modeled 


as shown in Figure 15. 


Figure 15: Per Phase Model 


The armature resistance (R,) is the resistance of the windings of the machine and it is usually 
relatively small. The synchronous inductance (Ls) of the machine comes from the inductance of 
the windings and is composed of the air gap inductance, the slot leakage inductance, and the end- 
turn inductance. The back voltage (E;) is produced through the flux linkage in the windings 
from the rotating magnetic field in the machine. Lastly, V, is the terminal voltage and is found 


using basic circuit analysis once the other parameters are known. 


3.3.2 Winding Resistances 


The stator coils in the machine are made of copper and therefore have some resistance to 
the current flow. This resistance of the copper phase windings is calculated using Eqn 3-8. 
R= 
G.A 
where |= length of conductor 
o = winding conductivity 
A = winding cross-sectional area 
Egn 3-8 
The length of the conductor comes from the windings traveling twice the length of the machine 


and twice around the end turns of the machine. It is assumed that the end turns follow roughly a 
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circular path from one side of the machine to the other where the radius of the circle is the 
distance to one half the stator slot height. The cross-sectional area of the conductor is obtained 
from the slot area and slot fill factor as shown in Eqn 3-9, assuming form-wound windings. 
Aac = = 
2N, 
where A,- slot area 
N, - turns per coil 


Eqn 3-9 


3.3.3 Winding & Magnet Factors 

As discussed in section 3.2.7, windings are normally not full-pitched or concentrated but 
rather are short-pitched and have breadth associated with them. To account for these effects, a 
winding factor (ky) is utilized which is the ratio of flux linked by an actual winding to the flux 
linked by a full-pitch, concentrated winding having the same number of turns. The winding 
factor is the product of a pitch factor (kp) and a breadth/distribution factor (ky) as shown in Eqn 


3-10. 


Eqn 3-10 
The pitch factor accounts for the windings spanning a electrical degrees vice spanning a 


full 180" E as shown in Figure 16 [26]. 


Figure 16: Short-Pitch Coil 
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The pitch factor is the ratio of the flux produced by a short-pitch coil to the flux produced by a 
full-pitch coil. Short-pitching is an important means for eliminating harmonics and improving 
the power quality of the machine. The pitch factor can be derived with the final result shown in 


Eqn 3-11. 


where n = harmonic number 
Eqn 3-11 
The breadth factor explains the effect of the windings occupying a distribution or range 


of slots within a phase belt. A phase winding normally consists of numerous coils connected 


together linking flux slightly out of phase with each other as shown in Figure 17 [26]. 


Figure 17: Winding Breadth 


The breadth factor can be derived either magnetically or geometrically to obtain Eqn 3-12. 


n nm y : 
sin 2 where n = harmonic number 
m - slots per pole per phase 


Kon Ë ny = coil electrical angle 
m so 21) Y= g 


Eqn 3-12 
In addition to estimating different winding effects, the geometry of the magnetic air gap 
must be represented. Field methods are utilized along with vector potential analysis to develop 


expressions that account for different magnetic gap geometries. Reference [19] contains detailed 
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A 


derivations of the magnetic gap factor (kgn) for several magnet and slot configurations. The 


equation for the slotted stator, surface magnet configuration is shown in Eqn 3-13. 


RAP 
E I np Je _ np) å np Ry ú v 


gn 2np 2np | np + I np ~ 1 
RR 
where R, = outer magnetic boundary R, = outer boundary of magnet 
R, = inner magnetic boundary R, = inner boundary of magnet 


Eqn 3-13 


3.3.4 Flux and Voltage 

The primary significance of the magnetic flux linkage in a machine is that it induces 
voltage across a winding whenever the flux varies with time as explained through Faraday’s 
Law. The first step in the process is to determine the air gap flux density. The flux from the 
magnet poles crosses the air gap to the stator windings but some flux leaks along the way and 
this is accounted for using a leakage factor (Ki ~ 0.95 for surface magnets). In addition, the flux 
path is normally dominated by the air gap reluctance since the reluctance of the stator steel is 
much less than that in the air gap. However, a reluctance factor (K, ~ 1.05 for surface magnets) 
is used to compensate for the small effects of the steel reluctance on the air gap flux. 

The presence of the slots in the stator also affects the air gap flux density because of the 
difference in permeance caused by the slots. The flux crossing the air gap in a slot region travels 
farther before reaching the highly permeable stator back iron. Carter's coefficient (Ke) is used to 
account for this effect [17]. The air gap flux density is also affected by the magnet geometry in 
the air gap as previously described by Eqn 3-13. Since the magnet poles rotate north/south, the 
air gap flux density shape can be approximated as shown in Figure 18. This can be represented 
as a Fourier series using only odd components because of half-wave symmetry as shown in Eqn 


3-18. Overall, the air gap flux density is calculated using Eqn 3-14 through Eqn 3-18. 


44 


enerator Top Maenetic Generator Magnet Generator DIY 


Figure 18: Air Gap Flux Density 


-1 


K,=|1- p where w, = average slot width 
=| s J w, = tooth width 
Ws Ws x.=W.+W 
s s t 
Be= Keg 


where g, - effective air gap 


hm 


PC = 
ge C, 


where PC = permeance coefficient 
C = flux concentration factor (A m/Ag) 


Kr Cs 
By = ‘By where Hrec = recoil permeability 
De Kj B, - remnant flux density 

PC 


45 


AC Generator Top Magnetic Generator Magnet Generator DIY 


Eqn 3-14 


Eqn 3-15 


Eqn 3-16 


Eqn 3-17 


B(0) E y B,-sin(np6) 


n=1 4 HPD: nz 
where B = —.B .k. .sin| —— |-sin| — 
n odd n ax 8 Bn 2 2 


Om = magnet physical angle 
n = harmonic number 
Eqn 3-18 
Once the flux density is known, the flux must be calculated. Given a machine uniform in 
the axial direction, the flux linked by a single, full-pitched coil which spans an angle from O to 
n/p is represented by Eqn 3-19. Assuming Baux is sinusoidally distributed, the peak flux for this 
ideal coil is given by Eqn 3-20. 


x 


p 
p | Brux Rs La 48 where B flux = radial flux through coil 
0 


Eqn 3-19 


Eqn 3-20 
Given that there are N, coils in a stator phase winding and including all the real winding effects, 
the total flux linkage is shown in Eqn 3-21. Through Faraday's Law, the back EMF for the 
machine is given by Eqn 3-22. | 


Mo) = y As sin(np8) 


n=1 where 4 = 2 RS Toc Ny By kn Kan 
n odd n I pP 


Eqn 3-21 
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E, = » Vu: sin (np 0) 
Í S where v, = on = Ay 
Eqn 3-22 

With permanent magnet excitation, the field cannot be controlled as in a wound-rotor 
machine so the number of turns in the stator phase windings must be chosen so that the machine 
EMF is close to the nominal system voltage. Another option is to use power electronics to | 
convert the machine EMF to the system voltage thereby providing steady regulation. 

Another effect that is sometimes an issue is armature reaction. It is caused by current 
flowing in the stator windings which creates a magnetic field that tends to distort the magnetic 
field established by the permanent magnets. For surface-magnet machines, because the magnet 
recoil permeability is approximately one and the magnet height is large compared to the air gap, 
the armature reaction flux density is small. As long as the stator teeth are not highly saturated 


due to the permanent magnets, armature reaction is negligible [17], [18]. 


3.3.5 Machine Inductances 
In a slotted permanent magnet machine, there are three distinct components of 

inductance: air gap inductance, slot leakage inductance, and end-turn inductance. The most 
accurate means for calculating these parameters is finite element analysis but analytical methods 
provide almost as good results (within a couple of percent) and are used in this paper [33]. 

. The air gap inductance is usually the largest portion of the total inductance and it is due 
to the interaction of the stator windings with the flux crossing the air gap. To calculate the air 
gap inductance, a full-pitch, concentrated winding carrying a current I is initially examined 


which leads to an air gap flux density shown in Eqn 3-23. 


Brjux = > B,,:sin(np6) 
n=1 4 Ho NI 
where B = —.—————.—— 
n odd nnn (g + hm) 2p 


Eqn 3-23 
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When this concept is expanded to polyphase windings with balanced operation, the air gap flux 


density becomes Eqn 3-24. 


Bflux = > B,,:sin (np@) 


n-l q 
where p = +.—.——_—.. 
n odd n 5 an (g + hm) 2p 


Eqn 3-24 
The flux can be found using equation Eqn 3-19 and the total flux linkage is) = N,®. With all 


real winding effects included, the air gap inductance is then given by Eqn 3-25. 


Eqn 3-25 
In addition to the air gap, the coil currents generate a magnetic field that crosses from one 
side of the slot to the other producing a slot leakage inductance. For calculating the slot leakage 
inductance, it is assumed that the slot is rectangular with slot depressions (Figure 13) which 


results in a slot permeance per unit length shown in Eqn 3-26 [17], [18], [19]. 


h h 
1 d 
Perm = Fe. K — 
W st Wd 


Eqn 3-26 
Assuming m slots per pole per phase and a standard double layer winding, it can be shown that 


the slot leakage inductance is given by Eqn 3-27 through Eqn 3-29 [19]. 
2 2 
ie 2-p-Ly-Perm| 4 No (m— Ny) + 2-Nyp'Ne | (self) 
Eqn 3-27 
Lim = 2p: L Perm NN. (mutual) 


Eqn 3-28 
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slot = Las ^ “am (3phase ) 


2 
Laia = Las — sial =) (higher odd phases ) 


Eqn 3-29 
The end turn inductance is the smallest of the three components. It is created by the 
magnetic field that surrounds a coil after it leaves one slot and before it enters another slot. 
Since it is extremely difficult to accurately determine because of complex winding patterns, a 
rough approximation is used. It is assumed that the end turns are semi-circular with a radius 
equal to one-half the mean coil pitch. Using reference [17], the total end turn inductance per 
phase is shown in Eqn 3-30. 
Bg NeNSA, TET 


———————— s lD 
Y 2 ZA 


Eqn 3-30 
The total inductance for the phase is the sum of the three inductances, ignoring other small 


factors.. 


L; = Lag + Lslot f Le X, = Og Lg 


Eqn 3-31 


3.3.6 Basic Losses 
Losses in a machine consist of core losses, conductor losses, friction and windage losses, 
and rotor losses. Rotor losses include magnet losses and retaining can losses, they require 


detailed waveform and harmonic analysis, and therefore are discussed later in Chapter 5. 


3.3.6.1 Core Losses 
High-speed generator stator core losses (per weight) can be greater than normal machines 
because of the higher frequencies. These losses are minimized by using laminated steels in the 
stator construction as discussed in 3.1.2 and by not generating frequencies that are too high. 
Core losses consist of hysteresis and eddy current losses. Hysteresis loss results from the 
steel not wanting to change magnetic state. As the flux density varies, the material traverses the 


B-H curve and energy is lost. Eddy current loss is also caused by the variation in flux density. 
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Electrical currents are induced within the ferromagnetic material and circulate dissipating power 
because of the resistivity of the material. 

Because there are usually various imperfections in materials, the best way to approximate 
core losses is to use empirical loss data. If the flux density is estimated for each part of a 
machine and the mass of the steel calculated, empirical core loss data can be used to estimate the 


total losses. Empirical data for M-19, 29 gauge material is obtained as shown in Figure 19. 


Mon-Orieniod Silicon Steele 
M-19 FP 29 Ge 100 Bz; 
— s = MELD FR 29 Ga 200 Hz pu 
— — — MADR AN GRADE o 
M-19 FP 29 Ga 600 Hz 
--— 4 M-19 FP29 Ga 1000 Hz ° 
` HA9 PP 29 Ga 500 Hz L 


MAGNETIC MATERIALS COMPARISON" 


Proto Laminations, Inc., 13666 East Bora Drive, Santa Fe Springs, California 90670 USA 
tel: (562) 926-4777 fax: (562) 404-1650  www.protolam.com 


Figure 19: Core Loss Data 


An exponential curve fit is then applied to the data to obtain an equation for estimating the core 
losses (Eqn 3-32) [9], [35], [36]. 


Eqn 3-32 
The values for the bases and exponents in Eqn 3-32 are listed in Table 8. 
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Table 8: Core Loss Parameters 


NUTUS 
. 


3.3.6.2 Conductor Losses 


Conductor losses arise from the current flowing in the stator windings. The resistance of 
the windings is calculated using Eqn 3-8 from section 3.3.2. The conductor losses are then found 


using the traditional power equation for a resistance (Eqn 3-33). 


Eqn 3-33 
3.3.6.3 Friction & Windage Losses 
For rotors operating at high-speed, friction and windage in air can cause losses which 


result in inefficiency and heat production. These losses are calculated using the power necessary 


to overcome the drag resistance of a rotating cylinder as given by Eqn 3-34 [32]. 


3 4 
Pwind = Cf p air 9 R La 


where C: = friction coefficient 


pair = density of air 


Eqn 3-34 
The friction coefficient depends on numerous factors such as surface roughness and flow region. 
Since the air gap is a small annulus and the rotor is spinning at high speed, it is assumed that the 
air in the gap is in the turbulent region. Therefore, the coefficient of friction can be 
approximated by Eqn 3-35 [34]. 
Cy = 0.0725 Rey 020 
where Rey = Reynold's Number 


Eqn 3-35 
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3.4 Machine Sizing Methods 


3.4.1 Basic Sizing Method 


Whenever a machine is being designed, it is important to perform some back-of-the- 
envelope calculations to gain insight into initial sizing estimates. Most generators are 
constrained by two competing design parameters, air gap magnetic shear stress and rotor tip 
speed. Air gap magnetic shear stress (t) is the magnetic shear force developed per unit gap area 
and is constrained by magnetic design and thermal management [8]. It is proportional to the 
product of the surface current density and magnetic flux density as shown in Eqn 3-36. 

T e Kz- B, 

where 7 = shear stress (psi) 
Kz = surface current density 
B; = air gap flux density 


Eqn 3-36 
Typical values for air gap shear stress for different types of generators are shown in Table 9 [8], 
[11], [16], [18], [31]. For the basic sizing calculations, 15 psi is assumed since the generator is a 


large liquid-cooled machine. 


Table 9: Air Gap Shear Stress Values 


i 1 
Large Air-Cooled 
Large Liquid-Cooled 10 — 20 psi 
High-Temperature Superconducting 30 — 50 psi 


Small Air-Cooled 


Rotor tip speed is discussed previously in section 3.2.2 and given by Eqn 3-1 with the 
assumed limit of 200 m/s. The fundamental machine power equation is utilized to derive the 


rotor radius and stack length of the machine (Eqn 3-37). 


P=2ÆR Lg tVip 
where R = rotor radius 
L4 = stack length 


Eqn 3-37 
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In order to simplify the equation, the L/D ratio is substituted in for Ls. Using the air gap shear 
stress, the rotor tip speed limit, and the power rating of the machine (16 MW ), the basic power 
equation is iterated to obtain a final L/D ratio, rotor radius, and stack length while matching the 

. desired rotational speed of the machine (13,000 RPM). Using an air gap flux density of 0.8 T, a 
pole pair value of 3, a slot height of 15 mm, and a slot fill fraction of 0.5, the frequency and 
current density of the machine are found. The detailed MATLAB code is contained in Appendix 
B with the results shown below. 


Basic Machine Design 


Input Parameters: 

Power = 16000.0 kW Shear Stress = 15.0 psi 
L/D Ratio = 2.85 Tip Speed = 200.0 m/s 
Pole Pairs = 3.0 Air Gap Bg = 0.80 T 
Output: 

Rotor Radius = 0.147 m Stack Length = 0.838 m 
Speed = 13000 RPM Frequency = 650.0 Hz 
Ja = 1757.67 A/cm2 


3.4.2 Detailed Sizing Method One 


Once basic sizing of the machine is complete, in-depth analysis is conducted to ascertain 
the overall performance and scale of the 16 MW generator. Two detailed sizing methods are 
developed using MATLAB code. To help gain an understanding for current PM machine 
designs, numerous articles, transactions, and proceedings are examined and a database is 
developed containing comprehensive specifications for over 20 different PM machines (see 
Appendix C). This database is used to help develop input parameters for the first method (Table 


10), although many of these values are obtained after numerous iterations through the code. 


Table 10: Input Parameters for Sizing Method 1 


rameter - [^ "Parameter > - 


Number of Phases 
Rotational Speed 13,000 RPM Number of Slots 


| Power Factor Angle — | Odeg |SlotsShort-pitched — | 
Stator Back Iron Ratio 
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Winding Conductivit 60x10 S/m | 
7700 kg/m 8900 kg/m 


7400kghmm | — 5 0 0 0 [| | 


Once the input parameters are entered, the first step in sizing is to generate the geometry 
of the machine. This includes determining the number of slots per pole per phase (Eqn 3-5), 
number of armature turns (Eqn 3-38), tooth width, slot dimensions, stator back iron dimensions, 
coil pitch, and winding end turn geometry. 
N, = 2:p-mN, 
where N. = Turns per coil 
N. assumes each slot has 2 half coils 


Eqn 3-38 
Next, the electrical frequency and rotor surface speed are determined using Eqn 3-1and Eqn 
3-39. 
p:N 
fs—— oz2Tf @ = p-@ 
60 Pm 


where œ = electrical frequency (rad/sec) 
Om = mechanical frequency (rad/sec) 
Eqn 3-39 

Winding, skew, and magnetic gap factors are then estimated as discussed in sections 3.2.3 and ` 
3.3.3. The magnitude of the air gap magnetic flux density (Bg) is determined accounting for 
slots, varying reluctances, and flux leakage per Eqn 3-14 through Eqn 3-17. 

The magnetic flux and back voltage magnitudes can then be calculated. The fundamental 
component of the magnetic flux links the stator windings to create useful voltage. Therefore, 
only the fundamental components of Eqn 3-21 and Eqn 3-22 are used to determine the internal 


voltage of the generator as shown in Eqn 3-40 and Eqn 3-41. 


4 T pD m 
mo 5 | 


Eqn 3-40 
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Az 2-Re Ley Ny kyrkg By E, = ya 
p 
Eqn 3-41 

Utilizing equations from section 3.3.5, the machine inductances and reactances are 
established. The lengths, volumes, and masses of components and the overall generator are 
calculated using basic geometric equations and relationships. A 15% service mass fraction is 
added to the total mass estimate to account for the additional services associated with large 
liquid-cooled machines [31]. Once the mass of each of the stator parts is known, the core losses l 
are estimated in accordance with section 3.3.6.1. 

The terminal voltage and current of the machine must then be calculated accounting for 
conductor losses and windage losses. The vector relationship (Figure 20) between terminal 
voltage (Va), internal voltage (Ea), and the synchronous reactance voltage drop is utilized to 
obtain Eqn 3-42. The armature resistance is usually ignored because it is much smaller than the 


synchronous reactance. 


Figure 20: Voltage Vector Relationship 


2 2 š 
Va? IE, - [X 1 cos) -X ysiny 


Eqn 3-42 
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The machine efficiency (n) is then easily obtained. Appendix D contains the detailed MATLAB 


code and the results are shown below. 


3.4.3 


PM Machine Design, Version 1: Surface Magnet, Slotted Stator 


Machine Size: 
Machine Diameter = 
Rotor radius = 
Slot Avg Width = 
Back Iron Thick = 
Machine Ratings: 
Power Rating = 
Va (RMS) = 

Ea (RMS) = 

Synch Reactance = 
Stator Cur Den = 
Efficiency = 
Phases = 

Stator Parameters: 
Number of Slots = 
Breadth Factor = 
Tooth Flux Den = 
Slots/pole/phase = 
Rotor Parameters: 
Magnet Height = 
Air gap = 

Magnet Remanence = 
Magnet Factor = 
Machine Losses: 
Core Loss = 
Windage Loss = 
Machine Weights: 
Core = 

Magnet = 

Services — 


Detailed Sizing Method Two 


0.472 m Machine Length = 
0.147 m Active length = 
16.493 mm Slot Height = 
34.300 mm Tooth Width = 
16000.0 kW Speed = 

2341 V Current = 

2925 V Arm Resistance = 
0.768 ohm Synch Induct = 
2199.7 A/cm2 Tip Speed = 
0.992 Power Factor = 
3 Frequency = 

36 Num Arm Turns = 
0.966 Pitch Factor = 
1.59 T Back Iron = 
2.00 

25.00 mm Magnet Angle = 
4.00 mm Pole Pairs = 

1.20 T Aig Gap Bg = 
0.949 Skew Factor = 
11.6 kW Armature Loss = 
30.7 kW Rotor Loss = 
396.49 kg Shaft = 

64.74 kg Armature = 
148.33 kg Total = 


1.003 m 
0.838 m 
25.000 mm 
15.403 mm 


13000 RPM 
2267.5 A 
0.00526 ohm 
0.188 mH 
200 m/s 
1.000 

650.0 Hz 


12 
0.966 
1.14 T 


50.0 degm 
3 


0.80 T 
0.989 


81.2 kW 
TBD kW 


438.05 kg 
89.58 kg 
1137.18 kg 


In order to provide a check on the methodology of the previous sizing procedure, a 


second MATLAB code is constructed. The second method is developed using a combination of 


processes from references [2], [17], [18], [19], [26], and [37]. 
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Table 11: Input Parameters for Sizing Method 2 


| Power Factor Angle | ^ Odeg [SlotvPole/Phase(m | 2 | 


Winding Conductivit 6.0 x 10° S/m 
8900 kg/m 


[A 
n : 
== 


7700 kg/m ° 
7400 kg/m 


Many of the same equations used in the first method are utilized but the process differs in some 
of the input parameters and calculations. The second sizing method inputs are shown in Table 
11 with different entries from method 1 being vip, Bsat, Ws, m, Ja, and PC. 

First, given the maximum Vtip (200 m/s) and the rotational speed (13,000 RPM), the 
electrical frequency and rotor radius are computed. Next, the winding and skew factors are 
determined similar to the first sizing method. The magnet dimensions (hm and Om) are then 
determined along with the air gap flux density. Using Bsa, the PC, and Eqn 3-43 through Eqn 
3-46, the magnet arc width and height are iterated until the tooth width (w;) equals the average 
slot width (ws). This differs from the first method where the peripheral tooth fraction is used to 


set the tooth width equal to the slot top width (wst). 


C, = 
$ I+ Om 
where am = magnet pitch coverage coefficient 
Eqn 3-43 
h m= 8 eC g PC 


Eqn 3-44 
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Kr C, 
=— B , xs 
Bg 7 r where yj = recoil permeability 


u 
1+ Ky rec B, = remnant flux density 

Eqn 3-45 

TD, B 
i —Ë where D, = 2*(R+h,,+9) 
t N, B 
S sat 
Eqn 3-46 


The permeance coefficient (PC) normally varies between 5 — 15 with higher PCs only 
used for high performance space and aircraft applications [18], [37]. A PC value of 4 — 6 is 
typical for large PM machines where the air gap flux density is approximately 60 — 80% of the 
remnant flux density [17]. The PC is input as 5.74 to coincide with the PC calculated in the first 
sizing method. 

Once the magnet dimensions are known, the geometry of the machine is generated using 
similar equations as the first method except for the stator core back iron depth. In the first 
method, the input variable “stator back iron ratio" is used whereas in the second method, the 
saturation flux density is applied as shown in Eqn 3-47. 

TD,.04 B 


d.= .— Š 
4p B 


sat 


c 


Eqn 3-47 
Utilizing a power balancing procedure, the stack length, terminal voltage, and terminal 
current are iterated to obtain a complete design. Core losses, conductor losses, and windage 
losses are calculated in accordance with sections 3.3.6.1, 3.3.6.2, and 3.3.6.3. The power 


crossing the air gap of the machine is determined using Eqn 3-48 [2], [37]. 


l f 
Paap = A7 ke krkyrkg sin(05) (Kz Bg). La 


where k, = electrical power waveform factor 
k; = current waveform factor 


Eqn 3-48 
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The electrical power waveform factor is defined in Eqn 3-49 and the current waveform factor is 
shown in Eqn 3-50. These factors depend on the back EMF and current waveforms generated in 


the machine and a detailed discussion of these factors is in reference [2]. 


T ` 
bat 20 
VOTI Elg 
0 
Eqn 3-49 
k.= lpk 
i I rms 
Eqn 3-50 


During the balancing procedure, the machine inductances and reactances are established 
per section 3.3.5. The lengths, volumes, and masses of components and the overall generator are 
calculated using basic geometric equations and relationships similar to the first sizing method. 
Once the balancing procedure is complete, the efficiency of the machine is determined. 
Appendix E contains the detailed MATLAB code and the results are shown below. 


PM Machine Design, Version 2: Surface Magnet, Slotted Stator 
Machine Size: 


AC Generator 


Machine Diameter = 0.484 m Machine Length <= — 0.979 m 
Rotor radius = 0.147 m Active length = 0.813 m 
Slot Avg Width = 16.000 mm Slot Height = 25.000 mm 
Back Iron Thick = 39.587 mm Tooth Width = 15.863 mm 
Machine Ratings: 

Power Rating = 16000.0 kW Speed = 13000 RPM 
Va (RMS) = 2604 V Current = 2064.1 A 
Ea (RMS) = 3007 V Arm Resistance = 0.00534 ohm 
Synch Reactance = 0.723 ohm Synch Induct = 0.177 mH 
Stator Cur Den = 2064.1 A/cm2 Tip Speed = 200 m/s 
Efficiency = 0.993 Power Factor = 1.000 
Phases = 3 Frequency = 650.0 Hz 
Stator Parameters: 

Number of Slots = 36 Num Arm Turns= 12 

Breadth Factor = 0.966 Pitch Factor = 0.966 
Tooth Flux Den = 1.70 T Back Iron = 1.05 T 
Slots/pole/phase = 2.00 

Rotor Parameters: 

Magnet Height = 26.17 mm Magnet Angle = 47.7 degm 
Air gap = 4.00 mm Pole Pairs = 3 

Magnet Remanence = 1.20 T Aig Gap Bg = 0.85 T 
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Magnet Factor — 0.958 Skew Factor = 0.989 
Machine Losses: 


Core Loss = 11.9 kW Armature Loss = 77.5 kW 
Windage Loss = 29.7 kW Rotor Loss = TBD kW 
Machine Weights: 

Core = 438.62 kg Shaft = 424.57 kg 
Magnet-  . .62.89 kg Armature = 85.50 kg 
Services = 151.74 kg Total = 1163.32 kg 


3.4. Comparison of Methods 

Both sizing methods produce similar sized generators with the results agreeing within 1 — 
5% on most parameters. The main reasons for the differences are twofold. First, the calculations 
for the magnet dimensions are different. In the first method these parameters are input whereas 
in the second they are determined through an iterative process using material properties. Second, 
the procedure for determining the slot and tooth widths differs slightly. 

Overall, both methods underestimate the overall dimensions (length and diameter) and 
weight. This is because the sizing programs do not include calculations for portions such as 
structure, frames, mounts, and maintenance access. These pieces are added in section 6.1 to get 
the total sizes and weights of the PM generators, specifically including: 

e 50% weight factor added to total weight 
e 0.8 m added to overall length and 0.4 m added to overall diameter 
In later chapters when waveforms are examined and machine optimization occurs, the first sizing. 


method is mostly employed with the second method used for verification. 
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Chapter 4 Power Electronics and Conversion 


41 Background 


In conjunction with employing a PM generator, a DC bus architecture is one of the 
preferred schemes for the future [5]. DC power distribution systems can offer an advantage in 
size and a weight over high-power AC systems [3]. Therefore, a power electronics module is 
required to rectify the AC generator output and then convert it to the appropriate DC distribution 

voltage of 700 VDC. However, conversion of the high-frequency generator AC output to DC 
generates harmonic voltages and currents on the AC side and ripple on the DC side. These 
harmonics can be reduced through proper design of the conversion module. 

Since the generator is a permanent magnet machine, the field cannot be controlled like in 
a wound-rotor machine. The power electronics module (PEM) is needed not only to convert the 
machine EMF to the DC system voltage but also provide steady regulation. Also, since the 
frequency of the generator output is higher than the normal 60 Hz systems, the filter 
requirements for the PEM are less since the components can be smaller. 

Many options exist for the type and model of the PEM, and a combination of a rectifier 
and DC-DC buck converter are chosen. This design provides a two-step conversion from the 
generator to the DC bus and allows for good DC bus voltage regulation while reducing 
harmonics reflected back to the generator. A block diagram of the system is shown in Figure 21. 
The calculations for the rectifier and buck converter are performed for the ideal case and ignore 


parasitics such as equivalent series resistance, inductance, and capacitance. 


Rectifier 


PM Generator 


Figure 21: Basic System Layout 
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42 Rectification 


The rectifier stage converts the AC voltage from the PM generator to a DC voltage. 
There are two main choices when designing the rectification stage: active or passive rectification 
and series or polyphase rectification. Active rectification uses controllable components, such as 
thyristors or MOSFETs, to actively rectify the AC voltage. This increases the cost and 
complexity of the rectifier but provides voltage regulation and control. Passive rectification 
consists of diode bridges to convert the AC to DC while relying on some other circuit to perform 
the voltage regulation. Å passive rectification scheme is used in this paper because of its 
simplicity, potential for lower input harmonics, and lower cost. 

Series rectification is when each phase is rectified independently from the others and the 
output DC voltage is obtained by a series connection of each single phase rectifier. Each phase 
conducts during the full period resulting in high power losses and higher stresses on components. 
Also, the output DC voltage is extremely high since each rectified phase voltage is added 
together. Parallel rectification is when the lines with the most positive and negative 
instantaneous voltages provide the forward bias to turn on two diodes. The output voltage 
waveform corresponds to the instantaneous difference between two line voltages. For sinusoidal 


voltage supplies, the equation for the output DC voltage is approximately Eqn 4-1 [38]. 


2.n T 
Vpc —'sin| — |: V. 
DC " ( ü J S 
where V, = peak amplitude of phase voltage 
Eqn 4-1 
The losses are lower with parallel rectification and its output voltage provides a better 


match with the buck converter. Therefore, a parallel rectification scheme is used in this paper. 


For the initial PM generator from section 3.4.2, the rectifier is shown in Figure 22. 
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Neutràl 


Figure 22: Basic 3-phase Rectifier 
The DC-side inductor and capacitor are necessary to obtain a good DC output voltage 
with low ripple (large C) and an acceptable AC line current waveform (large L) [39]. The goal is 
to produce maximum power factor with minimal line-current harmonics. A power factor of 0.96 
with low line harmonics can be achieved by selecting the inductor and capacitor values using 
Eqn 4-2 and Eqn 4-3 [40]. 


V. 

f 

L.= Lon: - 
Lef "ref 


where Lon < 0.10 
V ef = nominal |-n RMS source voltage 
ler = P ret! V ret (RMS current from V ep) 
P ef = nominal output power 


La = nominal source frequency 


Eqn 4-2 


- [scc 
Vref fret 
where Cau = 100.0 
V ef = nominal l-n RMS source voltage 
La = P ,efV ret (RMS current from V ef) 
P, = nominal output power 
La = nominal source frequency 
Eqn 4-3 
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For the initial PM generator from section 3.4.2 where Var = 2925 V, Prep = 16 MW, and fjer = 650 
Hz, the values for the DC-side inductor and capacitor are L, = 82.3 uH and C, = 287.7 mF. 


These values will change as the system is optimized. 


43 DC-DC Conversion 

After the rectifier converts the AC voltage to DC voltage, the DC voltage must be 
stepped down to the bus voltage of 700 VDC. A high-power buck converter is ideal for this 
application. The power circuit has a basic topology to which other circuit components are added. 
to perform functions such as filtering and over-voltage protection [41]. The overall goal is to 

| alter electrical energy provided by an input system to that required by an output system. 

The system consists of a two-stage input filter, a converter, and a low-pass output filter as 

shown in Figure 23. Duty ratio control using a feedback controller could be used to provide DC 


bus voltage regulation during normal operation. 


FROM 
‘RECTIFIER 


INPUT FILTER CONVERTER OUTPUT FILTER/OUTPUT. 


Figure 23: Basic Buck Converter 


4.3.1 Buck Converter 

The buck converter is considered a high-frequency DC-DC switching converter. The 
difference between the switching frequency and the frequency of the external waveforms allows 
the use of low-pass filters in the input and output to help improve performance [41]. The buck 
converter is designed to take a high DC voltage and reduce it to a lower DC voltage under 
varying loads and over a range of input voltages. The level of voltage reduction is controlled by 
the high frequency switching. The switch opens and closes at a specific frequency with the ratio 


of on-time to the period defined as D, the duty ratio. 
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For the buck converter, the voltage conversion ratio is dictated by Eqn 4-4. 


Vout lou 


SIL 


t 
in n 
Eqn 4-4 

It is assumed that the required output power from the generator and PEM can be as high as the 


rated 16 MW and as low as 0.1 MW. Based on this, output currents and resistive load values are 


determined and are listed in Table 12. 


Table 12: Buck Converter Load Values 


M 
i Ses nins 


| 129A | 4900 | 
16MW 700 VDC 22.8571 Å 0.0310 


The switching frequency for the converter is selected to be 50 kHz (T = 20 usec). This is 
sufficiently fast to ensure passive components are not excessively large but also not too high 
speed so that there are not excessive switching losses. Additionally, current switching ratings for 
advanced, high-power IGBT devices are limited to the hundred kilohertz range. By selecting 50 | 
kHz, this provides a margin below this maximum switching frequency and ensures devices are 


presently available or will be in the near future. 


4.3.2 Output Filter 

A low-pass output filter is required to help reduce the output voltage ripple to within the 
required specification listed in Table 2. The output voltage from the converter normally 
fluctuates between 0 and Vour, so the filter, consisting of an inductor and capacitor, is used to 
minimize this problem. It is important that the corner frequency of the filter be much lower that 
the switching frequency (fc «« fsw) so that the switching frequency ripple is eliminated in the 
output voltage [42]. Normally, a factor of at least 100x is used [43], so to ensure an adequate 
margin for parasitics and fulfill the strict DC bus voltage ripple requirement, a factor of 400 is 


selected so that fc ~ 125 Hz. For an LC-filter, Eqn 4-5 holds. 


l 1 
f c = ———— 
E 2l LC 


Eqn 4-5 


65 


AC Generator Top Magnetic Generator Magnet Generator DIY 


Therefore, given the value selected for fc: 


LC= 16210 ° 


Eqn 4-6 
With this equation the output filter inductor and capacitor are sized using a buck 
converter sizing program [44]. It is assumed that an inductor current ripple of 0.1% at maximum 
power (16 MW) is satisfactory for sizing the inductor and Eqn 4-7 is used where Vg = 1.9 V, the 
forward-bias diode voltage [44]. 


| E 1 Your * VF 1 
HI EAS Vo max * YR | AL 


in max 


Eqn 4-7 
For the initial generator of section 3.4.2, this produces an inductor size of 551 uH. Using this 
inductance along with Eqn 4-6, a capacitor value of 2.94 mF is calculated. 
In order to validate the filter's performance, a rough estimate for the voltage ripple is 
calculated. It is assumed that the ripple component of inductor current flows through the 
capacitor and the average component of the output current flows through the load resistor [42]. 


Therefore, the peak-to-peak voltage ripple can be written as shown in Eqn 4-8. 


Eqn 4-8 
The output voltage ripple is estimated to be 0.019 V which is within the specification listed in 
Table 2. This overall process is used later to design the converter output filter as the PM 


generator system is optimized. 


4.3.3 Input Filter 

The input filter is needed to prevent the converter switching fluctuations from affecting 
the voltage and current from the rectifier. The filter selected is a two-stage filter with an L-C 
filter and an R-C damping leg. From the switching frequency, the angular frequency is œ = 
3.142x10? rad/sec. In order to examine worst case, the 16 MW, 22,857 A-output is used, and 


with the corresponding duty ratio of 0.1023 for the generator of section 3.4.2, an average current 
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of 2,338.3 A is calculated. Since the fundamental of the ripple is attenuated much less than the 
higher harmonics, the goal is to suppress the fundamental of the ripple along with some margin. 


The fundamental of the average current is calculated in Eqn 4-9. 


iconv! = 5 < icon >= i. (2338.3) = 2977.2Amps 
m m 


Eqn 4-9 
It is assumed that the buck converter input current ripple must be kept very small (« 10 A at 
maximum power) to ensure no detrimental effects on the rectifier. Therefore, at the 50 kHz 


switching frequency, the maximum gain of filter is G as calculated in Eqn 4-10. 


10.0Amps 
. ripple —T 
G= Shu 2.2. -0.00168 
iconvi š 2977.2 


Eqn 4-10 
In order to account for the contribution of higher harmonics and non-idealities associated with 
actual components, a 50% margin is used for G resulting in G = 2.52x102. 

The output impedance of the input filter is also a concern because it must be sufficiently 
small so that it does not adversely affect the converter's performance. It is assumed that the 
output voltage ripple of the filter must be less than 10 V at maximum power. From this, a value 
of Zo at 50 kHz is calculated (Eqn 4-11). 

10.0 V 


Zo= = 0.00168 


9 29772 


Eqn 4-11 
Using the same margin of 50% for the same reasons, the output impedance is 1.12x10°. 
Initially, the L-C low pass filter portion is examined to calculate values for Lr and Cp. 
Then, a well-designed filter has Cy ~ 10*C;, and Rr is determined using reference [39]. The 


calculations are listed below. 
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iia ee 
1+ s?LC 
1.1210 3 + 1.1210 ^s. LCS s-L 
s-L= 0.444 
L< 1415uH 


2.5210 °= 
ida LE 


25210 > + 2.5210 ?.s”LC= 1 


Rr [2+ n)-(4 + 30) 
Ro 2.n2.(4 + n) 


R£ = 0.215 :R, =4.79-10 


Choose 10*Rf to be safe 


S-LC = 395.83 
s-C = 891.51 Cp = 2.84mF 
Eqn 4-12 
Lr 
R, = |— = 0.0223 
Cr 
3 
Rẹ = 479 mQ 
Eqn 4-13 


The performance of the input filter is determined by examining the filter's transfer 


function and verifying it has sufficient attenuation at the switching frequency and also minimal 


peaking. The transfer function for the filter is given in Eqn 4-14 [45]. 
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(RrCy:s) +1 


FREE + ls, " cd + (RpCy s) +1 


Eqn 4-14 
Figure 24 shows the Bode plot for the initial input filter and it operates satisfactorily (Appendix 
F contains the detailed MATLAB code). This overall procedure for designing the converter 


input filter is utilized later to as the PM generator system is optimized. 


System: H 
Frequency (radisec): 3.07e+005 
Magnitude (dB): -51 8 


System: H 
Frequency (rad/sec): 4.55e«003 UM ^ ve | 
. Magnitude (dB): 6.19 


Figure 24: Bode Plot for Converter Input Filter 


4.3.4 Converter Control 

One control technique that can be used is duty ratio control utilizing a feedback 
controller. This arrangement is presented as an illustration of one control methodology for the 
power electronics module but implementation is beyond the scope of this paper. 

The control scheme involves negative feedback control using linearized, averaged state- 
space models. Figure 25 shows a block diagram of this where K represents the controller and H 


represents the buck converter. It is designed so that the closed-loop system is stable and well 
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damped for the input voltages and output loads. Voltage-mode control (duty-ratio control) is 


implemented using a Proportional Integral Derivative (PID) controller. 


Figure 25: Block Diagram of Control Scheme 


The majority of control systems used in industry applications are based on this type of 
controller. In general, proportional and integral control can be applied independently or together 
(P, I, or PD, while derivative control can only be applied in combination with one or both of the 
other two (PD, PID). Testing of a PID controller normally involves plotting the system response 
to a simple input, such as a step input, and adjusting the gains based on the output. 

The first step in developing the control scheme is to obtain a linearized, averaged model 
for the buck converter and derive the transfer function from perturbation in duty ratio to 
perturbation in output voltage. The transfer function represents H in Figure 25. Next the PID 
controlled is developed, consisting of three gains, Kp, Ky, and Kp, and this represents K in Figure 
25. The H and K are then combined to form the complete transfer function for the system. 

Initially, the state variables are assigned (ir and vc) and the state-space equations are 


determined (neglecting the input filter). 


d. d ..9€ "HE. < 
p - (vi = vc)'a(t) - vod — q(t)) ire = Ë - ao == Ë _ lu - q(t)) 


d. 1 d 1f, "c 
— = —. . t)- — = — _ — 
ar L (vato - vc) ac Ë R J 
Eqn 4-15 
Then, the circuit is averaged and it is assumed that the state variables have slow variation and 
small ripple, and d = avg(q(t)). All of the variables are now averaged quantities and Eqn 4-15 


becomes Eqn 4-16. 
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d i -l(v d - vq) PEE CERE 
ae py! C ac C R 


Eqn 4-16 
The variables are linearized resulting in the equations shown in Eqn 4-17. 
YGS Vc + Vc d=D+d' 
H 271p ip v| = Vi +v] 
Linearized Equations 
V. 
d. 1 
Zh = —v' + —-d' — —- v' 
qt L l L c 
d v O2—dp- EN 
d C € Re © 
Eqn 4-17 


Since vc is the output voltage, utilizing the Laplace operator and solving for the transfer function 


yields Eqn 4-18. 


s 2 + L., + 1 y' D, + M d 
RC no) xc f' EC 
M: 
He Vout LC 
f 2 1 1 
s + —:s +t — 
R-C LC 
Eqn 4-18 
For the controller the transfer function is Eqn 4-19. 
E n E ux Kr where e'is Viet- Vout 
e p S 
K * K K 
t D's + p's + I 
s 
Eqn 4-19 


Overall, the total transfer function is then Eqn 4-20. 
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K (x g^ Kp:S + KI) 
Vou LC N D P I 


Eqn 4-20 
Using the total system transfer function, the gains on the controller are adjusted to obtain 
the desired steady-state and transient operation. This type of system can readily be designed and 


optimized to provide the required overall performance, but this is beyond the scope of this paper. 


44 Conversion Losses 


An integral part of a DC-DC converter's performance is how efficiently it operates 
during the energy conversion process. Losses that occur include winding losses due to 
conduction, switching losses, and inductor core losses. A detailed analysis of all of these losses 
is beyond the scope of this paper and they are only described qualitatively and estimated using 
analytical equations. 

Practical magnetic components, such as inductors, exhibit losses generated by winding 
conduction, hysteresis, and eddy currents. Resistive components in the converter as well as 
parasitics (equivalent series resistance) produce IR losses. Ideal inductors and capacitors are 
employed in this paper and parasitic losses are not calculated. 

The transistors and diodes have conduction losses based on their on-state resistances. 
The equations to estimate these losses are shown in Eqn 4-21and Eqn 4-22. Lastly, the 
transistors have switching losses associated with the transistors changing states. If the voltage 
and current waveforms are assumed to rise and fall approximately linearly, the switching losses 


are determined as shown in Figure 26. 


2 
Pcond = (u. D] Ras 
where I = average current 
D = duty cycle 
Ras = IGBT on-state resistance 


Eqn 4-21 
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2 
Pdiode = (1 - D) I, ‘Ra 


where |, = average current 
D = duty cycle 
Ry = diode on-state resistance 


Eqn 4-22 


Figure 26: Transistor Switching Losses 


4.5 Component Sizes and Weights 


The PEM detailed in the previous sections consists of IGBTs, diodes, capacitors, 
inductors, and resistors. The sizes and weights of the resistors are neglected because they are 
significantly smaller than the other components. To ascertain reasonable estimates for the other 
devices, numerous references are utilized [15], [46], [47], [48], [49], [50], [51], [52], [53], [54]. 
In addition, a 30% service fraction is added to the total component mass and volume calculations 
to account for heat sinks, cooling, etc. | 

The power transistor is the main switching device in the converter circuit. In many 


applications, IGBTs are used because they have the current density and low loss of bipolar 
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transistors with the high-speed and high input impedance of MOSFETs [54]. Power diodes are 


needed to complete the buck converter circuit and provide proper operation. Table 13 lists the 


assumed device parameters for the IGBTs and diodes used. The IGBT and diode values are not 


from exact devices from particular manufacturers but rather represent data assembled from 


numerous sources. In order to achieve the high power levels, the semiconductor devices must be 


connected in parallel. Specifically, since the IGBTs and diodes are rated for 1,200A each, 19 


modules (each rated for 845 kW) are connected in parallel to achieve the required 22.5 kA. This 
can be readily accomplished using IGBTs provided the switching is properly controlled [55]. 


The inductors and capacitors are needed to perform energy storage and filtering 


functions. In high-power applications they are usually the dominant components in terms of 


mass and volume. For estimating their contribution to the converter, component weight and 


volume energy densities are utilized. Table 13 lists the assumed characteristics for the inductor 


and capacitor sizing estimations. 


AC Generator 


Table 13: Power Electronics Module Component Characteristics 


Maximum Collector Emitter Voltage 


1200 A 
Turn Off / Fall Time 


Forward Voltage Drop 
i i | -100kHz | 


j Peak Blocking Voltage i E 6500 V 
Average Forward Current 1200 A 


0.1 kg/] 


Mass Densit 3 kg/J 
0.08 J/cc 
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In addition to using the component values to size the power conversion module, a weight factor 
(1.0x) and a volume factor (0.20x) are employed to account for portions such as structure, 


frames, mounts, and maintenance access. 
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Chapter 5 Waveforms, Models, and Machine/Module Optimization 


5.1 Initial Generator Waveforms 


The PM generator produces back EMF waveforms that are dependent on a number of 
factors as discussed in Chapter 3. The goal is to produce a voltage waveform that closely - 
resembles a sinusoidal waveform with a low total harmonic distortion (THD) because this results 
in minimal harmonic content which reduces losses in the machine. THD is a measure of the 


distortion in a waveform caused by undesirable frequency components. It is calculated as shown 


in Eqn 5-1. 


Eqn 5-1 
The back EMF waveforms are generated using Eqn 3-18, Eqn 3-19, and Eqn 3-22 which 


are shown again below. 


B(6) = y B, 'sin(np6) 


Om = Magnet physical angle 
n = harmonic number 


Eqn 5-2 
Mo) = Y; Aq sin(np8) 
n=1 where 3 = 2Ry Lar: Ny Ba Kun Kn 
n odd n 
P 
Eqn 5-3 
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E, = 2 V, sin(np9) 
n-i 
n odd dt n 
Eqn 5-4 
MATLAB code is developed to create the generator EMF waveforms as well as 
determine the harmonic content and THD of the generator output waveform. This code is 
contained in 0. For the initial PM generator designed in section 3.4, the air gap flux density, 
back EMF waveform, and harmonic content (line-neutral) are shown in Figure 27, Figure 28, and 
Figure 29. As a frame of reference, MIL-STD-1399 requires Type I 60 Hz AC systems to have a 
total THD of less than 5 % line-line [56]. This limit does not directly apply to the PM generator 
AC bus since it is being rectified to DC before supplying loads. However, it is used as a 
guideline because it provides a reasonable limit to ensure proper rectification and performance 


on the AC side. 


Figure 27: Initial Generator Flux Density Waveform 
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bu 


Above 10% of Fundamental P 
Below 1096 of Fundamental (s 


Figure 29: Initial Generator Harmonic Content 
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52 Rotational Stress and Retaining Sleeve 
Since the PM generator is spinning at high speed, the rotor and permanent magnets are 
subjected to extremely high centrifugal forces. These forces can cause significant amounts of 
damage if the magnets and rotor are not properly restrained. The rotational components can be 
strengthened by enclosing them in a retaining sleeve/can which also increases the air gap length. 
The centrifugal force on the magnets due to the rotor spinning is calculated in Eqn 5-5. 


2 
_ Ma mag 


F = 
cen 
Rocha 


where Mp = mass of magnets 


Vmag = Velocity of magnets (m/s) 


Eqn 5-5 
Using the inner surface area of the retaining sleeve, this force is converted to an outward 
pressure. Treating the retaining sleeve as a thin-walled vessel, the hoop stress felt by the sleeve 


is determined as shown in Figure 30 and Eqn 5-6. 


Figure 30: Retaining Sleeve Hoop Stress 


T 
* Fe=-2o;h'L+ | P-r-Lsin(@)-Ld@ = 0 
vert 0 


P-r 
o.=— 
h 


C 


Eqn 5-6 
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To provide a margin for mechanical tolerances and imperfections, a safety factor (SF) of 1.2 is 
applied to the hoop stress calculated in Eqn 5-6 to get a final SF stress for the retaining sleeve. 

The retaining sleeve can be made from many different types of materials including metal 
alloys and composites. Å disadvantage of a metallic sleeve is eddy currents are induced in the 
sleeve by variations in the flux density caused by the stator slots. The advantages of a metallic 
sleeve are that it shields the magnets from most of the flux density variations and it has a high 
thermal conductivity for heat removal. 

Å composite, such as carbon-fiber, provides reasonable strength while having lower 
losses since the lower conductivity reduces the eddy currents. However, the composite sleeve 
has low thermal conductivity and does not shield the magnets from the flux variations. This 
results in increased losses in the magnets themselves. In general, when the sleeve conductivity is 
low, the rotor losses are in the permanent magnets, while when the sleeve conductivity is high, 
the PM losses decrease and the sleeve losses increase. The materials considered for use as in the 
retaining sleeve are listed in Table 14 [10], [57], [58], [59], [60], [61], [62], [63], [64]. The 
numbers for the carbon fiber composite are in the middle of a range of typical values (0.5 — 18.0 


uQ-m) since wide variations exist depending on how the composite is made. 


Table 14: Retaining Sleeve Materials 


M th m) 
| Stainless Steel | 90 | — 072  — 


0 contains detailed MATLAB code to perform the retaining sleeve stress calculations. 


The initial PM generator from section 3.4.2 is analyzed with the results below. 


Retaining Sleeve Stress: 
Stress Limits: 


Stainless Steel = 90.0 ksi Aluminum Alloy 2 75.0 ksi 
Titanium Alloy = 110.0 ksi Carbon Fiber = 100.0 ksi 
Inconel = 132.0 ksi 


Actual Sleeve Stress: 
Sleeve Thickness Actual Stress SF Stress 


0.50 mm 1136.9 ksi 1364.2 ksi 

1.00 mm 568.4 ksi 682.1 ksi 

1.50 mm 379.0 ksi 454.7 ksi 
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2.00 mm 284.2 ksi 341.1 ksi 


2.50 mm 227.4 ksi 272.8 ksi 
3.00 mm 189.5 ksi 227.4 ksi 
3.50 mm 162.4 ksi 194.9 ksi 
4.00 mm . 142.1ksi 170.5 ksi 
4.50 mm 126.3 ksi 151.6 ksi 
5.00 mm 113.7 ksi 136.4 ksi 
5.50 mm 103.4 ksi 124.0 ksi 
6.00 mm 94.7 ksi 113.7 ksi 
6.50 mm 87.5 ksi 104.9 ksi 
7.00 mm 81.2 ksi 97.4 ksi 
7:50 mm 75.8 ksi 90.9 ksi 
8.00 mm 71.1 ksi 85.3 ksi 


It is evident that the initial 4 mm gap length is not sufficient once the retaining sleeve is 
considered (since the air gap dimension includes the retaining sleeve thickness). The machine 
must therefore be revised to allow for the retaining can. In this modification process, the 
aluminum alloy is removed as a sleeve material possibility since it has lower strength and much 
higher conductivity than the other materials which together would produce a larger air gap and 
higher eddy current losses. 

The PM generator will be redesigned taking the retaining sleeve and hoop stress limits 
into consideration. The updated PM generator will be a bigger machine with a much larger air 


gap, greater magnet height, lower Bg, lower voltage, and higher current density. 


5.3 Rotor Losses 


5.3.1 Model for Time Harmonics & Winding Space Harmonics 


The permanent magnets used in high-speed generators are electrically conductive and 
therefore support eddy currents. The retaining sleeves are sometimes made from electrically 
conductive material that also can carry eddy currents. These eddy currents are primarily caused 
by fluctuations in the magnetic flux density produced by time and space harmonics of the 
winding currents. The currents produce losses which can potentially cause excessive heating or 
demagnetization of the permanent magnets. An analytical model is developed using the winding 


and current harmonics and the surface impedance to estimate the rotor losses. 
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Figure 31 shows the PM generator geometry “flattened out" into rectilinear coordinates. 
This is an accurate representation provided the dimensions are on a radial scale that is much 


smaller than the radius of the machine so that curvature is not important [65]. 


| Retainin | Sleeve 
Maqnets 


Rotor Hub 


+ 


Figure 31: General Magnet Loss Model 


The direction of rotation is in the positive x-direction, the radial direction is y, and the armature 
current flows in the axial dimension, z. In addition to the arrangement of Figure 31, the 


geometry shown in Figure 32 is also utilized. 


Figure 32: Layer of Material 


The following assumptions are made in developing the rotor analytical loss model: 
e Layers of material extend to + infinity in the + x direction 


e Layers effectively extend to negative infinity in the negative y direction 
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e  Motion/rotation is in the + x direction 

e The physical constants of the layers are homogeneous, isentropic, and linear 

e The ferromagnetic material does not saturate 

e The machine is long axially so magnetic variations in the z direction are ignored (H and 
B only vary in x, y directions) 

° All currents flow in the z direction 

e The rotor and stator are constructed of laminated steel so their conductivity in the z 
direction is negligible 

e The time and space variations are approximately sinusoidal 

e Flux density at y = infinity is zero 

e A traveling flux wave harmonic can be represented by an equivalent traveling current 
sheet (K;) on the surface of the stator 

e The normal component of the flux density is continuous at all interfaces 

e The tangential component of the flux density is continuous at all interfaces except at the 
stator/air gap where it is increased by the current sheet density 

e The magnetic flux density crossing the air gap and magnets is perpendicular 

e The effect of magnet eddy currents on the magnetic flux density is negligible — this is 
accurate below 10 kHz [66] 

e The magnet flux density is constant over the magnet breadth 
One of Maxwell's Law states that there is no magnetic charge and therefore flux lines 


must close on themselves, represented by Eqn 5-7. 
> 
V.Bz0 

"e Eqn 5-7 
Knowing this, the flux density B can be represented as the divergence of the magnetic scalar 
potential A (Eqn 5-8). 

B=VxA 
Eqn 5-8 

Substituting in uoH for B, taking the gradient of Eqn 5-8, and using Ampere's Law for current 
density yields Eqn 5-9. 
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E 3 
HoH < Vx A 


3 > 
UyVXHzVxVxA 


ug VxH=vÅv-A) - v^à 


? 22 
ugJ2 -V-A 

? 22 

Up 0:E= -V -A 


Eqn 5-9 


Faraday's Law is then employed to determine a relationship between E and A as shown below. 


> > 
VxE=-43 (Faraday's Law) 
dt 
VxE=-Vx—A 
dt 
> > 
E= -d A 
t 
Eqn 5-10 
Combining Eqn 5-9 and Eqn 5-10 gives Eqn 5-11. 
27 d? 
V :'A=uryo—A 
O at 
Eqn 5-11 


Since H and B vary only in the x and y directions, A and E are in the z direction and A can be 


represented by Eqn 5-12. 


ET 
A E Rel (A4).e (at d 
where k = 2 x/), 
Eqn 5-12 


Performing the mathematical operations and solving, Eqn 5-12 becomes Eqn 5-13. 
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2 d 
V A, = Bo cna 


> > 
Since B=V xA andAisin z-dir, then B - A and p = Ån 
X z y Z 
dy dx 
2 2 2 2 
2 
V = aos TR + -— = -k Ta 
dx y dz dy 
2 d? ! 
-k -A,(y) + P iù = jog o: A (y) 
y 


a? 2 
—A;,(y) - (x * jongo] A0) -0 
dy 
| 2 . 
Let y= k + joo 
Solution of the form: 


Ay) = Aye! + Aye ” 
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Eqn 5-13 
Solving for the magnetic flux densities and magnetic fields produces Eqn 5-15. 
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Eqn 5-15 


The layer from Figure 32 can now be examined utilizing a “surface coefficient." A surface 
coefficient is defined as the ratio of the y-directed to x-directed magnetic field amplitude (a = 


Hy/Hx). At the bottom of the layer where y = 0, the surface coefficient is given by Eqn 5-16. 


Hy kp Du Dp 


d" 


Eqn 5-16 
At the top of the layer (y = h), the surface coefficient is shown in Eqn 5-17. 
-yh 
_k Hye” -H,e 
ET Sp n Hp yh_ -yh 
I -6 — e —6 
Hy e T He E H, 
Ot = L H, 
ET 
H, 
Eqn 5-17 


Using Eqn 5-16 to solve for the ratio H;/H,, a final expression for the top surface coefficient is 


determined (Eqn 5-18), and it is applicable to any uniform region. 
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Eqn 5-18 
If the region being examined is positioned on top of a ferromagnetic surface, such as the magnets 


on the steel rotor shaft, the boundary condition at the bottom of the layer (a — infinity as H, — 


0) produces Eqn 5-19. 


a= joan) 


Eqn 5-19 


In the case of the air gap where the conductivity is zero, Eqn 5-18 reduces to Eqn 5-20. 
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(jsinh(yh) + a cosh (yh) 
PE j cosh (yh) +O p sinh(yh) 


Eqn 5-20 


The electric field (E) is z-directed and the magnetic flux density and field (B and H) are 
y-directed so a relationship is determined between them using Faraday's Law. The surface 
impedance is then determined. Surface impedance (Z;) is the ratio of the z-directed electric field 


to the z-directed current. Eqn 5-21 and Eqn 5-22 show the detailed calculations. 


> 4> 

V x E= —B (Faraday's Law) 

dt 
> 
Vx E= =H 
dt 

d d d 
—E, — —E, = Une H 
dz X dx” Ñ dt Y 
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Eqn 5-21 
—(Q. 
n 
Z = Esurf E, k, 9 
s = — gg — € 
Ksurf -Hy -Hy 
Q 
n 

Zs = me 

Eqn 5-22 


Since electromagnetic power flow into the rotor is the desired quantity, Poynting’s vector 


is employed (Eqn 5-23). 


> > 
Ex H 


nv 
Il 


Eqn 5-23 


Dissipation in the rotor is in the negative y direction producing Eqn 5-24. 
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Sy E: +E, Hy - E,H,) 


There is no z-directed magnetic field (H , = 0) 


Sy = —E, Hy 
Eqn 5-24 
In time average form, Eqn 5-24 becomes Eqn 5-25. 
-< $ >= EL -Ho 
li, E; 
— < Sy >= E^ H. Re 

-<$ >= IK -Re(Z.) 

Eqn 5-25 


Eqn 5-25 yields the power dissipated at the stator surface. This is the correct result for 
the rotor because there is no mechanism for dissipating power between the stator and rotor. The 
power estimated by Poynting’s theorem flows directly from the stator to the rotor [65]. 

In order to calculate the rotor losses, the above analytical model is applied to the 
geometry in Figure 31. For this model, the stator is assumed to be a smooth surface without slots 
because the slot effects are considered later in section 5.3.2. The first step is calculating the 
surface coefficient at the bottom of the magnet layer. It is assumed that this is formed by the 
highly permeable rotor shaft below the magnets. This assumption allows the surface coefficient 
at the top of the magnet layer to be calculated using Eqn 5-19. 

Next, the surface coefficient at the top of the retaining sleeve is determined using the 
sleeve material’s conductivity and Eqn 5-18. Traveling across the air gap to the surface of the 
stator, the surface coefficient is estimated using Eqn 5-20. The surface impedance is then 
computed utilizing Eqn 5-22. 

Once the surface impedance is known, the rotor losses caused by the stator winding time 
and space harmonics are calculated using Eqn 5-25. Only the more significant harmonics are 
considered as shown in Figure 33. The higher order harmonic effects are ignored because the 


products of higher order time harmonics (small number) and higher order space harmonics 
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(small number) produce negligible contributions (small number * small number = very small 


number). Appendix I contains the detailed MATLAB code that performs all of the above 


calculations. 


Time Harmonics 
11 13 15 17 19 21 23 25 27 29 31 
X X X X X X X X X X X 


Higher order effects 
neglected because 
negligible 


Space Harmonics 
>x X X X X X X X x X X x x x x xi 


Figure 33: Relevant Harmonics - 


5.3.2 Model for Stator Slot Effects 

The stator slots cause variations in the magnetic field which produce losses in the 
retaining sleeve and magnets of the rotor. Accurate calculation of the losses in the retaining 
sleeve is extremely difficult. Several different methods have been developed and in this paper, 
the technique from reference [18] is employed. Figure 34 shows the geometry of the retaining 


can and the currents that are induced. 


Induced 
Currents 


Figure 34: Retaining Sleeve Induced Currents 
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As the rotor spins past the slot openings of the stator, the air gap flux density undergoes 
modulation due to the change in reluctance. The dip in B, (shown in Figure 35) travels along the 


B-waveform which is otherwise moving synchronously with the rotor. 


Stator 


Figure 35: Flux Density Variation 
The rotation of the rotor generates an E-field in the sleeve and a subsequent axial current 
density J = E/p. When this current density is integrated over the volume of the can, the average 


loss per unit area is determined (Eqn 5-26) [18]. 
2 [B.N(R +h Tt B 
E a [BN(R hy) where p = —. P 
3600 p V2 Jap 
Eqn 5-26 
It is evident from Eqn 5-26 that as the slot width increases, the width of the flux density dip (B) 
gets larger causing the sleeve losses to increase. The above equation only considers eddy 


currents flowing in the axial direction but there are also circumferential components. These 


portions are accounted for using a factor K, as shown in Eqn 5-27 where the total can losses are 


determined [18]. 
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Eqn 5-27 


One way to reduce the retaining can losses is to split the sleeve cylinder into separate 
rings. This reduces the effective flow path length for the eddy currents thus lowering the losses. 
There is a limit to how much this can be done since as the number of rings increases, the 
construction becomes more difficult. In this paper it is assumed that the maximum feasible 
number of sleeve sections is ten. 

When metal alloys are used as the retaining sleeve material, they effectively shield the 
magnets from the flux density variations, whereas there is little or no shielding with the carbon 
fiber composite because of its low conductivity. Therefore, with a metal alloy can, the rotor 
losses occur largely in the can compared to the carbon fiber sleeve, where the losses mostly 
occur in the magnets. These magnet losses are calculated using methods similar to Eqn 5-26 and 
Eqn 5-27 assuming that the eddy currents flow in the top 10% of the magnet volume. Appendix 
J contains the MATLAB code that performs the detailed rotor loss calculations caused by the 


stator slot effects. 


54 Complete System Model & Design Procedure 


The complete model for the power generation module consists of the PM generator, the 
rectifier, the filtering/energy storage components, the buck converter, and the load. The program 
PSIM (Version 6.0, Powersim, Inc.) is utilized to construct the circuit model, perform detailed 
time series and waveform analysis, and determine voltage and current values. An example of a 
7-phase version of this model system is shown in Figure 36. In each of the models, the converter 
and rectifier are simulated as one IGBT/diode module for simplicity but for actual 


implementation/construction, 19 modules are required as described in section 4.5. 
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GENERATOR RECT IFIER INPUT FILTER CONVERTER : OUTPUT: FILTERS OUTPUT- 


Figure 36: System Model 


The procedure for designing and analyzing the power generation module consists of the 
following: 

e Design a PM generator using the methodology from sections 3.4.2 and 5.2, iterating to 
ensure the sizes, current density, magnetic saturation, and sleeve thickness/stresses are all 
satisfactory 

e Verify the PM generator has proper output waveforms and THD per section 5.1 

e Calculate the rotor losses caused by stator slot harmonics as outlined in section 5.3.2 

e Determine the rectifier circuit parameters utilizing the procedures from section 4.2 

e Design the buck converter and associated input and output filters employing the methods 
of section 4.3 

e Construct the PSIM circuit using the correct generator parameters (number of phases, 
back voltage (Ea), and synchronous inductance (L,)), rectifier values, converter 
parameters, and filtering components 

e Simulate and analyze the PSIM model verifying the system meets the required 
specifications for power and output voltage 

e Determine the generator line harmonic current magnitudes from the PSIM model 

e Using the harmonic current values, calculate the rotor losses caused by the winding time 
and space harmonics as detailed in section 5.3.1 

e Calculate the switching losses and conduction losses of the power conversion module per 
section 4.4 


e Determine the overall losses and efficiency of the system 
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e Calculate the masses and volumes of the generator and power conversion module using 


the methodology of sections 3.4.2 and 4.5 


5.5 Optimization 

Optimizing the overall power generation module involves designing numerous PM 
generators and associated power conversion electronics to achieve an ideal result. Using the 
design procedure from section 5.4, machine parameters are varied producing different generator 
designs. For each unique generator, a power conversion module is devised so that its values are 
matched with the output of the generator. This results in a wide range of power generation 
modules. 

The main machine parameters that are varied are the number of poles, number of phases, 
and the retaining sleeve material as shown in Figure 37. This results in 80 different machines 
and associated power conversion modules. For each of these generators, the magnet thickness 
(hm), air gap (g), magnet angle (Om), and stack length (Ls) are iterated to ensure the sizes, current 
density, magnetic saturation, sleeve thickness/stresses, output waveforms, and THD are all 
satisfactory. 

Rotational speed (13,000 RPM) and magnet skew angle (10°) are held constant. The 
number of slots (N) is set to 36 for the 3 and 5-phase machines and to 72 for the 7, 9, and 11- 
phase machines to properly fit the windings and ensure a reasonably sinusoidal back voltage 
waveform. Other parameters including the number of slots/pole/phase (m), number of armature 
turns (Na), back voltage (Ea), synchronous inductance (Ls), electrical frequency (f), breadth 
factor (ky), pitch factor (kp), skew factor (ks), magnet factor (kg), and air gap flux density (By) all 


change as a result of the machine variations. 
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Number of Phases 
(3, 5, 7, 9, or 11) 


\ Kriba Titanium Inconel Carbon Fiber 
Sleeve MON ene SE 
Kim 


Figure 37: Machine Optimization Parameters 


Each machine is optimized to have minimal back voltage, size, and mass while ensuring 
it remains within the current density and magnetic saturation limits. The power conversion 
modules are then developed matching component ratings to those of the PM generators. The 
power electronics are designed to produce maximum power factor and minimal line harmonics 
as discussed in section 4.3. This yields the lowest possible losses in the PM generator thus 


maximizing efficiency and minimizing rotor heat generation. 
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Chapter 6 Results and Analysis 


6.1 General 


As discussed in section 5.5, the main machine parameters that are varied are the number 
of poles, number of phases, and the retaining sleeve material and resulting in 80 different 
machines and associated power conversion modules. These PM machines and power 
conversions modules are developed with varying results. Appendix K and Appendix L contain 
the detailed specifications for each of the generators and power conversion modules. Appendix 
M and Appendix N have the in-depth losses and weights for each design. Appendix O lists the 
mass and volume results for the rectifiers and Appendix P contains the mass and volume data for 


the converters. Å summary of the range of values for the designs is contained in Table 15 and 


Table 16. 


Table 15: General Module Specifications 


Back Synch 
Voltage E, | Inductance 
(RMS) L, (mH) 


Gen Volt 
THD 


"Minimum | 1,136.0 | 0.020 | 3.12% | 0.850 
3,380.0 | 0.348 | 23.05% | 1.080 
Average | 1,977.0 | 0.110 | 11.61% | 0.943 


Table 16: General Module Results 


Machine Total Machine PCM Total Line 
Weight Weight | Losses | Losses | Losses | Current 
ko kg kW kW kW THD 
4 


| Minimum | 1,388.6 | 2,302.7 | 3,900.1 1,688.3 | 2,126.7 | 6.8096 
| Maximum | 2,302.9 | 3,954.7 | 6,101.8 | 1,259.8 | 3,293.7 | 3,535.4 | 24.19% 
1,754.6 | 3,025.5 | 4,780.2 2,124.0 | 2,587.8 | 14.56% 


In general, the PCM losses are too high due to the IGBT switching losses. This is 
because high-frequency hard switching is used in the power conversion module as opposed to 
soft switching or a lower frequency. Soft switching involves using snubber circuits or resonant 
converters to minimize device switching losses by constraining the switching of the power 


devices to time intervals when the voltage across the device or the current through it is nearly 
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zero. Soft switching can reduce switching losses as much as 50-80% [67], [68]. In addition to 
lowering switching losses, soft switching reduces EMI, permits higher switching frequencies, 
and achieves higher efficiencies. Analyzing different soft-switching topologies is beyond the 
scope of this paper but the effects soft switching has on the PCM are studied. Using a lower 
switching frequency is not examined because this would produce a larger converter since the 
passive components would have to be bigger. 

In order to select a final design, the different variants are examined to ascertain which 


one is optimal. Parameters such as weights, machine losses, PCM losses, and THD are used to 


help facilitate the selection process. 


62 Number of Phases 

To facilitate studying the effects the number of phases has on the power generation 
module, several parameters are held constant. The retaining sleeve material and number of poles 
are arbitrarily set to inconel and 12. This does not affect the analysis because the trends related 
to varying the number of phases are consistent across all 80 power modules as indicated in the 


detailed results contained in the appendices. 
As the number of phases varies from 3 to 11, the weights of the PM generator and the 


power conversion module change, but only slightly. Figure 38 shows this and there is a 
maximum variation of approximately 242 kg, or 5.7%. The 11-phase machine produces the 


lowest weight power module but only by a small amount. 


Weights vs. Phase Number 


m 3 Phase 
m5 Phase 
07 Phase 
o9 Phase 

m 11 Phase 


Weight (kg) 


Machine PCM Total 


Figure 38: Weights vs. Number of Phases 
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The machine losses fluctuate significantly as the phase number changes with the 3-phase 
machine having 48.846 of the losses of the 11 phase machine (see Figure 39). The core, 
windage, and armature losses are approximately constant for the different phase generators, but 
the rotor losses vary considerably. This disparity is caused predominately by the AC line current 
harmonics from the rectifier being much greater in the higher phase machines. Their triple-n 
harmonics are higher order thus allowing lower order harmonics to have a greater effect, 
specifically the 3? harmonie: 

The higher rotor losses are also caused by the geometry of the machine. As the phase 
number increases, the number of slots/pole/phase (m) decreases which causes the winding 
factors to have less of an effect in reducing higher harmonic effects. The exception is the jump 
from 5-phases to 7-phases because the number of slots is doubled from 36 to 72 so that m 


increases initially. Overall, the 3-phase generator has the lowest rotor and total machine losses. 


Machine Losses vs. Phase Number 


3 Phase 
m5 Phase 


O17 Phase 
o 9 Phase 


Losses (kW) 


m 11 Phase 


Core Windage Armature Rotor Total | 


Figure 39: Machine Losses vs. Number of Phases 
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PCM Losses vs. Phase Number 
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Figure 40: PCM Losses vs. Number of Phases 


The PCM losses also differ greatly depending on the number of phases (Figure 40). The 
11-phase machine has the lowest total losses due to it having the lowest IGBT switching losses. 
This is because it has a smaller rectified voltage which reduces the stresses and switching losses 
in the IGBTs (see section 4.4). However, if soft switching is implemented, the PCM losses are 
similar for the different phase machines varying by a maximum of 63.5 kW instead of the 843 
kW associated with hard switching. Therefore, the lower phase machines become more 
attractive if soft switching is utilized. 

Lastly, as the number of phases increases, Figure 41 shows that in general, the generator 
voltage and line current THD increase. This is caused by the same effects as with the rotor 
losses, the slots/pole/phase and the AC line harmonics from the rectifier. As mentioned in 
section 5.1, the specification from MIL-STD-1399 for voltage THD is normally 5% so using 
these standards as guidelines, only the 3-phase and the 7-phase power modules are satisfactory. 
The line current THD should be made as low as possible because this directly affects the amount 
of losses and heating in the rotor. The typical limit for line current THD is 10% making the 3- 


phase generator the only one with acceptable THD performance [69]. 
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THD vs. Phase Number 
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Figure 41: THD vs. Number of Phases 


6.3 Retaining Material 


Similar to the analysis of the number of phases, several parameters are held constant 
when examining the various retaining sleeve materials. Specifically, the number of phases is 
arbitrarily set to 3 and the number of poles is fixed at 6. This does not affect the analysis 
because the trends related to varying the retaining material are consistent across the continuum of 
power modules as shown in the detailed results contained in the appendices. 

With respect to the weights of the PM generator and the PCM, the type of retaining 
sleeve has little effect. Using a stainless steel retaining can produces the largest power module 
but it is not appreciably larger than the others. When the retaining sleeve is constructed from 
titanium, inconel, or carbon fiber, the difference in total weight is less than 0.11%. Figure 42 
shows the effects the retaining sleeve material has on the weights of the power generation 


module. 
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Weights vs. Retaining Sleeve Material 
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Figure 42: Weights vs. Retaining Sleeve Material 


Machine Losses vs. Retaining Sleeve Material 
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Figure 43: Machine Losses vs. Retaining Sleeve Material 


The type of retaining sleeve material does affect the losses in the PM generator. The 
core, windage, and armature losses vary only slightly but the rotor losses fluctuate a great deal. 
Inconel has the lowest rotor losses (77.3 kW) while stainless steel has the highest (107.1 kW). 
The carbon fiber has the smallest losses in the actual can but larger losses in the magnets. This is 


because the carbon fiber provides no shielding for the magnets from flux variations and winding 
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current harmonics. Overall, the inconel sleeve produces the lowest rotor and total machine 
losses as shown in Figure 43. | 

The PCM losses are not influenced considerably by the kind of retaining can material 
(Figure 44). Depending on the type, the losses vary by about 157.2 kW, or 4.2%. When soft 
switching is utilized, the losses fluctuate by only 23.7 kW or 2.5%. Therefore, the PCM losses 


are not a major factor in selecting the retaining can material. 


PCM Losses vs. Retaining Sleeve Material " 
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Figure 44: PCM Losses vs. Retaining Sleeve Material 


Lastly, the selection of a retaining sleeve material is not a factor in the generator voltage 
THD or the AC line current THD. The voltage THD varies by only 0.546 and the line current by 
less than 0.1596 depending on the type of material as indicated in Figure 45. It is evident, 
however, that selected a 3-phase, 6-pole, 36-slot machine would not be wise since the voltage 


THD is above the 5% limit in all cases. 
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THD vs. Retaining Sleeve Material 
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Figure 45: THD vs. Retaining Sleeve Material 


64 Number of Poles 


The number of poles in a machine is an important factor because it affects numerous 
parameters including the electrical frequency, the magnet pole pitch, and the magnetic air gap 
flux density. To facilitate studying the effects of the number of poles, several parameters are 
held constant. The retaining sleeve material and number of phases are arbitrarily set to inconel 
and 3. This does not affect the analysis because the trends related to varying the number of poles 
are consistent across the range of machines as indicated in the detailed results in the appendices. 

First, the weights of the generator and PCM are examined to ascertain the effects of the 
number of poles. As Figure 46 indicates, as the number of poles increases the weights of both 
the machine and the PCM decrease. The generator weight goes down because assuming the 
same flux density and circumferential arc, a greater number of poles produces the same radial 
flux but requires less stator core back iron (see Figure 5). The PCM weight decreases because 
the higher frequencies and lower rectified voltage associated with the increased pole numbers 
reduce the size and energy requirements of the passive components, especially the capacitors. 
Therefore, the total power generation module weight decreases as the number of poles increases 


making the 12-pole machine the most attractive in terms of weight. 
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Weights vs. Pole Number 
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Figure 46: Weights vs. Number of Poles 


Next, looking at Figure 47, as the number of poles goes up, the losses in the machine 
increase, some more than others. The core losses increase because the higher frequencies cause 
the eddy currents and hysteresis to have a greater effect. However, the range in core losses is 
only 7 kW so this increase is not a major contributor. The windage losses are not significantly 


affected by the number of poles differing by only 3 kW across the pole variation. 


Machine Losses vs. Pole Number 
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Figure 47: Machine Losses vs. Number of Poles 


The armature losses do fluctuate some as the pole number changes with the 6-pole 


machine having the lowest (117.6 kW) and the 12-pole machine having the highest (137.9 kW). 
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This is due to the higher pole machines having higher stator current densities as the machine 
performance is optimized. The rotor losses also increase as the number of poles gets higher 
mainly because of the winding time and space harmonics. As the pole number increases, the 
number of slots/pole/phase (m) decreases which causes the winding factors to have less of an 
influence in reducing higher harmonic effects. 

Overall, all of the different pole machines have less than 160 kW of rotor losses with the 
6-pole machine having the lowest. However, as mentioned in section 6.3, the generator voltage 
THD is more significant in this machine causing the higher pole machines to become viable 
options. 

The number of poles also produces sizeable effects on the PCM losses. As the pole 
number increases, the PCM losses decrease due to much lower IGBT switching losses. This is 
because the higher pole machines have lower rectified voltages which reduce the stresses and 
switching losses in the IGBTs (see section 4.4). The difference between having 6-poles and 12- 
poles is 754.3 kW in PCM losses, a 22.9% reduction. If soft switching is implemented, the PCM 
losses are still lower for the higher pole machines, but the difference is now 98.1 kW, a 10.3% 
reduction. Overall, as shown in Figure 48, the 12-pole machine is the best option when I 


considering PCM losses. 
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Figure 48: PCM Losses vs. Number of Poles 


The generator voltage THD and the AC line current THD are both affected by the 


number of poles. Examining Figure 49, the 10-pole and 12-pole generators are the only variants 
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that produce voltage waveforms that have acceptable THD levels, 3.74% and 4.95% 
respectively. The AC line current THD values are satisfactory for all of the pole numbers with 
the 12-pole machine having the lowest value (6.82%). If the THD is examined in conjunction 
with the machine losses, the12-pole machine is the best option. It has the top overall THD 
performance, has 25.8 kW lower rotor losses compared to the 10-pole generator, and has 14.6 


kW lower total machine losses compared to the 10-pole design. 
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Figure 49: THD vs. Number of Poles 


Percent 


6.5 Final Power Module 


6.5.1 PM Generator 


Based on the previous analysis and the goal of producing an optimal power generation 
module, a final design is selected. The sizing codes and analysis from sections 3.4.1, 3.4.2, 
3.4.3, and 3.4.4 are also utilized in this process. The final PM generator is a 3-phase, 12-pole 
machine with an inconel retaining sleeve. The final design is not a fully optimized in terms of all 
parameters but represents the best selection using the methodology of this paper. Detailed 


characteristics for the PM generator are listed in Table 17. 
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Table 17: PM Generator Final Design Parameters 


0.848 m 
Armature Resistance 
Tip Speed 
[Phases [3 Frequency —  — | 13008z | 


WEM sama : quor Tm 
Back Iron Flux Densit 
[Slots/pole/phase(m) | ico, fT ET 


Magnet Height 


1.20 T 
0.961 


19.9 kW 137.9 kW 
344kW T Roer Loss — | 1589kW 


: Š anh: 
s 


g 522.13 kg 
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Toad] d L67375kg | GP 


A diagram of the PM generator is shown in Figure 50 with the stator winding layout 


given in Table 18. The numbers in the winding layout table refer to the slot number locations of 


Figure 50. 
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16 MW PM Generatór 


36 Slot, 12 Pole: 
13,000 RPM 


Figure 50: Diagram of PM Generator Final Design 


Table 18: Winding Layout 


6.5.2 Power Electronics Module 


In conjunction with selecting the final PM generator design, the associated power 
electronics module is determined. Figure 51 shows the complete power module, including 


individual component values. 
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Figure 51: Power Module Final Design Diagram 


Along with Figure 51, additional characteristics for the final power module are given in Table 19 


with detailed data contained in the appendices. 


Table 19: PCM Final Design Parameters 


r, m 


IGBT Switching Losses T 2,103.6 kW Diode Conduction Losses 168.2 kW 
IGBT Conduction Losses 267.6 kW Total PCM Losses 2,539.4 kW 
0.841 FEN STED 


e S TuS 


IGBT Switching Losses 420.7 kW Diode Conduction Losses 1 68.2 kW 
IGBT Conduction Losses 267.6 kW Total PCM Losses 856.5 kW 


CET RU GE Pee 


_ 3962 kg 420.9 ke 

IGBTs 51.8 kg 

| 1,159.2 kg 
NENNEN 


oss u 
2 


Structure 


JA LE 
Max Rectified Voltage 5,567 V 
Duty Cycle Output Voltage Ripple 
| (@16MW load) 9:1024 (@16MW load) sd 
Duty Cycle Output Voltage Ripple 
(& 100kW load) 01277 — | (2100 kW load) Wal 
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6.5.3 Performance/Waveforms 

As discussed in section 5.4, all 80 power modules are constructed and simulated using 
PSIM. In order to conserve space, only the simulation results for the final design are presented 
here. The module is simulated operating at the maximum load of 16 MW and at the minimum 


load of 100 KW. 


For the 16 MW loading, Figure 52 shows the generator voltage waveforms. Since the 


THD is low, the curves fairly closely resemble sine waves. 


Figure 52: PM Generator Voltage Waveforms 


Figure 53 through Figure 55 show plots for various parameters for the complete module loaded 
at 16 MW. The output voltage fluctuates approximately 0.09 V (0.001396) which is within 
specification, the output current varies by less than 22 A (0.009640), and the AC line currents 


have low harmonic content. 
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Figure 55: AC Line Current (16 MW) 


Figure 56 through Figure 58 display the module performance at the minimum load of 100 
kW. The output voltage ripple is 1.4 V, higher than at full load but still not significant. The 
output current also has some variation (about 25 A) but this is acceptable at the lower power 
levels. The AC line current has much higher harmonic content because of notches that occur in 
the waveforms caused by the DC link capacitor charging and discharging at the low load. 
However, the magnitudes of the harmonics are insignificant since the magnitude of the 


fundamental is low and therefore the rotor losses and heating are not adversely affected. 
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Figure 57: Output Current (100 kW) 
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Figure 58: AC Line Current (100 kW) 
6.6 Comparison 


The premise for designing and employing a high-speed PM generator and power 
conversion module is to offer a light weight, power dense alternative to existing wound rotor 
machines. In Table 1 in section 1.3.2, parameters for some typical wound rotor machines are 
presented. In order to compare them with the new design power module, these machines are 


scaled to 16 MW using the general relationship of Eqn 6-1. 


4 
Pwr = Vol? 


Eqn 6-1 


Table 20 then provides an evaluation of the different machines and power modules. 


Table 20: Comparison of Machines/Modules 
hi W Teig 


2 


4 | 37 | 
0 | 36 | 
5 | 32 | 

6 | 28 | 
0 


3.7 
3.6 
3.2 
2.8 
1.6 


Ë $, 
PCM N/A 


NA | 20 | 166 | 
13,000 
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Examining Table 20, the new design offers significant reductions in both weight and 
volume. Specifically, the new total power module is 11.6 — 16.8% the volume and 8.5 — 11.7% 
the weight of the other machines. This translates to about a 7x reduction in volume and a 10x 
reduction in weight. These reductions can provide flexibility to naval architects since power, 


weight, and volume are integral parts of the design and construction processes. 


116 


rator Top Magnetic Generator Magn enerator DIY 


Chapter 7 Conclusions and Recommendations 


7.1 Design Lessons Learned 


Numerous lessons are learned from both the generator and power electronics design 
processes. The entire 16 MW power module is a complex system with many interrelated 


components that affect the overall performance. 


7.1.1 PM Generator 

The final generator is the product of examining multiple parameters to determine an 
optimal design. Various lessons are learned relating to different portions of the machine 
including the number of slots, the number of poles, the number of phases, the magnetic 


geometry, the windings and output voltage, and optimization methods. 


7.1.1.4 Number of Slots 

Initially, it was thought that the number of slots could be held constant over the range of 
machine designs. However, this is not possible because as the number of phases increases, 
additional slots are necessary to fit the greater number of windings. One of the most important 
parameters that can be used to help reduce the harmonics in a rotating machine is the number of 
slots/pole/phase (m). Therefore, the number of slots can be used to help optimize this parameter 
and improve the shape of the output voltage waveform and THD. 

The number of slots also influences the electrical loading of the machine. The number of 
slots directly affects the slot size which helps determine the current density in the machine. This 
is usually one of the limits on the machine's performance because of cooling. Overall, the 
number of slots is an important parameter that must be carefully selecting when designing a 


machine. 


7.1.1.2 Number of Poles 


Traditionally, the number of poles is driven by the system frequency and prime mover 
speed. With a high-speed generator rectified to a DC distribution system, the number of poles 
becomes an important design variable. Therefore, throughout the design progression, careful 


consideration should be given to the number of poles. This proved to be true throughout this 
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paper since the number of poles affects the electrical, magnetic, and structural performance, 
including the electrical frequency, the voltage waveforms, the magnetic flux, the magnet volume, 
the magnet pole pitch, the air gap size, and the stator back iron thickness. 

It is recommended that initially the designer select a range of acceptable output electrical 
frequencies which puts bounds on the allowable number of poles. From there, the number of 
poles can be varied to help determine the size, waveforms, and performance of the generator. In 
general, if the rotational speed is held constant as the number of poles increases: 

e The number of slots/pole/phase (m) changes affecting the output waveforms and THD 
e The weight of the machine decreases 


e The electrical frequency increases 


7.1.1.3 Number of Phases 

An odd number of phases is usually the best option for rotating machines because it 
produces balanced operation with limited mechanical fluctuations. Initially, it was assumed that 
a higher number of phases (greater than 3) would produce a more optimal power module. It is 
discovered that this is not necessarily the case. As the number of phases increases, m decreases 
which can adversely affect the voltage THD. Also, the AC line current harmonics are more 
substantial in increased phase machines because the triple-n harmonics are higher order as the 
phase number increases. However, a big advantage of having a higher number of phases is that 
it reduces the requirements for the PCM because of lower rectifier voltage ripple and lower 
phase voltages and currents. | 

It is concluded that the benefits of 3-phases and a higher number of phases can be 
simultaneously realized through a modified design. Multiple 3-phase sets of windings can be 
located equidistant around the stator and then connected to the PCM to achieve the desired 
results. An example of this would be three 3-phase sets of windings, one set at 0°, 120°, and 


2409; one at 45°, 165°, and 285"; and one at 90°, 210°, and 330°. 


7.1.1.4 Magnetic Geometry 

Designing the magnetic geometry in a PM generator is an intricate process. Careful 
consideration should be given to the type of retaining material, the magnet material, and the 
desired voltage. It is determined that the magnet height and air gap length greatly affect the 


overall machine performance, and if they are not properly selected, the generator does not 
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operated properly. A typical starting point in a design is to have the magnet height be 
- approximately 5 — 10x the length of the air gap [17], [18], [30]. This ensures that the magnetic 


flux is adequate while not saturating the core material. 


7.1.1.5 Windings and Output Voltage 

As discussed in section 3.2.7, wye series connected windings are used in this paper. This 
produces a generator output voltage that is too high compared to the desired PCM output voltage 
of 700 V. It also unnecessarily stresses the PCM devices and produces a lower than optimal duty 
cycle. It is concluded that connecting some of the generator windings in parallel or using a lower 
number of turns to produce a lower output voltage might be viable options to alleviate this 
problem. However, care must be taken to account for possible circulating currents in the parallel 


windings in the event of a fault. 


7.1.1.6 Optimization Methods 


When designing the PM generator, the goal was to produce an optimal design that met 
the required performance criteria while minimizing the rotor losses, the size, and the weight of 
the machine. Optimizing the design involves trying to simultaneously maximize both the 
electrical loading (current density) and magnetic loading (air gap flux density). It is concluded 
that it is extremely difficult to do both at the same time since there are competing factors such as 
saturation limits, cooling limits, and rotational stress limits. 

In this paper, each of the 80 machines is optimized by trying to maximize the current 
density (within the limits) while ensuring an adequate air gap flux density. This is done because 
it seems to produce the smallest, high-performance machines. It is therefore recommended in 
any PM generator design process that careful consideration be given to both the electrical and 


magnetic loading. 


7.1.2. Power Electronics 

The power electronics module (PEM) designed in conjunction with the PM generator 
uses one of many possible topologies. The combination of the passive rectifier and buck 
converter is a stable, well-proven design. It is optimized to match the generator by selecting the 
passive components properly. However, the design does have limitations because it uses hard- 


switching and large passive components. 
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7.1.2.1 Rectification 

Passive rectification is utilized in the PEM because of its simplicity and potential for 
lower input harmonics. It is concluded that the rectification scheme is satisfactory and couples 
well with the PM generator. However, it is recognized that a more detailed analysis of the 


rectification scheme could produce a better topology. 


7.1.2.2. Conversion and Switching 

As discussed in Chapter 6, one of the main problems with the power electronics module 
is that the losses are too high due to the IGBT switching losses. This is because hard switching 
is used in the power conversion module as opposed to soft switching. Soft switching minimizes 
device switching losses, reduces EMI, permits higher switching frequencies, and achieves higher 
efficiencies. 

Analyzing different soft-switching topologies is beyond the scope of this paper but the 
effects soft switching has on the PCM are studied. It is concluded that if soft switching is 
implemented in the PCM, the efficiency increases by 10 — 20% with the losses reduced by 
overall I MW. Given the power requirements of the power generation module, soft switching 
should be employed. It is therefore recommended that in any high-power system that utilizes 


power conversion modules, soft switching should be examined to ascertain its usefulness. 


7.1.2.3 Control 

One control scheme for providing voltage regulation for the power generation module is 
discussed in section 4.3.4. This control methodology utilizing duty ratio control with a PID 
controller is relatively simple to implement. This ease of control is possible because the buck 
converter is a straightforward topology. It is concluded that if a more complex, soft-switching 
design is implemented, more extensive control schemes would be required. 

Even though it is not discussed in detail, the PM generator provides for easy control 
because it is a synchronous machine. This is a distinct advantage over an induction generation 


which requires extremely complicated control algorithms. 


72 Power Generation Module 


Overall, 80 different 16 MW machines and power conversion modules are sized, 


designed, and analyzed with a final design selected. The permanent magnet generation module, 
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consisting of the generator and associated power electronics, provides an excellent alternative to 
traditional wound rotor machines for naval applications. It is concluded that the new design 
offers significant reductions in both weight and volume. Specifically, it is estimated that the PM 
generation module has a 7x reduction in volume and a 10x reduction in weight compared to 
similarly rated wound rotor systems. These reductions can provide flexibility to naval architects 
since power, weight, and volume are integral parts of the design and construction processes. 
However, further study is necessary to verify the 16 MW PM generation module's thermal, 


structural, and mechanical performance. 


73 Recommendations/Further Study 


This paper discussed the electrical and magnetic design and analysis of a permanent 
magnet generation module including sizing, detailed machine design and analysis, analytical 
models for different rotor loss mechanisms as well as with other machine losses, power 
electronics and conversion for connecting the high-speed generator to a DC distribution system, 
and in depth simulation of the complete system. Numerous design issues are addressed and 
several issues are raised about the potential improvements a PM generation system can offer. A 
proposed PM generation module design is presented along with a detailed design methodology. 

However, this paper does not full cover all aspects associated with the design and 
analysis of a PM generation module. In order to fully ascertain whether the proposed design is 
technically sound, further study is required in both the machine and power electronics areas. 
Specifically, with the respect to the generator, the following is recommended: 

e Using the calculated losses, perform a detailed thermal analysis of the machine including 
development of the required cooling systems 
e Conduct an in-depth structural and mechanical study of the PM generator, including 
vibration and acoustic analysis 
For the power electronics portions, the further examination is recommended, including: 
° Examine different rectification schemes to see what effects they would have on the 
whole system 
e Explore soft switching topologies and perform detailed analysis on how to implement 
them to ensure acceptable switching losses while matching performance with the 


generator 
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For the complete PM generation system, the following is recommended: 
e Using the developed PM generator and power electronics module, implement a control 
scheme and conduct in-depth analysis with the system connected to an entire distribution 


system under changing load conditions 
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Glossary 


a Winding Pitch 
Ob Bottom Surface Coefficient 
Om Magnet Pitch Coverage Coefficient 
Ot l Top Surface Coefficient 
p Flux Density Dip Width 
£b Flux Density Exponent 
Ef Frequency Exponent 
y Winding Electrical Angle 
À Flux Linkage or Wavelength (depending) 
ÀB Flux Density Dip Interval 
Às Slot Fill Fraction 
n Efficiency 
Vair Kinematic Viscosity of Air at 20°C (m/s) 
y Power Factor Angle 
p Resistivity 
Pair Air Density at 20° C (kg/m?) 
Ps Steel Density (kg/m?) 
fm Magnet Density (kg/m?) 
Pe Conductor Density (kg/m?) 
Ost Stator Conductivity (S/m) 
05 Magnet physical angle 
Oa. Magnet skew angle 
T Air Gap Shear Stress 
Ts Stator Tooth Width + Stator Slot Width 
Lo Free Space Permeability 
Hrec Recoil Permeability 
(0 Electrical Frequency (rad/sec) 
Om Mechanical Frequency (rad/sec) 
A Magnetic Vector Potential 
Aac Armature Conductor Area 
Ag Air Gap Area 
Am Magnet Area 
As Slot Cross-sectional Area 
Bo Base Flux Density for Core Losses 
Bi Fundamental Flux Density 
Bp Back Iron Flux Density 
Ba Flux Density Dip 
Baux Ideal Radial Flux Density 
B, Air Gap Flux Density 
B, Remnant Flux Density 
B, Tooth Flux Density 
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Specific Heat Capacity of Air (J/kg*C) 
Flux Concentration Factor 
Friction Coefficient 


Buck Controller Average Duty Cycle 
Stator Core Back Iron Thickness 
Duty Ratio 

Overall Machine Diameter 


Buck Controller Error Signal 
Back Voltage (EMF) 


Electrical Frequency (Hz) 
Base Frequency for Core Losses 


Air Gap 
Effective Air Gap 


Slot Depression Depth 
Magnet Height/Thickness 
Slot Depth 


Phase Current 
Average Converter Output Current 


Current Density 


Wave Number 

Winding Breadth/Distribution Factor 
Electrical Power Waveform Factor 
Magnet Factor 

Current Waveform Factor 
Winding Pitch Factor 

Winding Skew Factor 

Winding Factor 

Carter Coefficient 

Leakage Factor 

Reluctance Factor 

Retaining Can Loss Factor 
Surface Current Density 


Overall Machine Length 
Armature Conductor Length 
Air Gap Inductance 

Self Slot Leakage Inductance 
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AC Generator 


Top Magnetic Generator 


Mutual Slot Leakage Inductance 
End Turn Inductance 

Overall Machine Length 

Slot Leakage Inductance 
Synchronous Inductance 

Rotor Stack Length 


Slots per Pole per Phase 
Mass of Armature Conductor 
Mass of Core 

Mass of Core Back Iron 
Mass of Core Teeth 

Mass of Magnets 

Mass of Shaft 

Total Mass 


Rotational Speed (RPM) 

Number of Armature Turns 

Turns per Coil 

Number of Stator Slots 

Actual Coil Throw (in slots) 
Full-pitch Coil Throw (in slots) 
Number of Stator Slots Short Pitched 


Pole Pairs 

Power Factor 

Base Power for Core Losses 
Total Core Loss 

Core Iron Loss 

Retaining Can Loss 

Tooth Loss 

Input Power 

Friction & Windage Losses 
Permeance Coefficient 
Power Conversion Module 
Power Electronics Module 
Permeance 

Machine Power 


Number of phases 


Rotor Radius 

Inner Bound of Magnet — 
Outer Bound of Magnet 
Armature Resistance 
Core Inside Radius 
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Core Outside Radius 
Diode On-State Resistance 
IGBT On-State Resistance 
Inner Magnetic Boundary 
Outer Magnetic Boundary 
Reynold's Coefficient 


Stator Back Iron Ratio (Yoke Thickness/Rotor Radius) 


Retaining Sleeve Thickness 
Peripheral Tooth Fraction 
Total Harmonic Distortion 


Magnet Speed (m/s) 
Rotor Tip Speed (m/s) 
Terminal Voltage 


Retaining Can Loss per Unit Area 
Slot Depression Width 

Slot Average Width 

Slot Bottom Width 

Slot Top Width (closest to rotor) 
Tooth Width 


Synchronous Reactance 


Surface impedance 
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Appendix A. Detailed Power Requirements 


In order to properly size the PM generator being designed, a typical load list is developed 
for an IPS naval ship. For the propulsion portion, it is assumed that the ship contains a single - 
shaft requiring 30,000 SHP. Shaft horsepower (SHP) is the power measured in the shafting 
within the ship by a torsion meter as close as possible to the propeller or stern tube [10]. This 
SHP then translates to 32,552.1 BHP (24,275 kW) assuming a 96% efficient electric propulsion 
motor and a 96% efficient power control module for the motor. Brake horsepower (BHP) is the 
power required at the engines/generators. It is also assumed that the ship contains two equal- 
sized generators. Overall, Table 21 lists the systems and their associated connected (maximum) 


electrical power requirements. 


Table 21: Ship Connected Loads 
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Classic U.S. Navy design practices (DDS 310-1) determine generator sizes by first 
estimating the expected load during the preliminary design phases. Propulsion load is a function 
of the hull form, maximum speed, and drive technology. Ship service loads are a function of the 
ship's mission, crew size, and ship's size. Once this is complete, the estimates are adjusted for 
anticipated construction and life cycle growth. Propulsion load is not expected to change but 
ship service loads are given a 20% margin for construction growth and a 20% margin for life 
cycle growth [11]. The generators are then sized using Eqn A- 1, where Pgen is the power 


provided by one generator. 


_ (1.2: 1.2 MaxLoad ) 


8n "^  [09.(n- 1)] 
where n = number of generators 
Eqn A- 1 
Applying this method to the total electrical load from Table 21 produces two generators that are 


each capable of 47.9 MW. 

Since IPS ships require electrical power for ship loads and propulsion, efficient power 
management should be implemented. This means that power is shared among all systems so that 
at top speed, more power is diverted to the propulsion motor, and at lower speeds power is 
available for other uses. It is assumed that rated ship service power and full propulsion power 
are possible at the same time only if both PM generators are available. If one generator is lost, 
the rated ship service load is still provided with the remaining power sent to the propulsion 
motor. In addition, full power for both the main motor and auxiliary propulsion is not required 
simultaneously. Similar to the classical method, a 2076 margin for construction growth and a 
20% margin for life cycle growth are allotted. Based on these ideas, Eqn A- 2 is developed to 
size the generators, where Picu is the total generation capacity (both generators). 


_ (1.21.2 MaxNonPropulsionLoad ) + MainMotorLoad 
Protal = 
n-1 
where n = number of generators 
Eqn A- 2 
Utilizing this process produces a total generation capacity of 31.99 MW or two generators that 


are each capable of 16 MW. 
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A standard U.S. Navy ship spends 95% of its time at less than full speed [12]. Because 
of this, the classical method unnecessarily over-estimates the power requirements for the ship 
and results in generator sets that are too large. If the average demand for each system is also 
considered, then the disparity becomes even worse. Table 22 lists each system along with a 
demand factor which represents the average amount of power required for normal operation of 
the systems. Based on the 9.8 MW total average loads, the classical method generators would 
each be loaded to only 10.2 % capacity which is extremely inefficient. 

The improved generator sizing method does a better job of ensuring efficient operation 
while not excessively limiting the ship if one generator is damaged or down for maintenance. If 
only one of the 16 MW generators is available, the connected ship service load (5.4 MW) can be 
supplied while still having 10.6 MW for propulsion, giving 44 96 propulsion capability which is 


acceptable when operating in reduced status. 
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Table 22: Ship Average Loads 
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Appendix B. | MATLAB Code: Basic Sizing Method 


Jonathan Rucker, MIT Thesis 

% May 2005 

96 Program: pml basic 

% Program performs basic sizing and parameter calculations 
% for generators. 


o Fo Yo Fo Yo Yo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Yo Fo Fo % To % Yo Yo Fo Yo Yo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Go Fo Fo 
% Definition of variables 


% Name Variable 

% General variables f 

% Pwr Required power 

9b rpm Speed (RPM) 

% psi Power factor angle 

96 f Electrical frequency (Hz) 

% omega Electrical frequency (rad/sec) 

% vtip Tip speed (m/s) 

% LovD L/D ratio 

% stress Gap shear stress (psi) 

% Rotor variables 

96 R Rotor radius (m) 

PD > Rotor diameter (m) 

% Lst Rotor stack length (m) 

%p Number of pole pairs 

% Bg Expected air gap flux density (T) 
% Stator variables 

% Kz Surface current density (A/m) 

% Jz Current density (A/m2) 

% hs Slot height (m) : 
Fo Yo Fo Fo Vo Yo Fo Yo To Yo Yo Fo Fo Fo Fo Fo Fo Fo Vo Fo Fo Yo Yo Vo Fo Fo Yo Fo Fo Fo Yo Yo Fo Fo Yo Fo Fo Fo Yo Fo Fo Fo Fo Yo Fo Yo Fo Fo Po 
clear; 

% Constants & conversion factors 

hs = .015; % Assume slot depth of 15 mm 
lams = 0.5; % Assume slot fill fraction 
convi = 9.81; % 9.81 W per Nm/s 

conv2 = 703.0696; % 103.0696 N/m2 per psi 


Fo Fo Yo Fo % Fo o Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Fo Fo Fo No Fo Fo Yo Fo Fo Yo Vo Fo Fo Fo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Po 


% INPUTS 
Fo Yo Vo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Vo Fo To Fo Fo Fo Fo Yo Fo Vo Fo Fo Fo Vo Yo Fo To Fo Fo Fo Fo Fo Fo Vo Fo Fo Vo Fo To Fo Yo Fo Fo To Po % 


Pwr = 16e6; % Required power 
vtip = 200; % Max tip speed (m/s) 
LovD = 2.851; % Wound rotor usually 0.5-1.0, PM 1.0-3.0 


% Shear stress usually 1-10 psi small machines, 10-20 large liquid 
% liquid cooled machines 


stress = 15; 
p=3; % Pole pairs 
Bg = 0.8; % Tesla 


o Yo Yo Yo Yo Yo Yo Yo Fo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Yo Yo Fo Vo Yo Yo Fo To Fo Fo Fo Yo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo % Po 
% Calculations 
% Size 
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% Initially use Pwr = 2*pi*R*Lst*stress*vtip 
% Lst = 2*LovD*R 

hscm = hs*100; 

R = sqrt(Pwr/(2*pi*(LovD*2)*vtip*stress*conv1 *conv2)); 

D z2*R; 

Lst = LovD*D; 

% Speed j 
omega = (p*vtip)/R; 

f = omega/(2*pi); 

rpm = (60*fyp; 

% Current densities 

Kz = (stress*conv2)/(Bg*100); 

Ja = 10*Kz/(hscm*lams); 


% Output 

fprintf( Basic Machine Design Wn"); 

fprintf( Input Parameters:\n’); 

fprintfPower =%10.1f kW Shear Stress = %10.1f psi\n',Pwr/1e3,stress); 
fprintf('L/D Ratio = %10.2f Tip Speed < %10.1f m/s\n',LovD,vtip); 
fprintf('Pole Pairs = %10.1f Air Gap Bg — %10.1f T\n',p,Bg); 


fprintf('Qutput:An'); 
fprintf('Rotor Radius = %10.3fm Stack Length = %10.3f m\n',R,Lst); 
fprintf('Speed =%10.0f RPM Frequency = %10.1f Hz\n'rpm,f); 
fprintf(‘Ja = %10.2f A/cm2\n' Ja); 
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Appendix C. PM Machine Database 


Specific Machines (Motor or Generator - specify w/ M or G 
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[2 j| 2 J s | 4 |J 5 | 6 | 7 | 8 | 9 j| 10 | 11 | 12 | 
Type (M-Motor, G-Generato) | G | G | M | M | M | M J MI GI G j G | G | G | 
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ower factor 
otational Speed (rpm 
otor radius (cm 


v 


= 
M 
E] 


105 


Magnet thickness (mm 

ctive length (m 
Pole pairs 
Number of phases 
umber of Slots 
lots per pole per phase 
lots short pitched 
lot depth (mm 
lot depression depth (mm 
lot depression width (mm 
eripheral Tooth Fraction 
ack iron thickness (mm 
urns per coil 
oils per phase 
lot fill fraction 
agnet remant flux density (T 
agnet electrical angle (deg e 
ase frequency (Hz 

ss (MT 
olume (cubic m 


ii 
å 
8 


3 
m 

o 

MNS 

o 


E 
N 


e 
N 


1.128 


= 
= 
m 
_ 
a 


118. 


œ 


14 112. 71 


N 
ion 
= 
A 


ength (m 4.87 NEN 
eight (m [ 15 | 11 | 335 | 0229 | 53 | 54 | 0. | 0.463 | 0.155 | | 
Woltage  — — Í  leovoc| | jaevac| | 15000 VAG 42 VDC [800 VDC] ^ [i10 VAC[2200 VAC] 


L =—— 


Ls | 14 | 15 | 16 ]| 17 | 18 | 19 [| 20 | 21 | 2 | 23 | 24 | 

Type (M-Motor,G-Generator) | M | G | G | G | G | M | M J| M [ G Jj G | G | G j 
[Power (kW) 10 | 200 | 200 | 250 | 250 | 3800 | 5000 | 4300 [ | 1400 | 25 | 110 | 
Powerfacor JPP 
Rotational Speed (rpm) å 1265 | 3600 | 90000 | 18000 70000 

otor radius (cm 

ir gap (mm 

agnet thickness (mm 
ctive length (m 
ole pairs 
umber of phases 
umber of Slots 
lots per pole per phase 
lots short pitched 


o de 
Sl ISIS 
[^] [^] 


A 
= 
N 


= 
= 
a 
o ys 
Ë o 


m 
Ë 


pl 


ke 
m 
o 
e 
© 


io 
En 
nN 
o 


D 


1.1 


| 584 [| 
0.097 | 0.16 | 


o 
Ln 
A 


k 
e 
© 


eo 
"m 
° 
N N 


kk 
^ 


lot depression depth (mm 

lot depression width (mm 

eripheral Tooth Fraction 

ack iron thickness (mm 

urns per coil 

oils per phase 

lot fill fraction 

agnet remant flux density (T 

agnet electrical angle (deg e 
Base frequency (Hz 

ss (MT 
Radial | Radia! | Radia 
18VAC| | 730 VAC 


L 180 | 


o 
eT RETE 
m D 


N 
° 
o 
he) 
I3 
E 8 


o 
8 e 
3 


° 
N 


= 


37 0.10 
Radial 1 Radial | Radial 
iC} 705 VDC. 


I 


NENNEN 
L 
Radial 
500 voG 


a ` 


= 
< 
GIS 


Radia 


Radial a 
600 VDC] 600 


eo 
E 
M 
M nm 
D [oz] 
NISIN 


141 


AC Generator Top Magnetic Generator f Magnet Generator DIY 


Page Intentionally Left Blank 


142 


AC Generator Top Magnetic Generator Magnet Generator DIY 


Appendix D. | MATLAB Code: Sizing Method 1 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pmlinput 

% Program used as input file for pm1calc 

% All necessary input parameters entered here. 


clear; 
o Fo Fo To Fo Fo Fo Fo % Fo Yo Fo Yo Fo Fo Fo Fo Yo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Yo Fo Yo Fo Fo Fo Fo Fo Po Fo Yo Vo Wo 


% Definition & Entry of variables 
% General variables 


Pwr = 1666; % Required power (W) 

rpm = 13000; % Speed (RPM) 

psi = 0; % Power factor angle 

% Rotor variables 

R = 0.147; % Rotor radius (m) 

hm = 0.025; % Magnet thickness (m) 

Lst = 0.838; % Rotor stack length (m) 

p=3; % Number of pole pairs 

Br=1.2; % Magnet remnant flux density (T) 
thm = 50; % Magnet physical angle (deg) 
thsk = 10; % Magnet skew angle (actual deg) 
% Stator variables 

q-3; % Number of phases 

Ns = 36; % Number of slots 

Nsp = 1; % Number of slots short pitched 

g = .004; % Air gap (m) 

tfrac = 0.5; % Peripheral tooth fraction 

hs = .025; % Slot depth (m) 

hd = .0005; % Slot depression depth (m) 

wd = 1e-6; % Slot depression width (m) 

syrat = 0.7; % Stator back iron ratio (yoke thick/rotor radius) 
Ne = 1; % Turns per coil 

lams = 0.5; % Slot fill fraction 

sigst = 6.0e+7; % Stator winding conductivity 

% Densities 

rhos = 7700; % Steel density (kg/m3) 

thom = 7400; % Magnet density (kg/m3) 

rhoc = 8900; % Conductor density (kg/m3) 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pmlcalc 

% Program performs sizing and parameter calculations 

% for permanent magnet machines with surface magnets and 

% slotted stators. 
% Program developed from J.L. Kirtley script with permission 


% MUST RUN pmiinput PRIOR TO RUNNING pmicalc 
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o Fo Yo Yo Fo Fo Fo Fo Yo Fo Fo Fo Yo Fo Fo Fo Fo Yo Fo Fo Yo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo To Fo 
% Definition of variables 


% Name Variable 

% General variables 

% Pwr Required power (W) 

% rpm Speed (RPM) 

% psi Power factor angle 

% f Electrical frequency (Hz) 

% omega Electrical frequency (rad/sec) 
% vtip Tip speed (m/s) 

% lambda Flux linkage 

% Ea RMS Internal voltage (V) 

% Rotor variables 

% R Rotor radius (m) 

% hm Magnet thickness (m) 

% Lst Rotor stack length (m) 

% p Number of pole pairs 

% Br Magnet remnant flux density (T) 
% thm Magnet physical angle (deg) 
% thsk Magnet skew angle (actual deg) 
% Stator variables 

%q Number of phases 

% m Slots per pole per phase 

% Ns Number of slots 

% Nsp Number of slots short pitched 
% g Air gap (m) 

% ge Effective air gap (m) 

% tfrac Peripheral tooth fraction 

% hs Slot depth (m) 

% hd Slot depression depth (m) 

% wd Slot depression width (m) 

% syrat Stator back iron ratio (yoke thick/rotor radius) 
% Nc Turns per coil 

% lams Slot fill fraction 

% sigst Stator conductivity 

% Kc Carter coefficient 

% Loss Models 

% PO Base power for core losses 

% FO Base frequency for core loss 
% BO Base flux density 

% epsb Flux density exponent 

% epsf Frequency exponent 

% rhos Steel density 

% rhom Magnet density 

% rhoc Conductor density 


Fo Yo Fo o Fo Fo Fo Fo Yo Fo Yo Vo Yo Fo Yo Fo Fo Yo Fo Yo Yo To Fo Vo Fo Fo Yo Fo Yo To Yo Fo Fo Fo Fo To Fo Vo Fo Fo Yo Fo Fo Yo Fo Yo To Yo Fo Yo Fo Yo o %o 


% Constants to be used 


mu0 = 4*pi*1e-7; 
tol = 1e-2; 

cpair = 1005.7; 
rhoair = 1.205 ; 
nuair = 1.5e-5; 


% Free space permeability 

% Tolerance factor 

% Specific heat capacity of air (J/kg*C) 

% Density of air at 20 C (kg/m3) 

% Kinematic viscosity of air at 20 C (m2/s) 


PO = 36.79; % Base Power Losss, W/lb 
FO = 1000; % Base freuency, 60 Hz 
BO = 1.0; % Base flux density, 1.0 T 
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epsb = 2.12; 
epsf = 1.68; 


% Generate geometry of machine 


% Number of slots/pole/phase 

m = Ns/(2*p*q); 

% Number of armature turns (each slot has 2 half coils) 
Na = 2*p*m*Nc; 

95 Tooth width 

wt = 2*pi*(R+g+hm+hd)*tfrac/Ns; 

% Slot top width (at air gap) 

wst = 2*pi*(R+g+hm+hd)*(1-tfrac)/Ns; 
% Slot bottom width 

wsb = wst*(R+g+hd+hs)/(R+g+hm+hd); 
% Stator core back iron depth (as p increases, dc decreases) 
dc = syrat*R/p; 

% Full-pitch coil throw 

Nsfp = floor(Ns/(2*p)); 

% Actual coil throw 

Nsct = Nsfp - Nsp; 

% Estimate end turn length 

% End turn travel (one end) 

laz = pi*(R+g+hm+hd+0.5*hs)*Nsct/Ns; 
% End length (half coil) 

le2 = pi*laz; 

% End length (axial direction) 

lel = 2*le2/(2*pi); 


% Calculate electrical frequency & surface speed 
f = p*rpm/60; I 

omega = 2*pi*f; 

vtip = R*omega/p; 


% Winding & skew factors 

gama = 2*pi*p/Ns; 

alfa = pi*Nsct/Nsfp; 

kp = sin(pi/2)*sin(alfa/2); 

kb = sin(m*gama/2)/(m*sin(gama/2)); 

kw = kp*kb; 

ths = ((p*thsk)+1e-6)*(pi/180); % skew angle (elec rad) 
ks = sin(ths/2)/(ths/2); 


% Calculate magnetic gap factor 

Rs = R+hm+g; 

Ri=R; 

RI-R; 

R2 = R+hm; 

kg = ((Ri^(p-1)/(Rs^Q*p)-Ri^Q*p)))*((p/(p- D) *(R2^(p- D-R L^(p--1))... 
*(p*Rs^(2*p)/(p-1)) (R1^(1-p)-R2^(1-p))); 


% Calculate air gap magnetic flux density 
% Account for slots, reluctance, and leakage 


WS = (wst+wsb)/2; % Average slot width 
taus = ws + wt; % Width of slot and tooth 
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Ke = 1/(1-(1/((taus/ws)*((5*g/ws)+1)))); 


ge = Kc*g; 

Cphi = (p*thm)/180; % Flux concentration factor 
Kl = 0.95; % Leakage factor 

Kr = 1.05; % Reluctance factor 

murec = 1.05; % Recoil permeability 


PC = hm/(ge*Cphi); % Permeance coefficient 
Bg = ((KI*Cphi)/(1+(Kr*murec/PC)))*Br; 


% Calculate magnetic flux and internal voltage 
thmrad = thm*(pi/180); 

BI = (4/pi)*Bg*kg*sin(p*thmrad/2); 

lambda = 2*Rs*Lst*Na*kw*ks*B 1/p; 

Ea = omega*lambda/sqrt(2); % RMS back voltage 


% Calculation of inductances/reactances 


% Air-gap inductance 

Lag = (q/2)*(4/pi)*(mu0*Na^2*kw^2*Lst*Rs)/(p^2*(g--hm)); 
% Slot leakage inductance 

perm = mu0*((1/3)*(hs/wst) + hd/wst); 

Las = 2*p*Lst*perm*(4*Nc^2*(m-Nsp)42*Nsp*Nc^2); 

Lam = 2*p*Lst*Nsp*Nc^2*perm; 


ifq ==3 

Lslot = Las + 2*Lam*cos(2*pi/q); % 3 phase equation 
else 

Lslot = Las - 2*Lam*cos(2*pi/q); % multiple phases 
end 
% End-turn inductance (Hanselman) 
As = ws*hs; % Slot area 


Le = ((Nc*mu0*(taus)*Na*2)/2)*log(wt*sqrt(pi)/sqrt(2*As)); 
% Total inductance and reactance 

Ls = Lag+Lslot+Le; 

Xs = omega*Ls; 


% Lengths, Volumes, and Weights 


% Armature conductor length 
Lac = 2*Na*(Lst+2*le2); 

% Armature conductor area (assumes form wound) 
Aac = As*lams/(2*Nc); 

% Mass of armature conductor 
Mac = q*Lac*Aac*rhoc; 

% Overall machine length 
Lmach = Lst+2*le1; 

% Core inside radius 

Rei = R+hm+g+hd+hs; 

% Core outside radius 

Reco = Rei+dc; 

% Overall diameter 

Dmach = 2*Rco; 

% Core mass 


Mcb = rhos*pi*(Rco^2-Rci^2)*Lst; % Back iron 
Mot = rhos*Lst*(Ns*wt*hs+2*pi*R*hd-Ns*hd*wd); % Teeth 
Mc = Mcb + Met; 


% Magnet mass 
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Mm = 0.5*(p*thmrad)*((R--hm)^2-R^2)*Lst*rhom; 
% Shaft mass 

Ms = pi*R^2*Lst*rhos; 

% 1596 service fraction 

Mser = 0.15*(Mc+Ms+Mm+Mac); 

96 'Total mass 

Mtot = Mser+Mc+Ms+Mm+Mac; 


% Stator resistance 
Ra = Lac/(sigst*Aac); 


% Core Loss Calculations 


% Tooth Flux Density 

Bt - Bg/tfrac; 

% Back iron flux density (Hanselman) 

Bb = Bg*R/(p*do); 

% Core back iron loss 

Pcb = Mcb*PO*abs(Bb/BO)^epsb*abs(f/FO)^epsf; 
% Teeth Loss 

Pct = Mct*PO*abs(Bt/BO)^epsb*abs(f/FO)^epsf, 
96 'Total core loss 

Pc = Pcb + Pct; 


96 Start loop to determine terminal voltage and current 
notdone = 1; 
1-0; 
Ia = Pwr/(q*Ea); 
while notdone ==1 
i=i+]1; 
xa = Xs*Ja/Ea; 
% Conductor losses 
Pa = q*Ia^2*Ra; 


% Gap friction losses 

% Reynold's number in air gap 

omegam = omega/p; 

Rey = omegam*R*g/nuair; 

9b Friction coefficient 

Cf = .0725/Rey^.2; 

95 Windage losses 

Pwind = Cf*pi*rhoair*omegam^3*R^4*T st; 


% Get terminal voltage 


Va = sqrt(Ea'2-((Xs+Ra)*Ta*cos(psi))'2)-(Xs+Ra)*Ta*sin(psi) 


Ptemp = q*Va*Ia*cos(psi)-Pwind; 
error = Pwr/Ptemp; 
err(i) = error; 
if abs(error-1) < tol 
notdone = 0; 
else 
Ja = Ia*error; 
end 
end 
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% Remaining performance parameters 


% Current density 

Ja = Ia/Aac; 

% Power and efficiency 
Pin = Pwr+Pc+Pa+Pwind; 
eff = Pwr/Pin; 

pf = cos(psi); 


fprintf(‘pm1calc complete: Ready An'); 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pmloutput 

% Program outputs values from pmlcalc. 

% Program developed from J.L. Kirtley script with permission 


% MUST RUN pmlinput and pm1calc PRIOR TO RUNNING pmloutput 


% Variables for output display 
Pout = Pwr/1e3; 

Jao = Ja/1e4; 

Pco = Pc/1le3; 
Pwindo = Pwind/1e3; 
Pao = Pa/le3; 

wso = ws*1000; 

hso = hs*1000; 

wto = wt*1000; 

dco = dc*1000; 

Lso = Ls*1000; 

hmo = hm*1000; 

go = g*1000; 


% Output Section: 
fprintf('nPM Machine Design, Version 1: Surface Magnet, Slotted Stator\n'); 


fprintf('Machine Size:n^; 

fprintf('Machine Diameter = %8.3f m Machine Length = %8.3f mW';Dmach,Lmach); 
fprintf(Rotor radius =  %8.3f m Active length = %8.3f m\n',R,Lst); 
fprintf(Slot Avg Width= %8.3f mm Slot Height = %8.3f mm\n',wso,hso); 
fprintf(Back Iron Thick = %8.3f mm Tooth Width 2 %8.3f mm\n',dco,wto); 


fprintf(Machine Ratings:\n'); 


fprintf(Power Rating 2 %8.1fkW Speed = %8.0f RPM\n’, Pout,rpm); 
fprintf('Va (RMS) = %8.0f V Current = %8.1f A\n’, Va,Ia); 
fprintf('Ea (RMS) = %8.0f V Arm Resistance = %8.5f ohm\n',Ea,Ra); 


fprintf(Synch Reactance = %8.3f ohm Synch Induct = %8.3f mH\n',Xs,Lso); 
fprintf(Stator Cur Den = %8.1f A/cm2 Tip Speed = %8.0f m/sW', Jao,vtip); 


fprintf(Efficiency = %8.3f Power Factor = %8.3f\n’, eff,pf); 

fprintf(‘Phases = %8.0f Frequency - %8.1f Hz\n',q,f); 

fprintf('Stator Parameters:Vn"); 

fprintf'Number of Slots = %8.0f Num Arm Turns = %8.0f \n',Ns,Na); 
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%8.3f 
%8.2f T 
fprintf('Slots/pole/phase = %8.2f\n',m); 


fprintf('Breadth Factor = 
fprintf("Tooth Flux Den = 


fprintf('Rotor Parameters:\n’); 
fprintf(Magnet Height - —968.2f mm 
fprintf('Air gap = %8.2f mm 
fprintf(Magnet Remanence = %8.2f T 
fprintfMagnet Factor= %8.3f 


fprintf('Machine Losses:\n’); 
fprintf('Core Loss = %8.1f kW 
fprintf('Windage Loss = %8.1f kW 


fprintf('Machine Weights:\n’); 
fprintf(‘Core = %8.2f kg 
fprintf(‘Magnet = %8.2f kg 
fprintf(Services = %8.2f kg 


AC Generator 


Top Magnetic Generator 


Pitch Factor = %8.3f \n', kb,kp); 
Back Iron = %8.2f T\n', Bt,Bb); 


Magnet Angle = %8.1f degm\n',hmo,thm); 


Pole Pairs= %8.0f \n',go,p); 
Aig Gap Bg= %8.2f TW',Br,Bg); 
Skew Factor = %8.3f W',kg,ks); 


Armature Loss = %8.1f kW\n', Pco,Pao); 


Rotor Loss = TBD kW\n’, Pwindo); 

Shaft = %8.2f kgin',Mc, Ms); 

Armature= — 968.2f kg, Mm, Mac); 

Total = %8.2f kgW',Mser,Mtot); 
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Appendix E. MATLAB Code: Sizing Method 2 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pm2input 

% Program used as input file for pm2calc 

95 All necessary input parameters entered here. 


clear; 
o Fo Fo % Fo Fo Yo Fo Fo Vo Fo Fo Fo Fo Fo Yo Yo Fo Fo Fo Fo Fo Yo Yo Fo Fo Fo Fo Yo Vo Fo Fo Yo Fo Fo Fo Yo Fo Fo Fo Fo Yo Fo Fo Yo Yo Fo Yo Yo Yo Yo No Yo 


% Definition & Entry of variables 
% General variables 


Pwr = 16e6; % Required power (W) 

rpm = 13000; % Speed (RPM) 

psi = 0; % Power factor angle 

Bsat = 1.65; % Stator saturation flux density 
% Rotor variables 

vtip = 200; % Tip speed limit (m/s) 

pz3 % Number of pole pairs 

Br = 1.2; % Magnet remnant flux density (T) 
thsk = 10; % Magnet skew angle (elec deg) 
PC =5.74; % Permeance coefficient for magnets 
% Stator variables 

Ja = 2200; % Initial current density (A/cm2) 
q=3; % Number of phases 

m=2; % Slots/pole/phase 

Nsp = 1; % Number of slots short pitched 
g = .004; % Air gap (m) 

hs = .025; % Slot depth (m) 

hd = .0005; % Slot depression depth (m) 

wd = le-6; % Slot depression width (m) 

ws = .016; % Avg slot width (m) 

Nc = 1; % Turns per coil 

lams = 0.5; % Slot fill fraction 

sigst = 6.0e+7; % Stator winding conductivity 

% Densities 

rhos = 7700; % Steel density (kg/m3) 

rhom - 7400; % Magnet density (kg/m3) 

rhoc = 8900; % Conductor density (kg/m3) 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pm2calc 

% Program performs sizing and parameter calculations 

% for permanent magnet machines with surface magnets and 
% slotted stators. 


% MUST RUN pm2input PRIOR TO RUNNING pm2calc 
% Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Po 909596 


% Definition of variables 
% Name Variable 
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% General variables 
% Pwr 

% rpm 

% psi 

% f 

% omega 

% vtip 

% lambda 

% Ea 

% Rotor variables 
%R 


Required power (W) 

Speed (RPM) 

Power factor angle 

Electrical frequency (Hz) 
Electrical frequency (rad/sec) 
Tip speed (m/s) 

Flux linkage 

RMS Internal voltage (V) 


Rotor radius (m) 


% hm Magnet thickness (m) 

% Lst Rotor stack length (m) 

% p Number of pole pairs 

% Br Magnet remnant flux density (T) 
% thm Magnet physical angle (deg) 
% thsk Magnet skew angle (actual deg) 
% Stator variables 

% q Number of phases 

9b m Slots per pole per phase 

% Ns Number of slots 

% Nsp Number of slots short pitched 
% g Air gap (m) 

% ge Effective air gap (m) 

% tfrac Peripheral tooth fraction 

% hs Slot depth (m) 

% hd Slot depression depth (m) 

% wd Slot depression width (m) 

% syrat Stator back iron ratio (yoke thick/rotor radius) 
% Nc Turns per coil 

% lams Slot fill fraction 

% sigst Stator conductivity 

% Kc Carter coefficient 

96 Loss Models 

% PO Base power for core losses 

% FO Base frequency for core loss 
% BO Base flux density 

% epsb Flux density exponent 

% epsf Frequency exponent 

% rhos Steel density 

% rhom Magnet density 

% rhoc Conductor density 


Fo Yo Yo Fo Fo Fo Fo Fo Yo Fo Fo Fo Yo Fo Yo Fo Fo Fo Fo Fo Vo Fo Yo Vo Fo Fo Yo Fo Yo Fo Fo Fo Fo Fo % Fo Fo Yo Fo Fo Fo Fo Vo Yo Fo Yo Fo Fo Fo Yo Fo Yo Po Fo 


% Constants to be used 


mu0 = 4*pi*1e-7; 
tol = 1e-2; 

cpair = 1005.7; 
rhoair = 1.205 ; 
nuair = 1.5e-5; 
PO = 36.79; 


epsb = 2.12; 
epsf = 1.68; 


% Free space permeability 

% Tolerance factor 

% Specific heat capacity of air (J/kg*C) 

% Density of air at 20 C (kg/m3) 

% Kinematic viscosity of air at 20 C (m2/s) 
% Base Power Losss, W/Ib 

% Base freuency, 60 Hz 

% Base flux density, 1.0 T 
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% Calculate electrical frequency & rotor radius 
f = p*rpm/60; 

omega = 2*pi*f; 

R = p*vtip/omega; 


% Winding & skew factors 

Ns = floor(2*q*p*m); % Number of slots 
gama = 2*pi*p/Ns; 

Nsfp = floor(Ns/(2*p)); 

Nsct = Nsfp - Nsp; 

alfa = pi*Nsct/Nsfp; 

kp = sin(pi/2)*sin(alfa/2); 

kb = sin(m*gama/2)/(m*sin(gama/2)); 

kw = kp*kb; 

ths = ((p*thsk)+1e-6)*(pi/180); 76 skew angle (elec rad) 
ks = sin(ths/2)/(ths/2); 


% Calculate magnet dimensions, tooth width, & air gap flux density 


thme = 1; % Initial Magnet angle (deg e) 
notdone = 1; 
ge-g % Initial effective air gap 
while notdone == 1 
alpham = thme/180; % Pitch coverage coefficient 
Cphi = (2*alpham)/(1+alpham); % Flux concentration factor 
hm = ge*Cphi*PC; % Magnet height 
Ds = 2*(R+hm+g); % Inner stator/air gap diameter 
Kl = 0.95; % Leakage factor 
Kr = 1.05; % Reluctance factor 
murec = 1.05; % Recoil permeability 


Bg = ((KI*Cphi)/(1+(Kr*murec/PC)))*Br; 
wt = ((pi*Ds)/Ns)*(Bg/Bsat); % Tooth width 


taus = ws + wt; % Width of slot and tooth 
Kc = 1/(1-(1/((taus/ws)*((S*g/ws)+1)))); % Carter's coefficient 
ge = Kc*g; 


eratio = ws/wt; 
if abs(eratio - 1) « tol 
notdone = 0; 
else 
thme = thme + 1; 
end 
end 


% Set final values 


thm = thme/p; % Magnet physical angle 
thmrad = thm*(pi/180); 

hm = ge*Cphi*PC; % Magnet height 

Ds = 2*(R+hm+g); % Inner stator/air gap diameter 


% Generate geometry of machine 


% Peripheral tooth fraction 

tfrac = wt/(wt+ws); 

% Slot top width (at air gap) 

wst = 2*pi*(R+g+hm+hd)*tfrac/Ns; 

% Slot bottom width 

wsb = wst*(R+g+hm+hd+hs)/(R+g+hm-+hd); 
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E 


% Stator core back iron depth 

dc = (pi*Ds*thmrad/(4*p))*(Bg/Bsat); 
% Core inside radius 

Rci = R+hm+g+hd+hs; 

% Core outside radius 

Rco = Rci+dc; 

% Slot area 

As = ws*hs; 

% Estimate end turn length 

% End turn travel (one end) 

laz = pi*(R+g+hm+hd+0.5*hs)*Nsct/Ns; 
% End length (half coil) 

le2 = pi*laz; 

% End length (axial direction) 

lel = 2*le2/(2*pi); 


% Calculate magnetic gap factor 

Rs = R+hm+g; 

Ri=R; 

R1 =R; 

R2 = R+hm; 

kg = ((Ri^(p-1)/(Rs^(*p)-Ri^(2*p)) *(p/(p--1)) (R2^(p-1)-R1^(p*-1))... 
*(p*Rs^(2*p)/(p- 1))(R1^(1-p)-R2^(1-p))); 


26 Core loss calculations (per length) 


% Core mass per length 


McbperL = rhos*pi*(Rco^2-Rci^2); % Back iron 
MctperL = rhos*(Ns*wt*hs+2*pi*R*hd-Ns*hd*wd); % Teeth 
McperL = McbperL + MctperL; 

% Tooth Flux Density 

Bt = Bg/tfrac; 


95 Back iron flux density (Hanselman) 

Bb = Bg*R/(p*dc); 

% Core back iron loss per length 

PcbperL = McbperL*PO*abs(Bb/BO)^epsb*abs(f/FO)^epsf; 
% Teeth Loss per length 

PctperL = MctperL*PO*abs(BUBO)^epsb*abs(f/FO)^epsf; 
% Total core loss per length 

PcperL = PcbperL + PctperL; 


% Current and surface current density 

% Armature turns (each slot has 2 half coils) 
Na = 2*p*m*Nc; 

% Arm cond area (assumes form wound) 
Aac = (As*lams)/(2*Nc); 


% Power & Current waveform factors (Lipo) 
ke = 0.52; 
ki = sqrt(2); 


% Initial terminal current 
Ia = Ns*lams*As*Ja* 1e4/(2*q*Na); 


notfin = 1; 
Lst = 0.1; % Initial stack length 
i=]; 
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% Start loop to determine Lst, Ea, Va, and Ia 
notdone - 1; 
k=0; 
while notdone == 
k=k+1; 
% Surface current density 
A = 2*q*Na*Ia/(pi*Ds); 


% Calculate stack length of machine 
% Loop to get stack length 
while notfin == 
% Gap power 
Pgap = 4*pi*ke*ki*kw*ks*kg*sin(thmrad)*(f/p)*A*Bg*(DsA2)*Lst; 
% Length of conductor i 
Lac = 2*Na*(Lst+2*le2); 
% Stator resistance 
Ra = Lac/(sigst* Aac); 
% Copper Loss 
Pa = q*Ia*2*Ra; 
% Core losses 
Pc = PcperL*Lst; 


% Iterate to get length 
Ptemp1 = Pgap-Pa-Pc; 
error = Pwr/Ptemp1; 


err(i) = error; 

if abs(error-1) < tol 
notfin = 0; 

else 
Lst = Lst*error; 
i=i+1; 

end 

. end 


% Calculate magnetic flux and internal voltage 
thmrad = thm*(pi/180); 

B1 = (4/pi)*Bg*kg*sin(p*thmrad/2); 

lambda = 2*Rs*Lst*Na*kw*ks*B 1/p; 

Ea = omega*lambda/sqrt(2); 96 RMS back voltage 


% Calculation of inductances/reactances 


% Air-gap inductance 

Lag = (q/2)*(4/pi)*(mu0*Na^2*kw^2*Lst*Rs)/(p^2*(g--hm)); 
% Slot leakage inductance 

perm = mu0*((1/3)*(hs/wst) + hd/wst); 

Las = 2*p*Lst*perm*(4*Nc^2*(m-Nsp)42*Nsp*Nc^2); 

Lam = 2*p*Lst*Nsp*Nc^2*perm; 


if q == 

Lslot = Las + 2*Lam*cos(2*pi/q); % 3 phase equation 
else 

Lslot = Las - 2*Lam*cos(2*pi/q); % multiple phases 
end 
% End-turn inductance (Hanselman) 
taus = ws + wt; % Width of slot and tooth 


Le = ((Nc*mu0*(taus)*Na^2)/2)*log(wt*sqrt(pi)/sqrt(2* As)); 
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96 Total inductance and reactance 
Ls = Lag+Lslot+Le; 
Xs = omega*Ls; 


% Lengths, Volumes, and Weights 


% Armature conductor length 

Lac = 2*Na*(Lst+2*le2); 

% Mass of armature conductor 
Mac = q*Lac*Aac*rhoc; 

% Overall machine length 

Lmach = Lst+2*le1; 
.% Overall diameter 

Dmach = 2*Rco; 

% Core mass 

Mc = McperL*Lst; 

% Magnet mass 

Mm = 0.5*(p*thmrad)*((R--hm)^2-R^2)*Lst*rhom; 
9b Shaft mass 

Ms = pi*R^2*Lst*rhos; 

% 15% service fraction 

Mser = 0.15*(Mc+Ms+Mm+Mac); 
% Total mass 

Mtot = Mser+Mc+Ms+Mm+Mac; 


% Gap friction losses 

% Reynold's number in air gap 

omegam = omega/p; 

Rey = omegam*R*g/nuair; 

% Friction coefficient 

Cf = .0725/Rey^.2; 

% Windage losses 

Pwind = Cf*pi*rhoair*omegam^3*R^4*T st; 


% Get terminal voltage 
xa = Xs*Ia/Ea; 
Va = sqrt(Ea*2-((Xs+Ra)*Ia*cos(psi))“2)-(Xs+Ra)*Ia*sin(psi); 


Ptemp = q*Va*Ia*cos(psi)-Pwind; 
Perror = Pwr/Ptemp; 
Perr(k) = Perror; 
if abs(Perror-1) « tol 
notdone = 0; 
else 
Ia = Ia*Perror; 
end 
end 


% Remaining performance parameters 


% Current density 

Ja = Ia/Aac; 

% Power and efficiency 
Pin = Pwr+Pc+Pa+Pwind; 
eff = Pwr/Pin; 

pf = cos(psi); 


156 


AC Generator Lop Magnetic Generator T Generator DIY 


fprintf((pm2calc complete: Ready. n"); 


% Jonathan Rucker, MIT Thesis 

95 May 2005 

% Program: pm2output 

% Program outputs values from pm2calc. 

% Program developed from J.L. Kirtley script with permission 


% MUST RUN pn input and pm2calc PRIOR TO RUNNING pm2output 


% Variables for output display 
Pout = Pwr/1e3; 

Jao = Ja/le4; 

Pco = Pc/1e3; 
Pwindo = Pwind/1e3; 
Pao = Pa/le3; 

wso = ws*1000; 

hso = hs*1000; 

wto = wt*1000; 

dco = dc*1000; 

Lso = Ls*1000; 

hmo = hm*1000; 

go = g*1000; 


% Output Section: 
fprintf(\nPM Machine Design, Version 2: Surface Magnet, Slotted Stator\n’); 


fprintf(‘Machine Size:\n’); 

fprintf'Machine Diameter = %8.3f m Machine Length = %8.3f m\n'’,Dmach,Lmach); 
fprintf('Rotor radius - “%#8.3f m Active length = %8.3f m\n',R,Lst); 
fprintf('Slot Avg Width 2 %8.3f mm Slot Height = %8.3f mm\n',wso,hso); 
fprintf('Back Iron Thick = %8.3f mm Tooth Width= %8.3f mm\n',dco,wto); 


fprintf(‘Machine Ratings:\n’); 


fprintf(Power Rating <=  %8.1f kW Speed = %8.0f RPM\n', Pout,rpm); 
fprintf('Va (RMS) = %8.0f V Current = %8.1f A\n', Va,Ia); 
fprintf('Ea (RMS) = %8.0f V Arm Resistance = %8.5f ohm\n',Ea,Ra); 


fprintf(Synch Reactance = %8.3f ohm Synch Induct = %8.3f mH\n',Xs,Lso); 
fprintf('Stator Cur Den = %8.1f A/cm2 Tip Speed = %8.0f m/s\n’, Jao,vtip); 


fprintf('Efficiency = “%8.3f Power Factor = %8.3f\n', eff,pf); 
fprintf('Phases = %8.0f Frequency - %8.1f Hz\n',q,f); 
fprintf('Stator Paramieters:n"); 

fprintf'Number of Slots = %8.0f Num Arm Turns = %8.0f \n',Ns,Na); 
fprintf('Breadth Factor 2 %8.3f Pitch Factor = %8.3f W', kb,kp); 
fprintf('Tooth Flux Den =  768.2£ T Back Iron = %8.2f TW', Bt,Bb); 


fprintf('Slots/pole/phase = %8.2f\n',m); 


fprintf('Rotor Parameters:\n’); 

fprintf(Magnet Height <= %8.2f mm Magnet Angle = %8.1f degm\n',hmo,thm); 
fprintf( Air gap = %8.2f mm Pole Pairs - %8.0f \n',go,p); 
fprintf(Magnet Remanence = %8.2f T AigGapBg= %8.2f TW',Br,Bg); 
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-—————————————————————————, il 


fprintf'Magnet Factor = %8.3f 


fprintf('Machine Losses:\n’); 
fprintf('Core Loss = %8.1f kW 
fprintf('Windage Loss =  %8.1f kW 


fprintf(Machine Weights:\n’); 


fprintf(Core = %8.2 kg 
fprintf(‘Magnet = %8.2f kg 
fprintf('Services = %8.2f kg 


enerator 


Skew Factor =  968.3f W',kg,ks); 


Armature Loss = %8.1f kW\n', Pco,Pao); 


Rotor Loss = TBD kWW', Pwindo); 

Shaft= — %8.2f kg\n',Mc,Ms); 

Armature= “#3.2f kgw', Mm, Mac); 

Total = %8.2f kg\n',Mser,Mtot); 
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Appendix F. MATLAB Code: Bode Plot 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: Buckfilter 

% Program calculates transfer function and outputs 
% Bode plot for buck converter input filter 


clear; 


% Input parameters 
R=4.79e-3; 

Cf = 2.84e-3; 

Cb = 28.4e-3; 

Lf = 1.415e-6; 


% Set up transfer function 

num = [R*Cb 1]; 

den = [Lf*R*Cf*Cb Lf*(Cf+Cb) R*Cb 1]; 
H = tf(num,den); 

bode(H) 


AC Generator Top Magnetic Generator 
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Appendix G. | MATLAB Code: PM Generator Waveforms 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pmwave 

% Program calculates and outputs different waveforms, 

9b calculates THD, and computes the harmonic content 

% for permanent magnet machines with surface magnets and 
% slotted stators. 


% MUST RUN pmlinput and pm1calc PRIOR TO RUNNING pmwave 


o Fo Fo Fo Yo Yo Fo Fo Yo Yo Yo Fo Fo Mo Yo Fo Fo Vo Fo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Yo Yo Fo Yo Yo Fo Fo Yo Yo Yo Fo Fo Fo Yo Yo Vo Po 
% Definition of variables 


% Name Variable 

% General variables 

% Pwr Required power (W) 

% rpm Speed (RPM) 

% psi Power factor angle 

% f Electrical frequency (Hz) 

% omega Electrical frequency (rad/sec) 
% vtip Tip speed (m/s) 

% lambda Flux linkage 

% Ea RMS Internal voltage (V) 

% Rotor variables 

% R Rotor radius (m) 

% hm Magnet thickness (m) 

% Lst Rotor stack length (m) 

96 p Number of pole pairs 

% Br Magnet remnant flux density (T) 
% thm Magnet physical angle (deg) 
% thsk Magnet skew angle (actual deg) 
% Stator variables 

% q Number of phases 

% m Slots per pole per phase 

% Ns Number of slots 

% Nsp Number of slots short pitched 
% g Air gap (m) 

% ge Effective air gap (m) 

% tfrac Peripheral tooth fraction 

% hs Slot depth (m) 

% hd Slot depression depth (m) 

% wd Slot depression width (m) 

% syrat Stator back iron ratio (yoke thick/rotor radius) 
% Nc Turns per coil 

% lams Slot fill fraction 

% sigst Stator conductivity 

% Kc Carter coefficient 


o Fo Fo Fo Fo Fo Fo Fo % Fo Yo Vo Vo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Yo Fo Fo Fo Fo Fo Fo Yo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo To To PoFo ' 


% Constants to be used 
mu0 = 4*pi*1e-7; % Free space permeability 
tol = le-2; % Tolerance factor 
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95 Harmonics to be evaluated 


n= 1:2:35; 

np=p.*n; % Use in kgn equation 
w-n.*omega; % Harmonic angular frequencies 
freq = w ./ (2*pi); % Harmonic frequencies 


% Harmonic winding and skew factors 

gama = 2*pi*p/Ns; 

alfa = pi*Nsct/Nsfp; 

kpn = sin(n .* pi/2) .* sin(n .* alfa/2); 

kbn = sin(n .* m*gama/2) ./ (m*sin(n .* gama/2)); 

kwn = kpn .* kbn; 

ths = ((p*thsk)+1e-6)*(pi/180); % skew angle (elec rad) 
ksn = sin(n .* ths/2) ./ (n .* ths/2); 


% Calculate magnetic gap factor 

Rs = R+hm+g; 

Ri=R; 

R1=R; 

R2 = R+hm; 

kgn = ((Ri.^(np-1))./(Rs.^(2.*np)-Ri.^(2.*np))).*((np-/(np--1))... 
*(R2.A(np+1)-R1.A(np+1))+(np.*Rs.4(2.*np)./(np-1))... 
A(RI.^(I1-np)-R2.^(1-np))); 


% Calculate magnetic flux and internal voltage 

thmrad = thm*(pi/180); 

thmerad = p*thmrad; 

Bn = Bg.*((4/pi)./n).*kgn.*sin(n.*thmerad/2).*sin(n.*pi/2); 
lambdan = ((2*Rs*Lst*Na).*kwn.*ksn.*Bn)./p; 

Ean = (omega.*lambdan); % Peak back voltage 


% Normalized values for plotting 
Eanorm = abs(Ean) ./ Ean(1); 


% Voltage THD 
Eah = 0; 
for r = 2:length(n) 
Eah = Eah + Ean(r)^2; 
end 
THD = 100*sqrt(Eah/(Ean(1)^2)); 


% Generate waveforms 
% Rotor physical angle goes from 0 to 2*pi - electrical to 2*p*pi 
ang = 0:pi/100:2*pi; 
angp = p*ang; 
Bout = zeros(size(angp)); 
Eaout = zeros(size(angp)); 
for i= l:length(n) 
Bout = Bout + Bn(i).*sin(n(i).*angp); 
Eaout = Eaout + Ean(i).*sin(n(i).*angp); 
end 


% Plot waveforms 
figure(1) 


plot(ang,Bout); 
title PM Generator: Flux Density"); 


162 


AC Generator Top Magnetic Generator Magnet Generator DIY 


ylabel( B (Tesla); 
xlabel('Rotor Angle (rad); 


figure(2) 

plot(ang,Eaout); f 
title((PM Generator: Back EMF'); 
ylabel(Peak Voltage (V)); 
xlabel('Rotor Angle (rad); 


figure(3) 
hold on 
title(('PM Generator: EMF Harmonics, THD = ',num2str( THD),' %']); 
ylabel(‘Normalized Back EMF’); l 
xlabel('Harmonic Number"); 
text(20,0.7,'Dark: Above 10% of Fundamental','FontSize',10); 
text(20,0.65, Light: Below 10% of Fundamental','FontSize',10); 
for z = 1:length(Eanorm) 
if Eanorm(z) « 0.10 
bar(n(z),Eanorm(z),'c'); 
else 
bar(n(z),Eanorm(z),'b'); 
end 
end 
hold off 
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Appendix H. MATLAB Code: Retaining Sleeve Stress Calculations 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pmcanstress 

% Program calculates and outputs retaining can stress 

% for permanent magnet machines with surface magnets and 
% slotted stators. 


% MUST RUN pmlinput, pmlealc, and pmwave PRIOR TO RUNNING pmcanstress 


% Fo % Fo Yo Fo Yo Fo Yo Fo Yo Fo Fo Vo Yo Fo Fo Yo Fo Fo Fo Fo Yo Fo Yo Fo Fo Fo Fo Fo Fo Yo Fo Fo Yo Vo Fo Yo Yo Fo Fo Yo Fo Fo Yo Fo Fo Yo Fo Fo Yo Fo Fo To Po 
% Calculate retaining sleeve stresses 

Fo Fo Fo Fo Yo Yo Fo Fo Fo Fo Yo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Yo Fo % Fo Yo Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo To Fo Yo Po 
% Conversion 

Patopsi = 1.45038e-4; % psi per Pa 

% Material yield stresses (ksi) 

Stain_str = 90; 

Alum, str = 75; 

Titan str = 110; 

CarFib. str = 100; 

Inconel str = 132; 

% Safety factor 

SF = 1.2; 


% Force on magnets/sleeves is centrifugal force 
% Magnet tangential velocity 

vmag = ((R+hm)*omega)/p; 

% Centrifugal force 

Fm = (Mm*vmag^2)/(R--hm); 

% Outward pressure 

Phoop = Fm/(2*pi*(R--hm)*Lst); 

% Hoop Stress (in general, str = P*R/t) 


stop = 22; 
for i = l:stop 
t(1) = i*.0005; % sleeve thickness t 


slev(i) = t(i)*1000; 
Sthoop(i) = (Phoop*(R--hm)/t(1))*Patopsi/1000; 
SFHoop(i) = Sthoop(i)*SF; 

end 


% Output results 

fprintf(‘Retaining Sleeve Stress:\n’); 

fprintf('Stress Limits:\n’); 

fprintf('Stainless Steel = %6.1f ksi Aluminum Alloy = %6.1f ksiln',Stain str, Alum str); 
fprintf(‘Titanium Alloy = %6.1f ksi Carbon Fiber = %6.1f ksi\n',Titan_str,CarFib_str); 
fprintf('Inconel = %6.1f ksi\n',Inconel_str); 

fprintf('Actual Sleeve Stress:\n'); 

fprintfSleeve Thickness Actual Stress SF Stress); 


for i= 1:stop 
fprintf(' %5.2f mm %6.1f ksi %6.1f ksi\n’,... 
slev(i), Sthoop(i), SFHoop(i)); 
end 
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Appendix I. MATLAB Code: Rotor Losses from Winding Time and 
Space Harmonics 


% Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pmharmloss 

96 Program performs rotor loss calculations caused by 

% winding time and space harmonics for permanent magnet 
% machines with surface magnets and slotted stators. 


% MUST RUN pmlinput, pm1calc, pmloutput, and get harmonic 
% current data from PSIM prior to running pmharmloss 


% Constants to be used 
mu0 = 4*pi*1e-7; % Free space permeability 
tol 2 1e-2; % Tolerance factor 


% Retaining sleeve/magnet material resistivity (ohm-m) 
Stain res = 0.72e-6; 

Titan res = 0.78e-6; 

CarFib res = 9.25e-6; 

Inconel res = 0.98e-6; 

Magnet res = 1.43e-6; 


% Retaining sleeve thickness set at 0.5mm less than air gap 
h_sl = g - 0.0005; % Sleeve thickness 
g_act = g- h sl; % Actual air gap 


% Retaining sleeve conductivities (S/m) 
cond, s = 1/Stain res; 

cond t = 1/Titan res; 

cond, c = 1/CarFib res; 

cond. i = 1/Inconel res; 


% Magnet & actual sleeve cond (S/m) 
cond m = 1/Magnet res; 
cond, sl = cond s; 


% Input time harmonic peak currents from PSIM 


I1 = 2895; 
B = 0; 

I5 = 209; 
T7 = 89.2; 
I9 = 0; 

I11 = 39.2; 
113 = 27.6; 
115 =0; 
117 = 17.0; 
119 = 12.8; 
Di =0; 
123 = 7.3; 
125 = 6.3; 
127 =0; 
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129 = 4.8; 
131 = 3.9; 


% Put currents in array 
Iharm = [I1 I3 I5 I7 19 I 11 113 115 117 119 [21 123 125... 


127 129 131]; 


% Calculate current THD 
Iah = 0; 
for r = 2:length(Iharm) 
Iah = Iah + Iharm(r)^2; 
end 
THDi = 100*sqrt(Iah/(Iharm(1)*2)); 


% Calculate current densities 

Iz = (1/sqrt(2)).*Iharm; 

Kz = ((q/2)*(Na/(2*pi*Rs))).*Iz; 

Iz 1 = (l/sqrt(2)).*I1; % Fundamental RMS current 
Kz 1 = ((q/2)*(Na/(2*pi*Rs))).*Iz 1; % Fundamental current density 


% Harmonics to be evaluated 


n=1:2:31; ` 

wzn.*omega; % Harmonic angular frequencies 
freq = w ./ (2*pi); % Harmonic frequencies 

lam = (2*(2*pi/(2*p)))./n; 

k = (2*pi)/lam; % Wavenumbers 


'% Eta values 
eta m = sqrt((j*mu0*cond_m).*w + (k.^2)); 
eta s = sqrt((j*mu0*cond_sl).*w + (k.^2)); 


% Surface coefficient at top of magnet layer 
alpha m = j.*(k/eta m).*coth(eta m.*hm); 


% Surface coefficient at top of retaining sleeve 

top1 = (j.*(k./eta s).*sinh(eta s.*h sl) + (alpha m.*cosh(eta s.*h sl); 
bot1 = (j.*(k./eta s).*cosh(eta s.*h sl) + (alpha m.*sinh(eta s.*h sl); 
alpha s = j.*(k./eta_s).*(top1./bot1); 


% Surface coefficient at surface of stator 

top2 = (j.*sinh(k.*g act)) + (alpha s.*cosh(k.*g act)); 
bot2 = (j.*cosh(k.*g act)) + (alpha s.*sinh(k.*g act)); 
alpha f = j.*(top2./bot2); 


% Surface impedance 
Zs = (mu0.*w./k).*alpha f; 


% Calculate losses due to time harmonics 

'% Use only fundamental space harmonic factors 

KZ t- kw.*Kz; 

Syt = 0; 

for i = 1:length(n) 
Sy_t(i) = 0.5*(abs(Kz t(1))^2)*real(Zs(1)); 
Syt = Syt + Sy_t(i); 

end 
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% Calculate losses due to space harmonics 
% Use only fundamental time harmonic current 
kpn = sin(n .* pi/2) .* sin(n .* alfa/2); 
kbn = sin(n .* m*gama/2) ./ (m.*sin(n .* gama/2)); 
kwn = kpn .* kbn; 
KZ s-kwn.*Kz 1 /n; 
Sysz 0; 
for i= 1:length(n) 
Sy_s(i) = 0.5*(abs(Kz_s(i))*2)*real(Zs(i)); 
Sys = Sys + Sy_s(i); 
end 


fprintf(\nRotor Losses Caused by Harmonics: n"; 
fprintf("Time Harmonic Losses = %6.1f kW\n',Syt/1000); 
fprintf('Space Harmonic Losses = %6.1f kW\n',Sys/1000); 
fprintf('Total Losses = %6.1f kW\n',(Syt+Sys)/1000); 
fprintf('Current THD = %6.2f %%\n',THDi); 
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Appendix J. MATLAB Code: Rotor Losses from Slot Effects 


96 Jonathan Rucker, MIT Thesis 

% May 2005 

% Program: pmcanloss 

% Program calculates and outputs rotor losses caused by 
% stator slot effects for permanent magnet machines 

% with surface magnets and slotted stators. 


% MUST RUN pmlinput, pmicalc, pmwave, and pmcanstress 
% PRIOR TO RUNNING pmcanloss 


Fo o Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Yo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo To % To To Vo Vo 


% Calculate retaining sleeve losses 
%%% % Fo Yo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Vo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Fo Yo Fo Fo Fo Fo Fo Fo Fo Vo Yo Vo Fo Fo Fo Fo Fo Fo Fo Fo Fo 5906907096 


% Retaining sleeve/magnet material resistivity (ohm-m) 
Stain_res = 0.72e-6; 

Titan_res = 0.78e-6; 

CarFib_res = 9.25e-6; 

Inconel_res = 0.98e-6; 

Magnet_res = 1.43e-6; 


% Calculate Bd as function of wst and wt (max 10% of Bg) 
Bd = (wst/wt)*0.1*Bg; 
% Calculate flux variation parameters 
beta = (wst/(2*pi*Rs))*2*pi; 
lamB = 2*pi/Ns; 
B = (Bd/sqrt(2))*sqrt(beta/lamB); 
% Calculate geometry and can loss factor for different rings 
% k is number of rings 
for k= 1:10 
A(k) = pi*2*(R+hm)*Lst/k; 
Ks(k) = 1 - ((tanh(p*Lst/(k*2*(R+hm))))/(p*Lst/(k*2*(R+hm)))); 
end 


% Input Stainless Steel sleeve thickness based on stress results 
for i= 1:stop 
if SFHoop(i) <= Stain_str 
t_Stain = t(i); 
break 
elseif t(stop) > Stain_str 
fprintf( Hoop Stress too high for Stainless Steel.\n‘); 
else 
dummy = t(i); 
end 
end 


% Input Titanium sleeve thickness based on stress results 
for i= l:stop 
if SFHoop(1) <= Titan str 
t Titan = t(1); 
break 
elseif t(stop) » Titan str 


171 


AC Generator Top Magnetic Generator Magnet Generator DIY 


fprintf('Hoop Stress too high for Titanium.\n’); 
else 
dummy = t(1); 
end 
end 


% Input Carbon Fiber sleeve thickness based on stress results 
fori- 1:stop l 
if SFHoop(i) <= CarFib_str 
t_CarFib = t(i); 
break 
elseif t(stop) > CarFib_str 
fprintf('Hoop Stress too high for Carbon Fiber An'); 
else 
dummy = t(i); 
end 
end 


% Input Inconel sleeve thickness based on stress results 
for i = 1:stop 
if SFHoop(i) <= Inconel str 
t Inconel = t(i); 
break 
elseif t(stop) > Inconel str 
fprintf(Hoop Stress too high for Inconel.An'); 
else 
dummy = t(i); 
end 
end 


% Calculate can losses 

w. Stain = (pi^2/3600)*((B*rpm*(R--hm))^2*t, Stain)/Stain, res; 

w. Titan = (pi^2/3600)*((B*rpm*(R--hm))^2*t Titan)/Titan res; 
w_CarFib = (pi^2/3600)*((B*rpm*(R--hm))^2*t CarFib)/CarFib res; 
w. Inconel = (pi^2/3600)*((B*rpm*(R--hm))^2*t Inconel)/Inconel res; 


for k = 1:10 
P_Stain(k) = k*w_Stain*Ks(k)*A(k)/1000; 
P_Titan(k) = k*w_Titan*Ks(k)*A(k)/1000; 
P_CarFib(k) = k*w_CarFib*Ks(k)*A(k)/1000; 
P_Inconel(k) = k*w_Inconel*Ks(k)*A(k)/1000; 
end 


% Calculate magnet losses (only with carbon steel) 

% Calculate geometry and can loss factor 

Am = pi*2*R*Lst; 

Ksm = 1 - ((tanh(p*Lst/(2*R)))/(p*Lst/(2*R))); 

% Calculate magnet losses 

% Assumes only 10% of magnet thickness sees effects 
w. Magnet = (pi^2/3600)*((B*rpm*R)^2*0.1*hm)/Magnet res; 
P. Magnet = w Magnet*Ksm*Am/1000; 

% Output Results 

z=[15 10]; 

fprintf('Retaining Can Losses:n"); 

for iz1:3 
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k = z(i); 

fprintf('%1.0f Rings:\n',k) 

fprintf('Material Thickness Can Loss\n'); 

fprintf('Stainless Steel 965.2£ mm %6.1f kW\n',t_Stain*1000,P_Stain(k)); 

fprintf(' Titanium %52fmm %6.1f kW\n',t_Titan*1000,P_Titan(k)); 

fprintf'Carbon Fiber %5.2fmm %6.1f kW\n',t_CarFib*1000,P_CarFib(k)); 

fprintf ^ Associated Magnet Loss %6.1f kW\n',P_Magnet); 

fprintf(‘Inconel %5.2fmm %6.1f kW\n\n',t_Inconel*1000,P_Inconel(k)); 
end 
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Appendix L. Results for Power Conversion Module Variants 


— 1 


o a L Conversion Module Component Calculations — p 
Rectfier Converter Input Filter Converter Output Filter 


Max 
i an Voltag 
! 


[ 6 [3255 | 76137 | ee RON ME MAE Ra mid tees sr ler =l. 
Stainless or EK EXT E reas S tate T 50% [aose MEC [reese eter | 2 Ses 
Steel a ser letes) sos | cose [saree] sose C tale | 3500 | 3508 | 432 (E sae [3.023 | 
[ 12 | 2263 | 52934 | 246 | 240.3 | 3,848.6 | 1.95E-09 | 6.66E-04 | 5.126+02| 1415 | 3.668 | 36.68 | 422 [1626-06 5327 | 3.043 | 
[ s | 3256 | 7,616.1 | 1019 | 2322 |2,674.8]|2.80E-03| 1.256-03 |3.56E«02| 1.415 | 2.547 | 25.47 | 507 [t.62E-06| 5575 | 2.908 | 
Titanium |—2—1 2260 1 64559 | 549 | 2423 |3,1555 | 2.586-03 | 1066-03 [4.206+02| 1.415 | 3006 | 3006 | 466 ]i62E-oe| 5473 | 2.962 | 
[40 pease [57405] 348 | 2045 [2522 L 00 fsted faren] 115 I aore I 596 I 440 MOE soot I 2007 | 
3 Phase [ 12 | 2352 | 55015 | | 55015] 266 | 2225 Pars} ame tte mec te set ais | TI RR | 430 [162E-oe| 535.8 [ 3.026 | 
[ 6 | 3380 | 79061 J 1099 | 2155 [2,576.7 | 2.91E-03 | 129E-03|3.43E«02| 1.415 | 2.453 | 2453 | 516 [1.62E-06| 559.6 | 2.897 
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| i2 | 2380 | 55670 | 272 | 2173 | 3,659.4 | 2.05E-03 | EXISCHEELSS RETE REC ster] aas heao CUN 
[ 5 | 3182 | 74430 | 974 | 2431 |2,737.1| 2.74E-03 | 1226-03 [3.64Es02| 1.415 | 2.607 | 2607 | 5.01 [1.62E-06[ 5562 | 2.915 | 
YY IY@TYC ETE EST BESSER eI ee iR] [5064 | 462 [1.6206] 546.1 | 2.969 | 
| 2,405 | 5,6255 | 334 | Pass ELSE FREE Ed a BE EE h SUR ser s] 

me EG ES e MESSE ER KE EE SE ME RE 


IE EEE LOO EE EE 


— eee LF eens sete rer | ses Dee pm [uetus 
Stainless | 8 | 1,564 | HE as es eE tae Laser asss Te fer a 
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[—310 | 1.452 | 3841.9 | 122 | 7005 [5.3025 | 1.41£-03 | 6206-04 |7.066+02| 1.415 | 5.056 | 5056 | 3.60 [.62E-06| 502.0 | $230 | 
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| Carbon [esse TN SE fais SEN NE SS MK 
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| 8 120 essa ara I seis forst fase E ta reset tds EE EE I 470 H 0 5483 | 2006 | 
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a a = S eae EE FTSE 4 IE ETE EE RETO UC BEE 
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Appendix M. Results for Power Module Losses 


Retaining 


Phase: Can 


Core 
[6 [| 20 | 
[| 10 ] i59 | 


Titanium 


m 
N 
i^] 


[12 | 
F s | 
aaa 
[12 | 
3 Phase [—6 —] 
[ 9 | 
[10 | Ë 

[ 2 ] 

= | 

[ 9 | 

[. 10 | 


Core 
Losses 


Retalning 


Can Potes 


Stainless 
Steel 


Titanlum 


Retalning 


° 
3 


| 165 | 
| 356 | 
[ 119 | 
[ 132 | 
[16.0 | 
| 180 | 
Core 


| 6 | 
Stalntess[ 8 | 
Steel | 10 | 


11.9 

[ 116 | 
[ 13.3 | 
| 15.0 | 
[ 2 J| 174 | 


| 139 | 

Titanium [160 ] 

| 144 | 

[ 165 | 

[| 13.4 | 
17.9 


Retaining 
Can 


Core 
Losses 


Stainies 
Steel 


Titanium 


9 Phase 


m 
N 
N 


Core 
Losses 


Retaining 


11 Phase 


AC Generator 


Bl 
Dp 
3 


Stainless 
Steel 


Titanium 


Power Conversion Module _ _ 
Rotor Losses Machine IGBT E 
Windage Slot Space | Winding | Winding | Total | switching 
Losses Tim Space 
Ham Harm | Harm | Losses 
[ 323 [| 1155 | 661 [| 193 | 217 [ 2669 | 2,253 | 1993 | 1765 | 3,201. | 
[ 298 [| 1396 | 741 j 399 | 922 | 3915 | 2,073.0 | 2717 167.7 | 25124 P 


579 | 3815 [920083 | 2805 | 1666 | 
28213 | 199.5 | 1765 | 31972 


Armature 
Losses 


| 


š 
3 Å 
a 
L1 
o 
9 


i^] 
m 


[ 604 | 3664 | 20774 | 2711 | 1678 [| 2,5163 | 

[ 221 | 250 | 2417 | 2,923.9 | 1925 | 1773 | 3,293: | 
[ 2840 | 24189 | 2326 | 1724 | 2,824.0 | 

[ 365.7 | 2.1840 | 2579 | 1694 | 2,6112 | 

.1 | 2,1036 | 2676 | 1682 | 25394 | 


| 255 | 390 | 425 [ 2882 | 2,509 | 2392 | 1716 | 27618 
[ 250 | 499 [| 1435 [ 4051 J 21053 | 2673 | 1682 | 25408 
e s o8 | 755 | 3585 TI 20380 | 2763 | 1671 | 25914 I 


Armature Winding | Winding 
Losses | SIL Spe Time | S 
larm Harm L. 
2871 | 15053 | 3735 | 1553 [| 20341 ES, 
579 | 3467 | 12245 | 4596 | 1448 | 18289 (2 


[20309 


N 


Losses 


31.0 


e| 
N 
œ 


[^l he MIR) 
3 Ë GIS 
N ES [7] 


co} cohen} on} a} EN 
RKS SIN 2 le 
ofi Njoo] 


NI 
a a 
MIA 


515.1 | 


Time 

Harm arm, 
| 575 [ 298 | 413 | 2918 ] 21172 | 2658 | 1684 | 25514 [2 
| 573 | 491 | 705 [ 3570 T 1791.8 | 3140 | 1625 J 252683 | 

1404 | 602 [| 676 | sos [ 4011 T 1,301 | 3451 | 1587 | 21339 E 2,5350 
| 1367 | $639 | 826 | 1253 | 4559 [1,564.7 | 3598 | 1570 | 2,0815 | 
| 1132 [| 414 | 319 [| 441 | 2727 | 21388 | 2630 | 1687 | 2,5705 F 2,8332 | 
[ 1305 [ 412 | soo J 729 | 3388 | 1,607 | 3024 | 1639 J 2,3271 |. 


| 465 | 856 | 1290 | 4489 | 1,6200 | 3475 | 1585 | 21259 | 
[ 273 [| 320 | 461 | 2524 [| 22419 | 2509 | 1702 | 2,663.0 | ES 15399] 
| 277 | 379 | 593 | 2825 [| 21462 | 2622 | 1688 | 2,5772 T 2959.7 | 


| 406.9 | 1,666.1 | 3377 | 1596 | 2,1635 | 
[ 1689 | 473.4 | 1.59048 | 3529 | 1578 | 2105.5 | 
Diode 
Cond 
Losses 


121.4 


[3871957] 
Pa 
PST 


29.4 


32.3 
32.1 


137.0 
101.3 


rm mr 
SISI] 
H NI |b] 


Ar 
Losses 


Winding 
Time 


Windage Winding 
Space 


150.0 
139.9 


148.6 | 3,890.6 | 2,18660. 


52 964. 


: 
[24703 1 


[5623 | 132.3 | 1.6947 | [ 25 8257 


179 


Top Magnetic Generator Magnet Generator DIY 


Page Intentionally Left Blank 


180 


AC Generator Top Magnetic Generator Magnet Generator DIY 


Appendix N. Results for Power Module Weights 
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Appendix O. — Rectifier/Input Filter Mass and Volume Calculations 
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Appendix P. ^ Converter/Output Filter Mass and Volume Calculations 
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[ 0.071 | 22,857 | 5562 | 145.3 | 

[ 2.969 | 073 | 0.073 [22,57| 5461 | 1426 [ 427.9 | 

[ 3016 | 074 | 0.074 | 2257] 5375 | 1404 | 4212 | 12 | 

[12 1 7000 | 3.035 | 074 | 0.074 | 22.857] 5341 | 1395 | 418.6 | 
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[ 6 j 7000 | 3.061 | 075 | oo75 |22857| 5296 f 138.4 | 4154 | 12 | 
Stainless| 8 ^ | 7000 | 3.172 | 0.78 | 0.078 | 22,857] 5100 | 1332 | 3997 | 12 | 
Steel | 10 ^[ 7000 | 3269 | 080 [ 0080 [22857] 4959 | 1295 | 3886 | 12 | 
[ 12 ]7000]| 3318 | 081 | 0.081 [22857| 4887 | 1276 | 3829 | 12 | 

7000 | 3059 | 0.75 [ 0075 [22857| 5299 | 1384 | 4153 | 12 | 

7000 | 3155 | 0.77 | 0077 [22,857[ 5139 ] 1342 | [ 12 | 
I [ 0.79 | 0079 | 22857] 5001 | 1307 | | 12 | 
7000 | 3294 | 081 f 0.081 L22857] 4922 | 128.6 | 3857 | 12 | 
J [ 12 | 

[ 12 | 


3997 | 0.078 | 228 | 0.45 | 4230 
0.080 | 228 | 045 | 4120 | 
[ 0.081 | 228 | 0.45 | 4063 | 
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o j 3038 | 074 | 0.074 | 22857 | 5336 | 1304 | 

7000 | 3.141 | 0.77 | 0.077 | 22,857 | 516.1 | 1348 | 404. 

[3230 | 075 | 0079 |-22,857 | 502.0 | 131.1 | 393.4 | 12 | 

7000 | 3279 | 0.50 | 0.080 | 22,857 | 4945 | 1292 | 5875 | 12 | 

7009 | 3054 | 0.75 | 0075 22,857 530.8 | 138.7 | 4160 | 12 | oas | —19 | 4169 | 0.075 | 228 | 0.45 [430.3 | 

7000 | 3.169 | 078 | 0078 | 22,857| 51.6 | 1936] 4009 | 12 | 045 | 19 | 400. | 0.078 | 228 | 045 | 4242 | 

7000 | 3239 | 079 | 0079 | ?2,857| 5005 | 1307 | 3922 | 12 | 045 | 19 | 022 0079 | 228 | 04 | 4155] 
700.0 | 3.203 | [0.081 -| 22,857 | 4923 | 128.6 | 5656 | +12 | 045 | 19 [5858 | 0.081 | 226 | 045 | 409.1] 


Output Average Total Total Total 
Kena Capacitance | Capacitor] Capacitor | Output | tnductance | Inductor | Inductor [IGBT Mass 
a Energy (J| Masa (kg) | Current | (uH) — |Energy (J| Mass (ko) 


[6 [700.0 | 2906 | 071 | 0071 J 22,857 [557.8 | i457 | 4371 | 12 [045 | 19 14374 | 0.071 | 228 [| 0.45 | 4605 | 
| 8 | 7000 | 2961 | 073 | 0.073 | 22,857| 5475 | 1430 | 4291 | 12 | 045 | 19 T 429.1 | o073 | 228 [ 0.45 | 4524 | 

[ 30 | 7000 | 3.002 | 074 | 0.074 [22,57 | 5401 | 1411 | 4233 | 12 | 045 | 19  [ 4233 | 0.074 | 228 | 0.45 | 446. | 

[ 12 | 7000 | 3.020 | 074 | 0.074 | 22857] 5353 | 1398 | 4195 | 12 J 0.5 | 19 | 4195 | 0.074 | 228 | 0,45 | 4428 | 

| 6 | 7000 | 2904 | 071 | 0071 [22857| 5582 | 1458 | 4375 | 12 | 045 | 19 4375 | 0071 | 228 | 0.45 $ 4608 | 

Titanium — 8 | 7000 | 2:947 | 072 [ 0072 | 22857] 5504 | 1437 | 4911 | 12 ] 045 | 19 f 431.1 L 0072 | 228 | 045 | 4544 | 
| 30 | 7000 | 2984 | 073 | 0.073 | 22857] 5433 | 1419 | 4258 | 12 | 045 | 19 | 4258 | 0073 | 228 [| 0.45 | 449.1 | 

[7000 | 3010 | 0.74 | 0.074 | 22,857] 5385 | 1407 | 4220 | 12 f 0.45 | 19 | 4220 | 0.074 | 228 J 0.45 | 445.4 | 

[ 2:591 | 0:71 [| 0.071 |22,57| 5608 | 1465 | 4394 | 12 | 0.45 | 19 f 4394 | 0071 | 228 | 0,45 | 4628 | 

| 2939 T 072 | 0.072 | 22,857[ 5516 | 1441 | 4323 | 12 | 045 | 19 | 4323 | 0.072 | 228 | 045 | 4559 | 

[ 2078 | 073 | 0.073 |22857| 5444 | 1422 | 4266 | 12 f 045 | 19 | 4266 | 0073 | 228 | 0.45 | 4499 | 

[ 3.002 [ 074 [ 0074 |22,857]| 5400 | 1411 | 4232 | 12 f 045 | 19 | 4232 | 0.074 | 228 | 0.45 | 4465 | 

| 2.902 [| 071 | 0071 | 22,857 | 5586 | 1459 [ 4378 | 12 f 045 | 19 | 4378 | 0.071 | 228 ] 045 | 4611] 
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[ 2.956 | 0.72 | 0.072 | 22,8577 5483 | 1432 | 4297 | 12 f 045 | 19 ] 4297 | 0072 | 228 | 045 | 453.0 | 
0.073 | 22,857 | "5418 | 1415 | 4246 | 12 jJ 0.45 | 19 4246 | 0.073 | 228 | 045 | 4479 
[3019 | 074 [ 0.074 [22857] 5370 | 1403 | 4208 | 12 ] 045 | 19 420.8 0074 | 228 | 045 [4442 | 
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0.073 | 22,857 | 5432 | 1419 | 4257 | 12 j 0.45 | 19 [425.7 | 0.073 | 228 f 0.45 | 
[ 138.4 | å 7 


[ 3.054 | 
Titanlum 
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[ 0.073 | 
| 39 [| 4151 | 0.075 | 
| 39 [4045 [ 0077 | 
| 19 [| 3965 | 0.079 | 
[ 19 14257 J 0075 | 
[ 319 [| 4160 | 0.075 | 
[ 19 [408.0 | 0.076 | 
[ 19 | 4001 | 0.078 | 
[ 9 [| 4169 | 0.075 | 
| 39 [24095 | 0.076 | 
[ 19 [ 4015 | 0.078 | 
| 19 4246 | 0.073 | 
[ 39 | 4143 | 0.075 | 
[ 19 | 406.4 | 0.077 | 
[ 19 ] 398.6 | 0.078 | 


m all 
M ho (of ho] 
bad bad bad be Ë 
ajajaja 

len 


3.087 0.076 | 22,857 
Stainless 3.191 22,857 
Steel .o | 3299 [| 081 [| 0.081 | 22857 E 


3.068 


[ 3.068 | 
| 3.182 | 
| 3281 | 
| 3.368 | 
| 3.062 | 075 | 

| 3164 | 078 [ 0.078 | 22,857} 5124 | 1338 | 
| 3.262 | o. 

[ 3367 | 

[ 3.076 | 

[ 3.173 | 

[ 3.280 | 
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| 0.80 | 0.080 | 22857] 4969 J 129.8 | 
3.367 | 082 | 0.082 [22,57] 4814 | 125.8 | 
3.076 
3.173 | 078 | 0.078 | 22,857 | 5109 | 1335 | 
0.83 | 0.083 | 22,857 | 4812 | 
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6.0106 | 00156 | 19  [ 179 | 0.0001 | 0.20 J 028 | 0.07 |. 2.92..| 
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19.0001 | 020 ] 026 | 067 |.290 | 
| 173 [00001] 020 | 026 | 068 |.2.85 | 
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